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Introduction 


The  goals  of  this  research  have  been  to  develop  in  situ  optical  instruments  and  computer 
models  of  the  apparent  and  inherent  optical  properties  of  the  submarine  light  field.  The 
radiative  transfer  function  in  seawater  is  a  complicated  mix  of  signals  caused  by  many 
components  some  of  which  are  covariant.  The  approach  taken  in  this  research  has  been  to 
measure  by  utilizing  custom  and  commercial  optical  instrumentation  as  many  of  the  components 
of  light  field  at  a  single  site  rather  than  combine  spotty  information  at  many  sites.  A 
fundamental  contribution  to  the  program  was  our  participation  in  the  Marine  Light  Mixed  Layer 
Experiment  in  the  North  Atlantic.  This  provided  an  opportunity  to  simultaneous  measure  the 
optical  components  of  the  light  field  in  a  multi-investigator  environment.  We  supported  this 
project  with  state-of-the-art  high  spectral  resolution  optical  field  instrumentation  and  computer 
modelling  of  the  submarine  light  field.  The  apparent  optical  properties  of  the  MLML  site  were 
measured  with  the  Spectral  Transmissometer  and  Radiometer  (STAR)  and  the  Average  Cosine 
Meter  (ACM). 


RESULTS 

There  were  three  major  efforts  this  program.  These  included  instrument  development, 
modelling,  and  cruise  support.  The  development  of  optical  instruments  and  sampling  technique 
are  synergistic  with  the  computer  modelling.  Advances  in  instrumentation  provided  more  known 
quantities  in  the  modelling  effort  which  indicated  the  direction  to  pursue  new  instrumentation. 


Instrument  efforts 

There  were  effectively  two  instruments  developed  under  this  project.  The  proposed 
instrument  was  the  Average  Cosine  Meter.  The  prototype  was  completed  using  what  was  at  the 
time  state-of-art  components.  It  was  deployed  on  a  small  floating  reel  and  electromechanical 
release,  this  novel  and  inexpensive  method  eliminated  shadowing  from  the  ship.  The  ACM 
measured  irradiance  at  490nm  using  cosine  and  411  collectors  as  well  as  detector  dark  current. 
The  ACM  was  utilized  in  conjunction  with  the  Spectral  Transmissometer  and  Radiometer  which 
provided  hyperspectral  downwelling  irradiance  and  upwelling  radiance.  During  the  summer 
experiment  at  MLML  site,  the  STAR  was  lost  when  the  hydrowire  on  the  USNS  Bartlett  failed 
due  to  age  and  fatigue.  A  new  hyperspectral  instrument  was  constructed  from  plans  and  spare 
parts  that  were  already  on  hand  to  upgrade  the  STAR.  The  new  instrument  due  to  its  reduced 
size  and  increased  capabilities  effectively  replace  both  the  ACM  and  the  STAR,  This  instrument 
has  been  named  the  Submersible  Upwelling  and  Downwelling  Spectrometer  SUDS.  It  is  a 
hyperspectral  upwelling  radiometer  and  downwelling  irradiometer  utilizing  fiber  optics,  and 
shutter  to  a  common  holographic  grating  spectrometer  and  512  element  linear  array  detector. 
(Steward,  et  al.  1994) 

Modelling  efforts 

During  this  project,  a  Monte  Carlo  type  model  of  the  submarine  spectral  irradiance  field  was 


produced  (Reinersman  et  al.  1994).  It  uses  measured  or  simulated  incoming  solar  and  sky 
irradiance  and  the  inherent  optical  properties  of  the  sea.  It  has  been  used  to  help  validate  Q- 
factors  for  a  remote  sensing  reflectance  model  (Lee  et  al.  1994)that  inverts  high  spectral 
resolution,  remote  sensing  reflectance  or  irradiance  reflectance.  These  reflectance  data  can  be 
interpreted  in  terms  of  absorption  that  when  divided  by  the  average  cosine  provides  a  good 
estimate  of  the  diffuse  attenuation  coefficient.  Specifically,  one  can  predict  the  spectral  vector 
irradiance  and  spectral  scalar  irradiance  fields  with  depth  for  variable  chlorophyll-a 
environments.  These  fields  specify  the  general  availability  of  light  for  passive  visual  or  other 
optical  observations  and  the  spectral  quanta  available  for  absorption  by  phytoplankton  for 
photosynthesis. 


Cruise  Support 

Four  persons  participated  in  two  cruises  to  the  MLML  site  in  the  North  Atlantic. 

The  first  cruise  was  aboard  the  R/V  Endeavor  from  May  14-June  2,  1991  (Reykjavik,  Iceland 
to  Narragansett,  Rhode  Island).  This  was  a  process-oriented  experiment  with  a  number  of 
investigators  working  on  concurrent  sampling.  Some  of  the  data  collected  on  this  experiment 
by  our  group  included: 

a.  Specific  absorption 

b.  Remote  sensing  reflection 

c.  CDOM  absorption 

d.  CDOM  absorption  and  fluorescence  on  concentrated  extractions 


The  second  cruise  was  the  Mesoscale  Mapping  experiment  aboard  the  USNS  Bartlett  August 
11-31,  1991  Cork,  Ireland  to  Newport,  l^ode  Island.  Due  to  restriction  in  access,  rough 
weather,  and  the  transect  nature  of  this  work;  fewer  samples  were  collected.  A  similar  suite  of 
samples  were  collected  to  make  seasonal  comparisons  with  the  spring  processes  cruise. 

RESULTS/IMPACT 

The  significance  of  the  work  performed  under  this  funding  is  best  described  by  summary  of  the 
peer-reviewed  publications  which  have  resulted.  These  publications  have  been  included  in  the 
Appendix. 

1.  Aircraft  flying  below  the  clouds  (e.g.  Unmanned  Air  Vessels  UAV’s)  will  be  able  to 
derive  accurate  absorption  coefficients  for  phytoplankton  and  gelbstoff  and  accurate 
chlorophyll  a  concentrations  from  hyperspectr^  remote  sensing  reflectance  data. 
(Carder,  et  al.  1992) 

2.  The  data  derived  above  can  be  used  to  estimate  primary  productivity  values  some  three 
times  more  accurately  than  using  traditional  methods.  A  mathematical  simulation  of  the 
absorption  coefficient  of  phytoplankton  derived  from  hyperspectral  remote  sensing 
reflectance  makes  more  accurate  estimation  of  the  primary  productivity  than  techniques 
that  use  chlorophyll  a  concentration,  irradiance,  average  chlorophyll  specific  absorption. 


and  quantum  yield  of  carbon  fixation.  (Carder  et  al.  accepted,  Lee,  et  al.  submitted) 


3.  Measurements  of  dissolved  organic  carbon  DOC  on  transects  to  and  from  the  MLML  site 
contributed  to  a  model  of  flux  of  dissolved  organic  matter  DOM  and  dissolved  nitrogen 
between  the  North  Atlantic  and  the  equator.  Results  of  the  modelling  suggest  how 
future  optical  measurements  by  remote  sensing  platforms  will  be  important  in 
discrimination  of  particulate  and  dissolved  components  of  the  organic  carbon  flux  and 
verify  the  importance  of  continental  margins  and  DOM  in  global  biogeochemical  cycles. 
(Walsh,  Carder,  and  Muller-Karger,  1992) 

RELATIONSHIP  TO  OTHER  PROJECTS 

1.  This  project  provided  high  latitude  data  helpful  in  validating  Seawifs  algorithms. 

2.  This  project  developed  methods  of  deriving  quantum  yield  data  for  primary  production 
useful  to  John  Marra  and  other  coinvestigators  that  deployed  moorings  and  drifting  arrays 
equipped  with  primary  productivity  incubators  and  simple  photosynthetically  available 
radiation  PAR  meters. 


PUBLICATIONS 

P  -  Lee,  Z.  P.,  K.L.  Carder,  S.K.  Hawes,  R.G.  Steward,  T.G.  Peacock,  and  C.O.  Davis. 
1994a.  A  model  for  interpretation  of  hyperspectral  remote  sensing  reflectance.  Applied  Optics. 
5721-5732. 

PI  -  Costello,  D.K.,  K.L.  Carder,  and  W.  Hou.  (submitted  to  a  special  issue  of  Deep-Sea 
Research).  Aggregation  of  a  diatom  bloom  in  a  mesocosm:  bulk  and  individual-particle  optical 
measurements. 


PI  -  Costello,  D.K.,  K.L.  Carder  and  R.G.  Steward,  (in  prep,  for  a  submission  to  Deep-Sea 
Research).  Design  and  performance  of  the  Marine  Aggregated  Particle  Profiling  and 
Enumerating  Rover  (MAPPER). 

PI  -  Costello,  D.K.  and  K.L.  Carder,  (in  prep,  for  submission  to  Journal  of  Geophysical 
Research).  Underwater  light  attenuation  by  the  3-dimensional  distribution  of  large  marine 
particles  in  Monterey  Bay. 

R  -  Costello,  D.K.,  W.  Hou  and  K.L.  Carder.  1994.  Some  effects  of  the  sensitivity  threshold 
and  spatial  resolution  of  a  particle  imaging  system  on  the  shape  of  the  measured  particle  size 
distribution.  In:  Ocean  Optics  XII.  Jules  S.  Jaffe,  Editor,  Proc.  SPIE  2258,  Bergen,  Norway. 

R  -  Costello,  D.K.  and  K.L.  Carder.  1994a.  New  instrumentation  and  platforms  for  subsurface 
optical  measurements.  In:  Ocean  Optics  XII.  Jules  S.  Jaffe,  Editor,  Proc.  SPIE  2258,  Bergen, 


Norway. 


R  -  Costello,  D.K.  and  K.L.  Carder.  1994b.  Using  unmanned  vehicle  systems  for  ground- 
truthing  oceanographic  satellite  data.  AUVS  ’94  -  Intelligent  Vehicle  Systems  and  Technology. 
Proc.  Asso.  for  Unmanned  Vehicle  Systems,  Detroit,  MI. 

R  -  Costello,  D.K.,  K.L.  Carder,  and  R.G.  Steward.  1994.  Particle  size  spectra  and  optical 
properties.  In;  Mesocosm  Experiment  Data  Report  to  the  Office  of  Naval  Research  under  the 
Accelerated  Research  Initiative,  Significant  Interactions  Governing  Marine  Aggregation 
(SIGMA).  51  pp. 

R  -  Costello,  D.K.,  K.L.  Carder,  W.  Hou  and  R.G.  Steward.  1994.  3-dimensional  particle  size 
spectra,  transmissometry  and  temperature  profiles,  and  filter  pad  absorption  data.  In:  Monterey 
Bay  Field  Experiment  Data  Report  to  the  Office  of  Naval  Research  under  the  Accelerated 
Research  Initiative,  Significant  Interactions  Governing  Marine  Aggregation  (SIGMA).  24  pp. 

ONR  SPONSORED  PRESENTATIONS 

C  -  Costello,  D.K.,  K.L.  Carder,  R.G.  Steward  and  A.  Alldredge.  1994.  The  Distribution  and 
Optical  Properties  of  Large  Marine  Particles:  Data  From  a  Culture  Tank  and  Field 
Experiments.  EOS  AGU/ASLO  (75),  3,  35. 

C  -  Hou,  W.,  D.K.  Costello,  K.L.  Carder  and  R.G.  Steward.  1994.  High-resolution  Marine 
Particle  data  from  MAPPER,  a  new,  in  situ  optical  ocean  instrument.  EOS  AGU/ASLO  (75), 
3,  21. 


C  -  Lee,  Z.  P.,  K.  L.  Carder,  and  T.  G.  Peacock,  1994b.  Hyperspectral  modeling  of  remote 
sensing  reflectance:  from  the  Florida  Shelf  to  the  Mississippi  River.  EOS  AGU/ASLO  (75)  3, 
193. 

C  -  Peacock,  T.  G.,  K.  L.  Carder,  P.  G.  Coble,  Z.  P.  Lee,  and  S.  K.  Hawes,  1994.  Long-path 
Spectrometer  for  measuring  gelbstoff  absorption  in  clear  waters.  EOS  AGU/ASLO  (75)  3,  22. 

C  -  Reinersman,  P.,  F.  Chen,  and  K.  L.  Carder,  1994.  Monte  Carlo  Simulation  of  Atmospheric 
Point  Spread  Function  With  Application  to  Coastal  Image  Processing.  EOS  AGU/ASLO  (75) 
3,  136. 

C  -  Steward,  R.  G.,  K.  L.  Carder,  and  T.  G.  Peacock.  1994.  High  resolution,  in  water  optical 
spectrometry  using  the  Submersible  Upwelling  and  Downwelling  Spectrometer  (SUDS).  EOS 
AGU/ASLO  (75)  3,  102. 

C  -  Young,  L.  R.,  K.  L.  Carder,  and  P.  Hallock,  1994.  Hyperspectral  remotely  sensed  colored 
dissolved  organic  matter  as  a  photoprotective  agent  against  bleaching  in  corals  and  foraminifera. 
EOS  AGU/ASLO  (75)  3,  103. 


APPENDIX 

The  following  publications  or  drafts  have  been  included; 

Carder,  K.L.,  Z.  Lee,  T.G.  Peacock,  R.G.  Steward,  C.O.  Davis,  1992.  Aircraft  algorithms  for 
water  absorption  coefficients:  rain  or  shine.  Western  AGU 

Carder,  K.L.  Z.  Lee,  J.  Marra,  R.G.  Steward,  and  M.J.  Perry,  (in  press).  Calculated  quantum 
yield  of  photosynthesis  of  phytoplankton  in  the  Marine  Light-Mix^  Layer  (59°N/21°W),  J. 
Geophys.  Res.  (Special  Edition). 

Lee,  Z.,  K.L.  Carder,  S.K.  Hawes,  R.G.  Steward,  T.G.  Peacock,  C.O.  Davis,  1994.  Model 
for  the  interpretation  of  hyperspectral  remote  sensing  reflectance.  Applied  Optics  33:24,  5721- 
5732. 

Lee,  Z.,  K.L.  Carder,  J.Marra,  R.G.  Steward,  M.J.  Perry,  (in  press).  Estimating  primary 
productivity  from  remote  sensing  reflectance;  a  new  approach.  J.  Geophys.  Res. 

Walsh,  J.J.,  K.L.  Carder,  F.E.  Muller-Karger,  1992.  Meridional  fluxes  of  dissolved  organic 
matter  in  the  North  Atlantic  Ocean.  J.  Geophys.  Res.  97-clO,  15625-15637. 


aircraft  Algorithms  for  Water  Absorption 
Coefficients;  Rain  or  Shine 

L  Carder.  Z  Lee,  T  G  Peacock,  R  G  Steward  (All 
at  Department  of  Marine  Science,  University  of 
South  Florida,  140  7th  Avenue  S.,  St.  Petersburg, 
FL  33701-5016;  813-893-9148 
:  0  Davis  (Jet  Propulsion  Laboratory,  Mail  Stop 
300-323,  4800  Oak  Grove  Drive,  Pasadena,  CA  91109) 

algorithms  for  calculating  the  absorption  and  the 
liffuse  attenuation  coefficients  of  surface  waters 
;ere  derived  for  use  in  aircraft  flying  near  the  sea 
-.urface  during  sunny  or  overcast  conditions.  Cloud- 
:ast  areas  are  typical  of  many  upwelling  or  high- 
latitude  regions  where  autonomously  piloted  vehicles 
^APVs)  may  be  needed  to  remotely  measure  water 
:larity  and  chlorophyll  pigments  to  supplement  ocean 
'.olor  satellite  data.  Data  from  forty  stations  were 
ised  to  develop  algorithms  for  the  absorption 
(a(490))  and  diffuse  attenuation  {kd(490)) 
coefficients  at  490  nm  that  are  based  upon  spectral 
ratios  of  remote -sensing  reflectance 
(Rr5(520)/Rrg(560) ) .  The  algorithms  provided 
estimates  with  correlation  coefficients  of  0.98  and 
accuracies  better  than  +  10%.  A  curvature  algorithm 
is  also  presented  that  flags  R„  data  which  contain 
effects  of  bottom  reflection.  The  diffuse 
attenuation  coefficient  is  used  to  classify  optical 
■/ater  types  and  to  parameterize  models  that 
characterize  the  vertical  structure  of  the 
chlorophyll  field  for  various  bio-optical  domains. 


1.  1992  Fall  Meeting 

2.  Member  #000575311 

3.  (a)  K  L  Carder 

Dept,  of  Marine  Science 
Univ.  of  South  Florida 
140  Seventh  Avenue  South 
St.  Petersburg,  FL  33701 

(b)  tel.  813-893-9148 

(c)  fax.  813-893-9189 

4.  0 

5.  (a)  Marine  Light  Mixed  Layer 
(special  session)  or 

008  Remote  Sensing  of  Ocean 
Color 

(b)  4275,  4855 

(c)  Ocean  color  algorithms 


7.  0% 

8.  Charge  K.L.  Carder 

AM  EX  AX  3787  020098  21  005 
Exp.  7/94 

9.  C 

10. 


11.  No 


Aircraft  Algorithms  for  Water  Absorption  Coefficients: 

Rain  or  Shine 


K.L.  Carder^  Z.  Lee^  T.G.  Peacock^  R.G.  Steward^ 

Department  of  Marine  Science,  University  of  South  Florida 
140  Seventh  Avenue  South,  St.  Petersburg,  FL  33701 
C.O.  Davis^ 

Jet  Propulsion  Laboratory,  California  Institute  of  Technology 
4800  Oak  Grove  Blvd.,  Pasadena,  CA  91109 


1.  Introduction 


In  1981,  Austin  and  Petzold^  developed  an  algorithm  to  estimate 
the  diffuse  attenuation  coefficient  K  at  490nm  using  the  in-water 
upwelling  radiance  (LJ  ratio  at  443  and  550nm: 

L  (443) 

/s:(490)=0.0883(-!^^^ - +0.022  d) 

\p50y 

The  importance  of  this  value  is  that  it  can  be  used  to  classify  optical 
water  types.  can  be  related  directly  to  the  water-leaving  radiance  (LJ 
and  remote-sensing  reflectance  see  below),  so  a  similar  expression 
may  be  written  using  spectral  ratios. 

By  applying  a  similar  algorithm  to  our  West  Florida  Shelf  (WFS) 
data  and  other  stations  (Table  1),  we  found  the  correlation  between  the 
absorption  coefficient  at  490nm  (^2(490))  and  the  remote-sensing 
reflectance  ratio  of  443  to  550nm  is  less  than  for  an  algorithm  using  the 
ratio  of  520  to  560nm  (Fig.l).  This  may  be  explained  in  part  by  the 
influence  of  bottom  reflectance  and  other  components  not  covarying 
with  chlorophyll  for  these  "Morel  Case-2"^  waters.  J^(490)  and  ^2(490) 
are  equivalent  in  near-surface  waters  (Eq.  6). 

For  studying  regions  that  include  shallow  coastal  and  "Morel  Case- 
2"  waters,  we  present  an  algorithm  for  the  calculation  of  a(490)  based 
on  the  study  of  39  stations,  including  blue-sky  and  over-cast  sky 


conditions.  A  curvature  algorithm  also  is  presented  to  "flag"  remote 
sensing  data  in  which  bottom  reflectance  represents  more  than  20%  of 
the  total.  We  also  studied  the  effect  of  the  20%  bottom  contribution  on 
the  calculation  of  pigment  (Chi  a  -f-  phaeo  a)  concentration  using  the 
"Morel  Case-l"  algorithm^. 


Table  1 


Data  source 

Date 

Number 

of 

stations 

Number  of 

overcast 

stations 

Lake  Tahoe 

Aug.  1990 

1 

0 

TBX  I 

Mar.  1990 

3 

0 

TBX  II 

May  1992 

11 

0 

MLML  I 

May  1991 

10 

6 

Monterey  Bay 

Apr.  Sep.  Oct. 

1989 

9 

2 

TTOlO 

July  1992 

5 

5 

TBX:  Tampa  Bay  experiment  (subtropical  Gulf  of  Mexico) 
MLML:  Marine  Light  Mixed  Layer  (high  lat.  North  Atlantic) 
TTQIO:  Thomas  Thompson  010  (Pacific,  San  Diego  to  Hawaii) 


2.  Background  and  Results 
A.  Absorption  Algorithm. 

Since  we  do  not  have  total  absorption  measurements  for  most  of 

I 

j 

the  stations,  we  used  the  remote-sensing  reflectance  model  of  Lee  et  at 
to  obtain  a(490).  We  call  this  modeled  value  "Mod.  ^(490)”.  Figure  1 
shows  the  relation  between  the  reflectance  ratio  and  Mod.  a(490). 
Figure  2  compares  the  reflectance  ratio  to  the  sum  of  absorption 
components  due  to  water  molecules  and  particles  (a^-hap),  which 
includes  all  in-water  absorption  except  a  small  gelbstoff  residual.  We 
find,  in  both  cases,  a  higher  correlation  when  520  to  560nm  reflectance 
ratios  are  used. 

Performing  regression  analysis  on  the  Rjs(520)/Rrs(560)  data  in 
Figure  1  yields  an  algorithm  for  calculating  a(490), 

/?  f520) 

Ca/c.a(490)=0.14(-!^ - (2) 

R^peo) 

with  r^=0.96  and  coefficient  of  variation  of  12.5%. 

In  order  to  compare  Eq.  2  to  Eq.  1,  we  note  that  remote-sensing 
reflectance  is  defined  by 


(3) 


where  is  the  above  surface  water-leaving  radiance,  and  is  the 


above  surface  downwelling  irradiance.  Methods  are  detailed  in  Lee  et 
at.  Also,  the  relation  between  and  is 


with  t  =  transmittance  of  the  air-sea  interface,  and  n  =  the  index  of 
refraction  of  seawater.  Applying  Eqs.  3  and  4,  Eq.l  can  be  rewritten 

^:(490)=0.0883(-!2:^ — +0.022  (5) 

/?„(550)*£^550)" 

For  oceanic  waters,  the  back-scattering  coefficient  b^,  is  small,  so  a  and 
K  can  be  related  by 

a(490)=^i*^(490) 

\ 

where  is  a  path-length  factor.  We  used  an  average  value^  of 
jLt«  1/1.5=0.67,  and  of  jEd(443)/Ed(550)«0.95,  so  Eq.  5  is  recast  as 

/?  (443)  -  -Q- 

Ca/c.a(490)=0.064(-^!^- - +0.015  (T) 

Figure  3  compares  the  results  using  Eq.2  and  Eq.7  to  estimate 
a(490).  If  we  apply  appropriate  corrections  to  the  reflectance  ratio  in 
Eq.2  to  simulate  Case  I  waters  (e.g.  remove  bottom  reflectance,  colored 
dissolved  organic  matter  (CDOM)  fluorescence,  and  water  Raman 


scattering^),  is  0.99  and  the  coefficient  of  variation  is  about  10%,  and 


/?„(520) 

a(490)=0.136(  " 


^„(560) 


) 


(8) 


B.  Curvature  vs.  Bottom  Reflectance 


Based  on  the  work  of  Gordon  et  al®,  Morel  &  Prieur^,  Carder  and 
Steward*  and  Lee  et  at,  the  remote-sensing  reflectance  for  optically 
shallow  water  is  expressed  as 

h  +h 

=0.  ^^[1  -e  +0. 173  pg  -2®“"  W 

Q*a 


where  and  are  back-scattering  coefficients  of  water  and  particles 
respectively,  a  is  the  total  absorption  coefficient,  Q  is  the  ratio  of  sub¬ 
surface  upwelling  irradiance  to  upwelling  radiance,  p  is  the  bottom 
albedo,  and  H  is  the  bottom  depth.  This  expression  excludes  the  terms 
for  CDOM  fluorescence  and  water  Raman  scattering,  since  they  play 
minor  roles  in  this  application.  Eq.9  is  then  used  to  define  two 
curvature  numbers: 


Cur.  1  = 


^,,(490)' 


RJ440)*RJ560) 


(10) 


rs' 


Cur. 2= 


RJ560')‘ 


R^^(490)*RJ620) 


(11) 


rs' 


For  SeaWiFS  usage  the  620nm  band  in  Cur,2  could  be  changed  to 


660nm,  without  altering  the  validity  or  usefulness  of  the  expression,  but 
we  suggest  the  addition  of  a  620nm  band  to  the  satellite  for  the 
following  two  reasons:  1)  Lv^(620)  has  a  higher  signal:noise  ratio  and 
less  influence  from  pigment  absorption  than  LJ660),  and  2)  with  this 
band,  the  remote-sensing  reflectance  curve  can  be  modeled  more 
closely  to  a  high  resolution  curve  (Figure  4). 

Using  Morel’s  modeF  to  provide  a  and  b^,,  a  Q  value^°  of  4.9,  and 
p=0.15*(X/400)^2  (a  general  assumption  based  on  the  measurement  of 
bottom  albedo  at  several  WFS  stations)  in  Eq.  9,  we  analyzed  the 
change  of  the  two  curvature  numbers  with  bottom  depth  and  pigment 
concentration.  Figures  5  and  6  are  the  results  of  this  analysis.  Figure  7 
shows  the  change  in  Cur.  1  for  optically  deep  water  versus  water  with  a 
20%  bottom  reflectance  contribution.  The  triangles  in  the  figure  are  the 
Cur.  1  values  of  the  optically  shallow  stations  in  the  WFS  data  set. 

The  20%  bottom  signal  is  defined  for  the  440,  490,  520,  560,  620 
and  660nm  wavelength  bands.  When  the  remote-sensing  reflectance  due 
to  the  bottom  is  20%  of  the  total  signal  for  any  of  these  bands,  then  we 
say  there  is  a  20%  bottom  signal.  Figure  8  illustrates,  at  low  and  high 
pigment  concentrations,  how  the  remote-sensing  reflectance  varies  when 
a  20%  bottom  signal  is  present. 

Figure  9  shows  how  the  derived  pigment  concentration  from  the 
"Morel  Case-l"  algorithm  varies  when  a  20%  bottom  reflectance 
contribution  is  included.  In  Figure  9,  percent  difference  is  calculated 


through 


where 


Percentdiff.  = 


C'-C 


xlOO 


/?(560) 


(12) 


(13) 


is  the  "Morel  Case-1"  algorithm^,  and  C  is  calculated  with  the  same 
equation,  but  with  a  reflectance  ratio  that  includes  a  20%  bottom 
reflectance  contribution. 


3.  Summary 


A.  Bottom  reflectance  in  optically  shallow  waters  will  cause 
significant  error  in  the  derived  pigment  concentration  through  its  effect 
on  the  radiance  or  reflectance  ratios. 

B.  For  shallow  coastal  waters  and  "Morel  Case-2”  waters,  an 
algorithm  that  uses  a  spectral  ratio  between  520  and  560nm  seems  a 
better  choice  to  remotely  derive  total  absorption  at  490nm. 

C.  This  remote-sensing  algorithm  (Equation  2)  performs  well  in 
both  clear-sky  and  overcast  conditions. 

D.  Curvature  analysis  of  the  remotely  sensed  reflectance  can 
provide  a  "flag"  for  bottom  reflectance  contributions  of  20%  or  more  of 
the  total  signal,  to  denote  pixels  likely  to  produce  significant  error  in 
pigment  concentration  calculations. 

E.  A  20%  bottom  signal  contributes  more  than  40%  error  for 
derived  pigment  concentrations  larger  than  2  mg/m^.  The  error  will  be 
even  greater  for  larger  bottom  signals. 
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Figure  1:  Reflectance  ratios  (from  data  uncorrected  for  bottom  reflectance  contribution) 
plotted  against  modeled  total  absorption  coefficient  at  490nm  (from  Lee  et  at 


Mod,  a(490) 


Figure  2.  ReflectHnce  ratios  from  Figure  1  plotted  against  the  sum  of  water  and  particulate 
absorption  at  490nm.  is  from  Smith  and  Baker'',  is  measured. 


("Y)"a/(\)“a 


+  a_  at  490nm 


Figure  4:  High  resolution  remote-sensing  reflectance  (solid  line),  and  SeaWiFs  spectral 
bands.  The  dotted  lines  are  water  column  and  bottom  renectance  components. 
Depth  is  7.5  m,  and  pigment  concentration  is  1.23  mg  m  ^ 
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Figure  6:  Surface  plot  of  Equation  1 1 .  Note  the  reversal  of  the  a(490)  axis  with  respect  to 
Figure  5. 


Figure  7:  Effect  of  bottom  reflectance  on  Equation  10,  as  a  function  of  total  absorption  at 
490nm. 


Figure  8:  Effect  of  bottom  reflectance  on  remote-sensing  reflectance.  Subscripts  rsw  and  rsb 
refer  to  the  remote-sensing  reflectance  signal  of  the  water  column  and  bottom, 
respectively.  [C]  is  the  pigment  concentration  (mg  m‘^). 
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ABSTRACT 

The  quantum  yield  <j>  of  photosynthesis  (mol  C  (mol  photons)  ')  at  6  depths  for  the  waters  of  the  Marine 
Light-Mixed  Layer  (ML-ML)  cruise  of  May  1991  were  calculated.  As  there  were  PAR  (photosynthetically  available 
radiation)  but  no  spectral  irradiance  measurements  for  the  primary  production  incubations,  three  ways  are  presented 
here  for  the  calculation  of  the  absorbed  photons  by  phytoplankton  (AF)  for  the  purpose  of  calculating  <t>.  The  first 
is  based  on  a  simple,  non-spectral  model;  the  second  is  based  on  a  non-linear  regression  using  measured  PAR  values 
with  depth;  and  the  third  is  derived  through  remote  sensing  measurements.  We  show  that  the  results  of  <j>  calculated 
using  the  non-linear  regression  method  and  those  using  remote  sensing  are  in  good  agreement  with  each  other,  and 
are  consistent  with  the  reported  values  of  other  studies.  In  deep  waters,  however,  the  simple  non-spectral  model  may 
cause  quantum  yield  values  much  higher  than  theoretically  possible. 
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INTRODUCTION 


Traditional  methods  for  estimating  the  light  absorbed  by  phytoplankton  at  depth  in  the  ocean  have  often 
used  nonspectral  PAR  (Photosynthetically  available  radiation)  measurements  or  calculations,  and  spectrally  averaged 
(without  light-field  weighting)  phytoplankton  specific-absorption  coefficients.  Two  alternative,  spectral  methods  are 
presented  here  for  use  with  remotely  sensed  data:  one  for  use  with  remote  floating  or  moored  arrays  equipped  with 
nonspectral  PAR  sensors,  and  one  for  use  with  satellite  or  aircraft-derived  remote-sensing  reflectance  data.  Both 
alternatives  utilize  spectrally  dependent  total  absorption  coefficients  to  model  the  depth  dependent,  spectral  light 
field.  These  light  fields  are  then  spectrally  integrated  at  depth  for  comparison  to  PAR  measurements. 

Primary  production  {PP)  in  the  ocean  is  dependent  on  the  absorption  of  incident  photons  by  phytopl.mkion 
and  on  the  efficiency  with  which  phytoplankton  use  this  absorbed  energy.  At  depth  z  for  a  time  period  of  T,  /’/’  is 

PP{z)  =  <i)(z)fJ^E^(zXt)[chl  a(z)]a*^(Zyk,t)ciXdt,  o 

(Symbols  are  defined  in  Table  1.  To  be  terse,  the  functionality  of  one  quantity  to  another  might  not  be  explicitly 
expressed  in  the  remaining  text.)  The  symbol  in  the  right  side  of  Eq.  1  describes  the  energy-conversion  efficicn..  y 
for  photosynthesis  and  is  defined  as  quantum  yield  of  photosynthesis  [Bannister,  1974;  Tyler,  1975;  Kiefer 
Mitchell,  1983;  Bannister  and  Weidemann,  1984;  Dubinsky  et  al.,  1984  and  Cleveland  et  al.,  1989],  The  revt  t 
the  right  side  of  Eq.l, 

AP(z)  =  fjE/z,XMchl  a(z)]a;,(zXt)dXdt,  >=• 

describes  the  absorbed  photons  by  phytoplankton  at  that  depth  for  the  time  period.  Eqs.  1  and  2  are  u>ujlly 
integrated  over  hours  or  a  day,  and  over  a  wavelength  range  from  400  to  700  nm.  So  quantum  yield  </>,  a  time 
average,  is 

In  order  to  estimate  primary  production  through  optical  properties,  using  either  moored  sensors  [Marm  i  t 
al.,  1992]  or  remote  measurements  [Cullen,- 1990,  Platt  et  al.,  1991;  Balch  et  al.,  1992],  the  magnitude  and 


variability  of  the  quantum  yield  of  photosynthesis  (^)  in  the  ocean  must  be  known.  Laboratory  studies  have  been 
carried  out  [Kiefer  and  Mitchell,  1983;  Langdon,  1988  and  references  cited],  but  there  is  a  paucity  of  in  situ 
estimates  of  <j>.  Ship-board  incubations  have  been  used  to  simulate  in  situ  studies,  but  most  have  used  white  light, 
which  can  cause  errors  as  large  as  a  factor  of  two  [Laws  et  al.,  1990].  Also,  m  many  cases,  non-spectral  models 
have  been  used  to  calculate  the  light  absorbed  by  phytoplankton  at  depth  for  ^  estimations  [Tyler,  1975;  Bannister, 
1979;  Marra  et  al. ,  1992].  However,  due  to  the  filtration  effect  of  the  overlaying  water  column,  non-spectral  models 
do  not  accurately  describe  the  light  absorbed  by  phytoplankton  at  depth  [Morel,  1978;  Kishino,  1986].  For  more 
accurate  in  situ  <f>  measurements,  the  spectral  interaction  between  the  light  at  depth  and  the  phytoplankton  absorj'Uon 
needs  to  be  considered. 

Bannister  and  Weidemaim  [1984]  and  Kishino  [1986]  presented  approaches  to  derive  in  situ  <i>.  Howcv  er. 
these  approaches  require  obtaining  knowledge  of  the  spectral  character  of  Efk,  z),  which  can  be  a  rather  expensive 
proposition:  either  i)  a  ship  must  remain  near  an  incubation  array  to  measure  the  spectral  light  field  at  ea..h 

incubation  depth,  or  ii)  a  spectral,  scalar-irradiance  meter  must  be  attached  to  the  floating  array  at  each  incubation 

\ 

depth. 

Neither  of  these  alternatives  were  foible  for  the  ML-ML  study,  although  spectrally  integrated  /’  \H 
measurements  (inexpensive)  were  made  at  four  depths  from  a  floating  array.  For  the  calculation  of  <^,  AP  at  i).  • 
incubation  depths  (different  from  the  depths  of  PAR  measurements)  were  calculated  by  three  simplified  methisJ^ 
The  first  is  a  typical  non-spectral  method  [e.g.  Marra  et  al.,  1992]  to  show  how  the  results  differ  from  spo-tul 
model  results.  The  second  is  through  use  of  a  non-linear  regression  method  for  interpolating  PAR  values  between 
and  below  PAR  measurements.  The  third,  based  upon  surface  PAR  and  absorption  coefficients  derived  from  surtax e 
remote-sensing  measurements,  is  applicable  to  ocean  color  data  collected  from  spacecraft  or  aircraft.  These  method' 
were  tested  at  four  Marine  Light-Mixed  Layer  stations,  and  the  differences  among  the  results  are  discussed  and 


evaluated. 
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MEASUREMENTS 


Data  of  PAR{z),  PP{z),  a;*(X)  (absorption  coefficient  of  phytoplankton)  and  /?„(X)  (remote-sensing 
reflectance)  were  collected  on  an  ML-ML  cruise  from  May  17  to  May  24,  1991  in  the  waters  south  of  Iceland.  The 
water  depths  where  PAR  and  PP  were  measured  were  not  the  same  (see  subsections  a  and  b  below),  although  they 
were  carried  out  on  the  same  floating  array  for  four  days:  May  17,  20,  22  and  24.  The  sampling  site,  conditions, 
wind,  mixing,  and  nutrients  are  presented  in  Marra  et  al.,  Jones  et  al.  and  Plueddemann  et  al.  [all  in  this  issue]. 

a.  Photosynthetically  available  radiation,  PAR{z) 

PAR(z)  at  four  depths  (0,  2,  12.5,  and  25  m)  were  monitored  and  averaged  for  every  10  minute  inierv  al 
throughout  the  day  with  a  Biospherical  PAR  sensor  attached  to  each  incubation  array  [see  Marra  et  al. ,  this 


b.  Primary  production,  PPiz) 

Dawn-to-dusk  (17  hours)  incubations  with  four  replicates  were  carried  out  in  situ  at  each  of  six  depths  (5, 
10,  15,  20,  30  and  40  m)  chosen  to  span  the  euphotic  zone.  Primary  production  measurements  were  made  using 
the  '^C  technique  [see  Marra  et  al.,  this  issue]. 

c.  Absorption  coefficient  of  phytoplankton  pigments,  a^(X) 

Measurements  of  a,*(X)  for  replicates  of  the  incubation  samples  were  obtained  in  the  following  manner. 
The  method  described  by  Mitchell  &  Kiefer  [/95S]  was  used  to  measure  the  particulate  absorption  coefficients,  with 
a  Lambertian  diffuser  added  between  the  light  source  and  the  filter  pad  to  minimize  any  dependence  of  the  diffuse 
transmittance  on  the  optical  geometry  of  the  light  which  is  incident  on  the  pad.  This  arrangement  is  very  similar 
to  that  of  Bricaud  &  Stramski  [1990],  who  illuminated  a  pad  in  front  of  a  diffusing  window  which  was  adjacent  to 
an  end-on  photomultiplier  tube.  Our  collector  was  a  Spectron  Engineering  Spectroradiometer  with  10°  acceptance 
angle.  We  filtered  about  300  ml  water  samples  onto  Whatman  GF/F  pads,  the  transmission  spectra  (380-780nm) 
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of  these  pads  were  ratioed  against  that  of  wetted  blank  pads  to  obtain  the  optical  density  of  the  particles,  and  this 
optical  density  was  corrected  for  large  particle  scattering.  The  particulate  absorption  coefficient  was  calculated,  using 
Eq.  2  from  Bricaud  &  Stramski  [1990]  to  correct  for  optical  pathlength  elongation  due  to  scattering  in  the  pad  (”/J 
factor").  After  this  measurement  each  sample  was  soaked  in  hot  methanol  [Kishino  et  al.,  1985;  Roesler  et  al., 
1989]  for  about  IS  minutes  to  remove  pigments,  and  its  transmission  spectra  was  again  measured  and  ratioed  against 
a  methanol  wetted  blank  to  obtain  the  detrital  absorption  coefficient  (a,XX)).  The  difference  between  the  particle  and 
detrital  absorption  coefficients  provided  the  absorption  coefficient  for  phytoplankton,  a^(X). 

Gelbstoff  absorption  coefficients  for  surface  waters  were  reported  in 
Walsh  et  al.  [1992],  who  filtered  seawater  through  0.2  /*m  Gelman  Supor-200  filters  and  measured  the  absorption 
on  a  Cary  model-2200  spectrophotometer  with  a  10  cm  cell.  The  approach  followed  that  of  Bricaud  et  al.  1/95/] 
where  absorption  at  370  nm  was  extrapolated  to  440  nm  as  a, (440)  =  a,(370)e‘°®'^^^™’  with  an  average  semi-log, 
spectral  slope  of  -  0.014  nm'. 

d.  Remote  sensing  reflectance,  /?„(X) 

Hyperspectral  remote-sensing  reflectance  /?„(X)  (the  ratio  of  water  leaving  radiance  to  the  downwelling 
irradiance  above  the  sea  surface)  was  measured  at  noon  for  each  station  by  the  method  developed  by  Carder  and 
Steward  [1985],  using  the  Spectron  Engineering  Spectroradiometer.  It  has  253  chaimels  and  covers  the  wavelength 
range  from  370  -  1100  nm.  With  this  instrument  the  water-leaving  radiance  and  downwelling  sky  radiance  v^e^e 
directly  measured,  and  downwelling  irradiance  was  measured  by  viewing  a  standard  diffuse  reflector  (Spectralon. 
~8.5%).  Sky  radiance  reflected  from  the  sea  surface  was  removed  by  the  method  of  Carder  and  Steward  1/955| 
when  calculating  R„{\). 

AP  and  <j>  CALCULATION 

a.  Simple,  non-spectral  model. 

Due  to  the  expense  of  obtaining  the  spectral  light  field  Eji\)  with  depth,  a  simplified  model  for  PP  is 
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usually  used  [Tyler,  1975;  Bannister,  1979;  Morel,  1991;  Marra  et  al,  1992], 

PP(z)  =  (|)(z)[cW  a(z)]a-^(z)PAR(z),  (4) 

where  a^*  (m^  (mg  chi  a)  ')  is  the  spectrally  averaged  chlorophyll-specific  absorption  coefficient: 


~  [cU  fl](700-400)  ’ 


PARiz)  =  f  r^E(k,t,z)dXdt.  (6) 

J  rJ400  * 

For  the  calculation  of  PAR  at  incubation  depths,  was  calculated  first  according  to 


'■fAJt 


MPARCZj)]  -  ln[PX/?(z^)] 


(7) 


and  then  the  value  between  2  and  12.5  m  is  used  to  calculate  PAR(5)  and  P/4f?(10),  and  the  value  between  12.5 
and  25  m  is  used  to  calculate  FAR(15),  PAR(20),  PAR(30)  and  PAR(40). 

Thus,  from  Eqs.  1  and  2,  we  have 


♦(z)  = 


PP(Z) 

[chi  a(z)]ayARCz)' 


(8) 


Since  PP,  [chi  a]  and  were  measured,  it  is  easy  to  calculate  this  way.  However,  as  has  been  pointed 
out  by  Morel  [1978]  and  Kishino  [1986],  a  simple  average  for  (Eq.  5)  does  not  accurately  describe  the 
phytoplankton  absorption  at  depth,  and  a  light-field-weighted  average  needs  to  be  applied  to  Thus,  a  better 
approximation  that  includes  the  character  of  light  at  depth  is  necessary.  The  following  two  sections  describe  two 
ways  to  model  the  spectral  character  of  the  light  field  at  depth  when  only  PAR(z),  or  i^(X)  and  z)  are 
available. 


b.  Non-linear  regression. 
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At  depth  z,  daily  PAR  can  be  expressed  as 

PAR{z)  -  I ^E^(Q,Xj)e-'^^'^^d\eb  ,  (9) 

where  £,(0,X,t)  is  the  sub-surface  scalar  irradiance  at  time  t  in  units  mol  quanta  m'^s  'nm  ',  {XT}  is  the  vertically- 
averaged  diffuse  attenuation  coefficient,  and  {^l}  is  approximated  by  {K)  =  1.08(a  +  [Kirk,  1984;  Gordon, 
1989].  In  Eq.  9,  since  most  of  the  stations  were  cloudy,  from  the  Sun  and  from  the  sky  are  combined.  For  clear 
sky  situations,  separation  of  these  two  may  be  necessary  since  the  average  cosine  values  for  sun  light  and  sky  light 
are  generally  different  [Platt  et  al. ,  1991].  Since  there  were  no  time  series  of  £^(0,X)  and  n^(X),  some  assumptions 
have  to  be  made  in  order  to  calculate  the  spectral  light  field  at  depth: 

i) .  Gelbstoff  absorption  was  constant  for  all  depths; 

ii) .  Average  of  a^,  for  dawn  and  noon  measurements  was  considered  to  represent  the  pariule 
absorption  coefficient  for  the  whole  day; 

iii) .  Due  to  the  presence  of  clouds,  it  is  difficult  to  estimate  fig  for  any  given  time.  So  an  average 
sub-surface  average  cosine  of  ^g  =  0.77  was  used  for  these  high-latitude  sites  for  the  whole  day,  including  sunny, 
cloudy  and  overcast  conditions; 

iv) .  bg  <<  a  [Morel  and  Prieur,  1977]. 

Then  we  have 


PARiz)  -  ,  <l'»i 

where  1.4  in  the  exponent  is  »  1.08/0.77,  and  integration  over  time  is  implicitly  combined  into  £,(0),  and  TA.- 
z,  and  -  z,). 

In  Eq.  10,  vertical  stratification  of  a^  is  considered  although  most  of  the  stations  were  quite  well  mixed  [  w  c 
Marra  et  al.,  this  feue].  Also,  for  Eq.  10,  PAR{z)  and  aj^z)  were  measured,  is  known  [Smith  and  Baker,  l'^si\. 
and  Eg{0)  and  a^  must  be  derived. 

aj(X)  can  be  expressed  as  fl,(X)  ==  aj(440)e'®®'^“**”  [Bricaud  et  al.  1981].  Since  this  variable  is  \cry 
difficult  to  measure  at  visible  wavelengths  for  open  ocean  waters  [Bricaud  et  al.,  1981],  Walsh  et  al. 
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measured  a, (370)  and  used  the  above  expression  to  extrapolate  to  a, (440)  for  waters  at  the  ML-ML  site.  The 
spectral  slope  0.014  nm  '  for  is  assumed  and  kept  as  constant  for  all  the  stations. 

£,(0)  is  expressed  as  £^(0)  =  PA/?(0)j„  *  £,*(0),  where  /Vl/?(0)j„  is  the  derived  surface  PAR  in  order  to 
correct  possible  errors  in  measured  PAR(z)  due  to  sea-surface  roughness,  and  £,*(0)  =  E^(0)/PAR(0).  £,(0)  and 
PAR(0)  are  computed  using  the  model  of  Gregg  and  Carder  {1990\.  There  were  typically  4  (3  for  May  22)  PAR 
measurements  for  each  day,  so  PAR(0)^,  and  a, (440)  for  each  day  can  be  derived  by  non-linear  regression:  i.e. 
minimize  the  variance  between  the  modeled  and  measured  PAR.  The  derived  0^(440)  values  are  shown  in  Table  3. 
After  the  derivation  of  PAR(0)^,  and  <1,(440),  quantum  yield  can  be  calculated  for  each  depth  through 


<l>(z)  = 


_ _ 

PAR(0) ' 


(Ml 


c.  Remote  sensing  reflectance 

With  recent  measurement  techniques  [e.g.  Carder  and  Steward,  1985]  and  models  for  interpreting 
hyperspectral  remote-sensing  reflectance  (£„)  data  [Carder  and  Steward,  1985;  Peacock  et  al.,  1990;  Lee  et  nl. . 
1992;  Lee  et  al.,  1994(a),  1994(b)],  it  has  been  shown  that  by  inversion  of  the  measured  R„,  the  total  absorption 
coefficient  for  the  surface  layer  can  be  derived  [Lee  et  al.,  1994(a),  1994(b)].  This  measurement  approach  a\oiJs 
the  inherent  problems  of  in-water  reflectance  measurements  such  as  instrument-self  shading  [Gordon  ami  Dm.:. 
1992]  and  ship  shadows  [Gordon,  1985]. 

To  use  this  approach,  i.e.  deriving  the  total  absoiption  coefficient  and  then  calculating  the  photons  ab.sorK  J 
by  phytoplankton  at  depth  for  a  day,  we  still  need  most  of  those  assumptions  made  for  the  non-linear  regression 
method.  However,  PAR(0)^,  is  set  to  equal  the  measured  surface  PAR,  and  the  total  absorption  a(X)  is  derived  by 
inverting  R„(\)  data  [see  Lee  et  al.,  1994(a)].  One  approach  for  deriving  a(\)  from  R„  is  suggested  by  Gordon  et 
al.  [79SS],  Morel  [1988],  and  Sathyendranath  and  Platt  [1988].  They  empirically  derive  the  chlorophvll 
concentration  first  and  use  an  empirical  specific  absorption  coefficient  to  calculate  the  particle  absorption  coeffit  u  ni 
The  gelbstoff  absorption  coefficient  is  assumed  to  covary  with  [chi  a].  This  method  only  works  well  for  "case  1 " 
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waters  [Morel  and  Prieur,  1977\. 

Recent  studies  [Gordon  et  al. ,  1988',  Morel  and  Gentili,  1993\  Lee  et  al. ,  1992,  1994(a)]  show  that  R„  can 
be  expressed  as 


G  f,  \  ft 

a+a*a^  k 


(12) 


r  t 


where  G  »  0.093(T/n)’,  the  air-sea  interface  transmittance  t  *  0.97  for  most  of  the  ocean  [Austin,  1974],  and  n 
»  1.341  is  the  refractive  index  of  sea  water.  The  two  terms  in  brackets  consist  of  scattering  contributions  due  to 
water  molecules  and  particles,  respectively.  The  back-scattering  coefficient  of  pure  water,  b^,  is  given  by  Smith 
and  Baker  [1981],  and  X  and  Y are  site  specific  and  can  be  derived  from  R„  [Lee  et  al.,  1994(a),  1994(b)]. 

By  adjusting  parameters  affecting  the  modeled  R„  curve  to  achieve  a  close  fit  with  the  measured  R„  curve, 
aj(X)  of  the  surface  layer  was  derived.  The  derived  a^(440)  values  are  shown  in  Table  3.  In  the  derivation  ot  a^ 
using  R„,  a  multiplication  factor /was  applied  to  the  surface  a,  as  in  Lee  et  al.  [1994(a)].  For  the  ML-ML  waters 
/varied  from  1.0  -  1.2.  This  factor  accounts  for  possible  errors  in  particle  absorption  due  to  plankton  patchines.s, 
vertical  stratification,  and  the  factor". 

After  the  inversion  of  R„,  the  sum  a„  +  f*a^  +  o,  was  used  as  the  average  total  absorption  coefficient 
for  the  water  column  to  calculate  PAR  and  AP  at  depth  using  measured  profiles  of  U;*(z).  Quantum  yield  values  were 
then  calculated  based  upon  measured  values  of  primary  production. 


RESULTS  and  DISCUSSION 

Figure  1  shows  the  pre-dawn  (open  symbols)  and  noon  (filled  symbols)  noeasured  chlorophyll  values  lor 
the  cruise.  Notice  that  there  was  no  strong  vertical  structure  in  the  water  column,  if  the  average  of  the  pre-daw  n 
and  noon  profiles  are  considered.  The  differences  between  pre-dawn  and  noon  chlorophyll  values  were  larger  lor 
May  17  and  20,  but  much  smaller  for  May  22  and  24,  however. 

Figure  2  shows  the  measured  primary  production  for  each  of  the  four  days.  The  surface  production 
decreased  from  Day  1  to  Day  3.  This  might  be  due  in  part  to  the  progressive  reduction  in  chlorophyll  with  time. 
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but  the  variation  in  PAR{0)  values  in  Table  3  suggests  a  more  complicated  explanation. 

Figures  3(a)  and  3(b)  show  the  measured  spectra  of  remote-sensing  reflectance  and  the  pre-dawn  particle 
absorption  coefficients  for  the  upper  water  column,  respectively.  May  17  had  the  lowest  R„  values  around  440  nm, 
consistent  with  the  highest  noon  chlorophyll  concentration  among  the  4  days  (Figure  1).  Also  notice  the  high 
/?„(685)  values  for  May  17,  likely  resulting  from  the  contribution  of  chlorophyll  a  fluorescence.  For  the  following 
3  stations,  the  R„  values  were  similar  as  their  surface  noon  chlorophyll  concentrations  did  not  vary  much  (Figure 
1). 

The  detritus-to-pigment  ratio  for  absorption  at  440  nm  was  in  the  range  of  15  to  20  %  with  an  average  of 
about  18  %.  The  spectral  curves  were  simulated  using  exponential  functions  with  an  average  semi-log  spectral  slope 
for  of  about  -  0.010  nm  ‘,  similar  to  that  found  by  Roesler  et  al.  [1989\. 

Figures  4(a)  to  4(d)  show  the  measured  and  calculated  PAR  for  the  four  days.  In  general,  PAR  values  at 
greater  depths  calculated  by  the  simple,  non-spectral  model  were  lower  than  those  determined  by  non-linear 
regression  and  by  R„  methods.  One  reason  why  PAR  estimates  at  greater  depths  were  lower  using  the  non-spectral 
model  is  that  Beer’s  Law  does  not  hold  for  Kp^  [Gordon,  1989].  The  PAR  values  determined  by  non-linear 
regression  and  R„  methods  were  quite  close  (average  difference  of  9%).  Compared  to  the  PAR  values  at  2  m,  12.5 
m  and  25  m,  the  average  difference  between  the  measured  and  /?„-method  calculations  was  19%.  Most  of  the 
difference  was  found  at  2  m  (see  figure  4),  where  the  measured  values  appeared  to  be  low,  especially  on  May  24 
However,  it  is  not  clear  what  caused  this. 

Figures  5(a)  to  5(d)  show  the  calculated  quantum  yield  <j>  values  for  the  four  days.  The  simple,  non-spectral 
model  generates  the  highest  <i>  values  at  30  and  40  m,  some  of  which  are  much  higher  than  the  0.12  mol  C  (mol 
photons)  '  theoretical  maximum  [Kok,  1960].  Possible  reasons  for  these  over-estimates  include:  1)  non-consideration 
of  the  spectral  character  of  the  light  field  [Morel,  1978;  Kishino,  1986]\  and  2)  F/1R(30)  and  P/1R(40)  were  deri\  ed 
by  extrapolation  using  Kp^,  which  does  not  follow  Beer’s  Law. 

The  differences  in  derived  <f>  values  for  the  non-linear  regression  and  R„  methods  are  2  to  28%  (average 
11%),  with  most  of  the  differences  occurring  at  40  m  near  the  bottom  of  the  euphotic  zone.  This  difference  is 
largely  due  to  the  different  derived  a,  values.  As  there  were  no  PAR  measurements  at  30  and  40  m,  it  is  difficult 
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to  judge  which  method  is  more  accurate.  The  values  determined  from  the  R„  model  inversion,  however,  are 
closer  to  the  measured  values  reported  by  Walsh  et  al.  [1992]  with  an  average  difference  of  only  about  18%.  The 
calculated  if>  values  by  both  methods  are  consistent  with  values  reported  elsewhere  [T)>ler,  1975;  Dubinsky  et  al., 
1984;  Kishino,  1986],  and  are  less  than  the  0.12  mol  C  (mol  photons)  '  theoretical  maximum. 

Part  of  the  error  in  4>  calculation  comes  from  estimates  of  PAR{z)  and  measurements  of  a^(X).  As  the  error 
of  PAR  was  19%  using  the  R„  method,  if  the  error  in  measuring  was  10%,  then  the  error  in  <f>  would  be  about 
22%. 

The  error  in  <l>  caused  by  errors  in  a,  varies  with  the  relative  contribution  of  to  the  total  absorption,  and 
also  varies  with  depth.  For  the  waters  studied,  a  100%  error  in  a, (440)  would  cause  a  factor  of  2  error  in  4>  at  the 
bottom  of  the  euphotic  zone,  whereas  it  would  only  contribute  an  error  of  about  10%  in  for  the  surface  layer. 
However,  a  recent  study  [Lee,  1994]  shows  that  the  average  error  of  /?„-derived  total  absorption  is  only  about  30%. 
Therefore  the  error  in  is  expected  to  be  much  less  than  a  factor  of  2,  especially  when  a^(X)  is  measured. 

Figure  6  shows  calculated  versus  PAR  values;  (a)  results  using  the  non-linear  regression  method;  and 
(b)  results  by  the  R„  method.  These  figures  are  similar  in  shape  to  those  reported  by  Tyler  [1975]  and  Kishino 
[1986],  If  we  ignore  the  <5  values  at  40  m  as  being  less  reliable,  <j>  values  for  both  methods  can  be  expressed  as 
suggested  by  Steele  [1962]-. 

(t)  =  ,  (13) 

where  <i>„  and  a,  derived  from  the  results  of  the  two  methods,  are  shown  in  Table  2  as  a  function  of  time.  All 
the  values  are  consistent  with  the  maximum  quantum  yield  values  reported  by  Welschmeyer  and  Lorenzen  [1981 ) 
(0.040  -  0.074  mol  C  (mol  photons)  '),  Kiefer  and  Mitchell  [1983]  (0.05  mol  C  (mol  photons)  '),  Langdon  [79551 
(0.034  -  0.10  mol  C  (mol  photons)  ')  and  Cleveland  et  al.  [7959]  (0.033  -  0.102  mol  C  (mol  photons)  ').  It  is  clear 
that  the  differences  in  and  a  as  calculated  by  the  two  methods  are  smaller  than  the  differences  due  to  changes 
in  the  water  column  over  the  four  days.  The  values  for  the  last  three  sampling  days  were  fairly  consistent  and 
about  40%  larger  than  on  the  first  day,  but  the  slope  value  a  increased  more  than  two-fold  between  the  first  two 
days  before  dropping  to  an  intermediate  level.  The  daily  change  of  <l>„  and  a  can  be  caused  by  changes  in 
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physiological  parameters,  such  as  nutrients,  adaptation,  species,  and  light  level.  For  this  cruise,  however,  nutrients 
were  always  abundant  [Marra  et  al.,  this  issue]. 

Table  3  summarizes  the  above  results  along  with  the  column-integrated  (within  the  euphotic  zone) 
chlorophyll  a  and  primary  production.  The  column-integrated  chlorophyll  a  decreased  from  the  beginning  to  the  end 
of  the  cruise.  Also,  for  the  four  days,  /VU?(0)  and  chlorophyll  a  changed  markedly  (*  4  X  and  3  X,  respectively), 
but  the  integrated  primary  production  did  not  vary  as  much  ( <  2  X).  This  was  especially  interesting  for  the  last 
three  days,  which  included  a  sutmy  day  (May  22)  when  the  productivity  near  the  surface  on  that  day  appeared  to 
have  been  significantly  reduced  by  photoinhibition. 

May  22  was  the  only  sunny  day  for  the  entire  cruise,  and  the  populations  present  would  have  been  l()^v  licht 
adapted.  The  three  to  four-fold  reduction  in  surface  <f>  values  for  the  rare  sunny  day  compared  to  values  tur  ihe 
previous  two  stations  suggests  that  great  care  must  be  taken  in  interpreting  satellite  estimates  of  primary  pr()du>.iii>n 
in  light-changing  environments.  The  occasional  view  of  the  ocean  a  satellite  observes  on  a  rare  sunny  day  may  be 
accompanied  by  perturbations  in  the  photosynthetic  performance  of  the  phytoplankton  because  of  a  light-hisiury 
factor  that  is  not  well  understood  at  this  time. 

SUMMARY 


1.  Three  methods  of  estimating  the  absorbed  photons  at  depth  were  evaluated,  one  of  which  was  n.  n 

spectral. 

2.  The  spectral  methods  (the  non-linear  regression  and  remote  sensing  reflectance  approaches)  both  pri)v  i  Jc 
good  alternatives  to  model  the  spectral  light  field  at  depth  for  reasonably  well  mixed  water  columns  such  as  the  M I.- 
ML  site.  The  PAR  values  for  the  non-spectral  method  were  too  low  and  resulted  in  theoretically  unrealistic  quantum 
yield  values. 

3.  Calculated  quantum  yield  values  by  the  two  spectral  methods  are  consistent  with  each  other  and  uuh 
values  reported  by  other  studies.  The  equation, 

provides  an  estimate  of  the  average  quantum  yield  with  PAR  for  the  ML-ML  site.  Combining  surface  PAR  valucN. 
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4i(z)  »  ,  (14) 

either  measured  or  modeled  [e.g.  Gregg  and  Carder,  1990;  Bishop  and Rossow,  1991],  with  Eq.  9  or  Eq.  10  allows 
PAR(z)  calculations  to  be  made  on  the  basis  of  remotely  sensed  data. 

4.  The  daily  PAR{0)  varied  four-fold,  and  the  column-integrated  chlorophyll  a  varied  three-fold  during  the 
cruise.  The  column-integrated  primary  production  varied  only  two-fold,  however.  This  argues  for  the  need  to  better 
understand  factors  affecting  the  quantum  yield  and  the  light  field  with  depth  in  order  to  improve  model  simulations 
of  primary  production. 

5.  The  methods  suggested  here  for  estimating  of  ij>  provide  alternatives  for  spectral  measurements  uhen 
only  data  of  PAR(z)  or  R„  are  available.  If  feasible,  time-series  measurements  of  the  spectral  light  field  ai  ihe 
incubation  depths  would  be  the  best  choice.  Also,  time  series  of/?„(X),  £,(0,X),  and  a^(z,\)  would  help  to  impo've 
the  suggested  methods  and  to  understand  in  situ  <f>. 

6.  The  reduction  exhibited  in  surface  0  values  for  a  rare  sunny  day  suggests  that  primary  produetmn 
estimates  based  upon  satellite  measurements  of  ocean  color  for  regions  with  a  history  of  cloud-cover  may  he  m 
serious  error  if  light  history  of  the  region  is  not  considered  in  assessing  effects  of  photoinhibition. 
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FIGURE  LEGENDS 


1.  Depth  profiles  of  chlorophyll  values. 

Open  symbols:  pre-dawn  measured,  filled  symbols:  noon  measured. 

2.  Depth  profiles  of  measured  production  values. 

3.  Measured  spectra  of  water  properties. 

3(a):  remote-sensing  reflectance 
3(b):  particle  absorption  coefficients 

4.  Depth  profiles  of  measured  PAR  vs.  calculated  PAR. 

4(a):  May  17,  4(b):  May  20,  4(c):  May  22,  4(d):  May  24. 

5.  Depth  profiles  of  quantum  yield. 

5(a):  May  17,  5(b)  May  20,  5(c):  May  22,  5(d):  May  24. 

6.  Quantum  yield  vs.  PAR. 

6(a):  calculated  by  non-linear  regression  method, 

6(b):  calculated  by  R„  method. 
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Table  1:  Symbols  and  units 


fl(X) 

a/X) 

a„(X) 

AP 

VX) 

£;(X) 

£/x) 

^PAK 

{m} 

n 

PAR 

PP 

PARiO)^, 

RrsfX) 

T 

a 

<t> 

X 

T 

1^0 


Total  absorption  coefficient,  m  ';  a(X)  =  aJX)  +  fl,(X)  +  a^fX) 

Absorption  coefficient  of  detritus,  m  ' 

Absorption  coefficient  of  gelbstoff,  m  ' 

Absorption  coefficient  of  phytoplankton  pigments,  m  ' 

Absorption  coefficient  of  pure  water,  m  ' 

Absorbed  photons  by  phytoplankton,  mol  photons/m?/day 
Backscattering  coefficient,  m  ' 

Specific  quantum  irradiance,  £^(X)//Vl/?(0),  nm  ' 

Quantum  irradiance,  /imol  photons/mVnm/s 
Attenuation  coefficient  for  PARy  m  ' 

Vertically-averaged  attenuation  coefficient  for  downwelling  irradiance,  m  ‘ 

Refractive  index  of  sea  water 

Photosynthetically  available  radiation  (integrated  from  400  -  700  nm),  mol  photons/m^/day 
Primary  production,  mol  C/mVday 
Derived  surface  PAR,  mol  photons/m^/day 

Remote  sensing  reflectance,  ratio  of  upwelling  radiance  to  downwelling  irradiance  above  the  sea 
surface,  sr'‘ 

Integration  or  incubation  period  of  time,  hours  or  day 
Slope  to  describe  ^  changes  with  PAR,  (mol  photons/mP/day)  ' 

Quantum  yield,  mol  C  (mol  photons)  ' 

Maximum  quantum  yield,  mol  C  (mol  photons)  ' 

Wavelength,  nm 

Air-sea  surface  transmittance 

Sub-surface  average  cosine  for  downwelling  irradiance 


Table  2.  Calculated  maximum  quantum  yield 


Table  3:  Summary  of  Results 
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Abstract 

A  mathematical  simulation  of  the  absorption  coefficient  of  phytoplankton  pigments  is  used 
to  estimate  the  primary  production  in  the  euphotic  zone  based  on  the  measured  remote-sensing 
reflectance  and  surface  solar  irradiance.  In  the  suggested  approach  here,  it  is  not  necessary  to 
know  the  pigment  concentration  and  the  chlorophyll-specific  absorption  coefficient,  which  are 
essential  in  traditional  methods.  This  approach  is  tested  against  measurements  made  in  the  high- 
latitude  North  Atlantic  in  May  1991.  A  comparison  of  the  calculated  and  measured  daily  primary 
production  in  the  euphotic  zone  had  an  r^  =  0.95  (n=24),  with  a  slope  of  1.26  and  -t-  32% 
difference.  Using  the  traditional  method  yielded  an  =  0.85  (n=24)  with  a  slope  of  0.34  and  - 
78%  difference.  These  results  mean  that  there  is  a  factor  of  3  improvement  in  estimation 
accuracy  for  the  phytoplankton  photosynthesis  of  waters  studied;  it  also  suggests  that  it  is 
preferable  to  shift  the  primary  production  model  from  traditional  pigment-concentration  base  to 
pigment-absorption  base. 
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Introduction 


In  the  past  decades,  primary  production  models  are  based  on  pigment  concentration, 
either  using  light  available  expression  (e.g.  P-I  relation  in  Platt  et  al.  1991)  or  light  absorbed 
expression  (e.g.  P-AQ  relation  in  Bidigare  et  al.  1992),  or  by  non-spectral  or  spectral  model  (see 
Bidigare  et  al.  1992).  Recent  popular  models  use  the  light  absorbed  approach,  i.e.,  chlorophyll 
a  concentration  multiplied  by  the  irradiance  and  two  factors:  averaged  chlorophyll-specific 
absorption  coefficient,  and  the  quantum  yield  for  carbon  fixation.  Thus  the  photosynthesis  rate 
at  depth  z  is  (Kishino  et  al.  1986;  Smith  et  al.  1989;  Cullen  1990;  Morel  1991;  Marra  et  al. 
1992;  Zaneveld  et  al.  1993): 


PPiz)  -  michl  a]fl;?A/?(z)  ,  (I) 

where  4>(z)  is  the  quantum  yield  for  carbon  fixation,  and  [chi  a]  is  the  chlorophyll  a 
concentration  in  mg/m^  (important  symbols  used  in  this  text  are  summarized  in  Table  1). 

~a^  is  the  spectrally-averaged  chlorophyll-specific  absorption  coefficient  over  400  -  700 
nm  (m^(mg  chi  fl)  ‘): 


P\;(X)E^{z,X)d\ 

— .  J400  ’ 

^  —  - 


(2) 


and  fl^*(X)  =  fl^(X)/[chl  a]  is  the  chlorophyll-specific  absorption  coefficient. 
PAR{i)  is  the  photosynthetically  available  irradiance,  computed  as 


PAm  -  f^E,(z,).)dX.  , 
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(3) 


where  E„(z.\)  in  equations  2  and  3  is  the  spectral  scalar  irradiance  in  quanta/mVnm/s.  From 
Sathyendranath  and  Platt  (1988), 

E,(z,X.)  -  .  (^) 

Recently,  Bannister  (1992)  pointed  out  that  Equation  4  is  a  rough  approximation  for  high- 
scattering  media.  However,  the  waters  studied  here  show  a  low  scattering-to-absorption  ratio, 
and  the  sub-surface  solar  zenith  angle  was  large.  That  means  there  was  not  much  vanation 
between  the  sub-surface  average  cosine  and  the  average  cosine  at  depth  (see  Fig.  3  of  Bannister 
1992).  So,  equation  4  can  still  be  used  for  the  situation  described,  especially  when  an  vertically- 
averaged  Kj  is  used. 

Quantum  yield  (Langdon  1988;  Cleveland  et  al.  1989;  Smith  et  al.  1989),  0,  varies  with 
light  intensity  and  physiological  status  of  the  phytoplankton  population.  We  chose  the  empirical 
formula  suggested  by  Kiefer  and  Mitchell  (1983)  to  express  how  <t>  changes  with  PAR, 


<t>(z)  =  4>„ - - - , 

”'PAR^^PAR(z) 

where  <{>„  is  the  maximum  quantum  yield,  PAR^  is  the  PAR  value  at  which  <f>  =  <I>J2. 

Therefore,  combined  with  the  photoinhibition  expression  suggested  by  Platt  et  al.  (1980), 
primary  production  at  depth  z  can  be  expressed  as 


PP(z)  =  APP(z) 


PAR^^PARiz) 


-VPAIK.Z) 

f 


(6) 


with 


APP{z)  -  p[q^;A4i?(z) .  (7) 

where  j/  is  a  parameter  to  describe  photoinhibition,  the  influence  of  which  will  be  discussed  in 
the  section  on  results  and  discussion. 

Eq.  7  accounts  for  the  absorbed  photons  by  the  phytoplankton  pigments,  in  which  p 
determines  the  fraction  of  chlorophyll  a  relative  to  chlorophyll  a  +  phaeophytin  a  concentration 

[q. 

Often,  estimation  of  primary  production  based  upon  remote  measurements  are 
accomplished  using  Eqs.  6  and  7.  Traditionally,  [q  is  estimated  first  (one  of  the  important 
reasons  for  the  development  of  pigment  algorithm  from  remote  sensing)  from  remotely  measured 
signals  (Platt  et  al.  1991;  Balch  et  al.  1992)  using  for  example  the  CZCS  (Coastal  Zone  Color 
Scanner)  algorithm: 


^  LJSSO)  '  RJ550) 


where  A  =  1.13  mg/m^  and  5  =  -1.71  (Gordon  et  al.  1983).  The  0.95  comes  from 
£’d(443)/£j(550)  »  0.95.  Using  [q,  the  diffuse  attenuation  coefficient  of  irradiance  can  be 
approximated  by  the  empirical  relationship  suggested  by  Morel  (1988),  which  is: 


K/X)  -  K^(Xhd(X)lC]^^^  ,  (9) 

in  which  values  for  K„,  d  and  e  are  found  in  tables  calculated  by  Morel  (1988). 

With  Kj  and  measured  or  calculated  £'„(0)  (e.g.  converted  from  Gregg  and  Carder  (1990) 
or  Bishop  and  Rossow  (1991)),  PAR  at  any  depth  can  then  be  estimated.  However  p,  p, 
PAR^  and  a^*  must  be  estimated  for  the  calculation  of  PP,  p  varies  from  0.3  to  0.9,  with  an 
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average  value  of  0.75  (Morel  and  Berton  1989;  Balch  et  al.  1992),  the  value  used  for  our 
calculation.  >,  and  PAR^  vary  with  phytoplankton  physiological  status,  which  at  this  point 
cannot  be  derived  based  upon  remotely  sensed  measurements,  but  must  be  estimated  from  other 
information.  <f>„  has  been  reported  to  range  from  0.03  mol  C  (Ein  absorbed)  ’  to  0. 1  mol  C 
(Ein  absorbed)  ’  (Bannister  and  Weidemann  1984;  Smith  et  al.  1989;  Morel  1991).  As  a  kind 
of  average  for  productive  waters,  is  assumed  to  equal  0.074  mol  C  (Ein  absorbed)  ’,  a  value 
suggested  by  Cullen  (1990).  In  Marra  et  al.  (1993),  however,  a  4>„  value  of  0.06  mol  C  (Ein 
absorbed)*’  was  used.  There  is  little  literature  information,  however,  about  the  values  of  v  and 
PAR^.  Platt  et  al.  (1980)  found  p  varied  from  0  to  0.0028  (W  m*^)*’.  If  we  choose  0.0028 
(W  m*^)*’  for  p,  its  equivalent  value  is  »>  *  0.01  (Ein/mVday)*’.  Kiefer  and  Mitchell  (1983) 
found  that  PAR^  equals  about  10  Ein/m^/day,  a  value  assumed  for  the  waters  in  this  study. 

The  only  unknown  variable  remaining  is  ~a^.  This  value  varies  with  phytoplankton 
condition  and  the  light  environment  (Kishino  et  al.  1986).  If  only  remote-sensing  information 
is  available,  must  be  estimated  from  other  data,  or  an  empirical  value  has  to  be  picked  for 
In  many  calculations,  an  value  of  ~  0.015  (m^(mg  chi  c)*’)  was  measured  (Marra  et 
al.  1992)  or  has  been  used  (Bannister  1974;  Dubinsky  et  al.  1984;  Smith  et  al.  1989)  and 
assumed  to  be  vertically  constant. 

As  we  know,  [C]  using  the  CZCS  algorithm  is  accurate  to  about  a  factor  of  2  (Gordon 
and  Morel  1983),  and  can  vary  by  a  factor  of  4  (Morel  and  Bricaud  1981;  Kishino  et  al. 
1986;  Bricaud  et  al.  1988;  Laws  et  al.  1990;  Carder  et  al.  1991)  because  of  the  pigment 
composition,  the  "package  effect,"  and  the  color  of  the  light  field.  On  a  global  scale,  the 
combined  variation  of  [Q  and  for  APP  could  potentially  range  over  a  factor  of  8  when  other 
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terms  are  known. 

Actually  Eq.  1  is  a  simplified  version  of  a  complete  expression  as  noted  by 
Sathyendranath  et  al.  (1989),  Morel  (1991),  Platt  et  al.  (1991)  and  Bidigare  et  al.  (1992);  i.e. 

PPiz)  -  Hz)f^a^a)E^(z,k)dX  ,  (10) 

and  the  absorbed  photons  by  phytoplankton  are 

APP^z)  -  f™a^(X)E„{z.X)dX  .  (II) 

However,  the  a^(K)  curves  in  those  full  spectral  models  were  directly  dependent  upon  fl^’(X)  and 
chlorophyll  a  concentration,  which  means  the  determination  of  [C\  and/or  [chi  d\  is  still  the  first 
step  in  the  PP  calculation,  and  empirical  curves  for  fl/(X)  must  be  selected.  Due  to  uncertainties 
in  both  [C]  and  a^*(X)  estimates  using  remote  sensing,  if  a^(\)  and  a(X)  can  be  obtained  directly 
and  analytically  from  remotely  measured  signals  (e.g.  remote-sensing  reflectance),  the  problems 
involved  in  the  estimation  of  [C]  and  choosing  values  for  p  and  a/  will  be  avoided.  The 
accuracy  of  estimating  APP  and  PP  can  then  be  improved  when  E„(0),  <)>„,  PAR^  and  n  are 
certain. 

For  a  study  carried  out  in  the  North  Atlantic,  where  both  primary  production  and  above¬ 
surface  remote-sensing  reflectance  were  measured,  two  sets  of  primary  production  calculations 
were  made.  One  used  Eqs.  6  and  7  (hereafter  referred  to  as  the  "pigment"  method),  and  the 
other  used  Eqs.  6  and  11  (hereafter  referred  to  as  the  "absorption"  method).  Comparisons  of  the 
calculated  to  the  measured  daily  primary  production  were  then  made,  using  identical  input 
parameters  for  the  light  field  and  the  PP  model. 
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The  new  approach 


For  remote-sensing  reflectance  R„  (ratio  of  the  water-leaving  radiance  to  the  downwelling 
irradiance  above  the  surface)  at  N  wavelengths,  recent  studies  (Gordon  et  al,  1988;  Morel  and 
Gentili  1993;  Lee  et  al.  1994a;  Lee  1994)  show  that: 


RrsiK)  ~~ 


0.17 

3.4 

0.17 

jif) 

3.4 

(12) 


where  a„  and  are  the  absorption  and  backscattering  coefficients  of  pure  sea  water, 
respectively,  and  can  be  found  in  Smith  and  Baker  (1981);  is  the  absorption  coefficient  of 
detritus  and  gelbstoff,  and  can  be  expressed  as  (Roesler  et  al.  1989;  Carder  et  al.  1991): 


.  <‘3) 

X  and  Y  describe  the  scattering  effects  of  suspended  particles. 

Eq.  12  is  a  series  of  N  equations,  which  in  total  have  at  least  N  -1*  4  unknowns  (N  for 
a^,  2  for  ajg  and  2  for  particle  scattering)  given  only  R„.  For  the  N  equations,  unless  we  can 
dramatically  reduce  the  unknowns  regarding  a^,  there  will  be  no  certain  solutions  for  a^  and  the 
total  absorption  a. 

Bidigare  et  al.  (1990)  suggested  that  can  be  reconstructed  by  knowing  concentrations 
of  the  pigments  and  the  specific  absorption  for  each  pigment,  but  these  analyses  cannot  be 
performed  remotely.  Hoepffner  and  Sathyendranath  (1991)  suggest  that  the  pigment  spectrum 
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can  be  simulated  by  a  sum  of  1 1  Gaussian  bands.  For  these  1 1  Gaussian  bands,  their  center 
wavelengths  and  half  bandwidths  vary  from  species  to  species.  Even  if  the  center  wavelengths 
and  half  bandwidths  are  fixed,  we  still  need  to  derive  11  parameters  to  obtain  a^.  Morel  (1980) 
suggested  the  use  of  averaged  chlorophyll-specific  absorption  coefficient  while  Roesler  and  Perry 
(1994)  suggested  the  use  of  a  single,  averaged,  spectral  shape  for  phytoplankton  absorption. 
Using  these  approaches,  only  one  unknown  (the  concentration  of  pigments  or  a  scale  factor)  is 
needed  to  model  a^.  However,  due  to  the  "package  effect"  and  variable  pigment  composition, 
it  is  well  known  that  the  shape  of  curve  can  vary  widely  from  region  to  region,  and  no  single 
curve  shape  can  be  used  globally.  So,  a  simpler  expression  with  adequate  accuracy  for 
describing  the  different  shape  of  curves  is  necessary. 

By  analyzing  surface  curves  collected  for  waters  from  the  Gulf  of  Mexico  with  a 
pigment  range  of  0.07  mg/m^  ^  [chi  a]  ^  40  mg/rrf,  a  mathematical  expression  was  derived 
from  a  combination  of  3  simple  functions  (Lee  et  al.  1994b),  which  involves  6  parameters. 
Figure  1  shows  examples  of  the  a^(X)  simulation.  Among  the  6  parameters,  2  parameters  vary 
only  slightly  for  different  waters,  and  only  2  parameters  have  strong  effects  on  APP  (Eq.  1 1 ) 
For  the  wavelength  range  of  400  nm  ^  X  <  7(X)  nm,  the  expression  is  (Lee  1994;  Lee  et  .il. 
1994b): 


-Wlailil 

400  <  X  <  590,  a^W  -  '  ■"  ■ 


(N) 


aA656)  -  aA590) 

590  <  X  <  656,  a.(X)  =  aA590)  +  ^ - (X  -  590),  0^) 

’  ♦  656  -  590 


and 


Parameter  F  describes  the  width  of  the  curvature  from  400  nm  to  ~560  nm,  X,  is  the 
reference  needed  to  obtain  the  blue  peak  wavelength,  100  +  X;,  for  a  log-normal  curve,  Xj  is 
the  red  peak  wavelength  of  a  Gaussian  curve,  and  2.3550^  determines  the  half  width  of  the  band 
around  the  red  peak.  For  the  examples  given  in  their  report,  F  varied  from  1.6  to  4.2,  X/  from 
338  -  342  nm  with  80%  at  340  nm,  Xj  ranged  from  672  to  675  nm,  and  2.355(T2  varied  from  2 1 

to  35  nm.  a^i  varied  widely  from  0.01  to  0.83  m  *  while  ranged  from  0.21  to  0.85.  The 

root-mean-square  (rms)  difference  between  the  measured  and  simulated  a^,  values  was  11%.  This 
difference  dropped  to  2%  when  comparing  the  calculated  .(4PP  using  measured  versus  simulated 
curves  (Lee  et  al.  1994b),  meaning  the  expression  works  well  for  the  calculation  of 
phytoplankton-  absorbed  photons.  Using  this  expression  in  the  inversion  to  analytically  derive 
and  a  allows  calculation  of  PP  at  any  depth  simply  from  R„  when  <l>„,  PAR^,  v  and  EJiSi)  are 
determined  for  a  euphotic  layer  that  is  well-mixed.  This  is  often  true  for  high-latitude  waters  and 
upwelling  regions. 

By  analyzing  the  measured  a^(440)  and  the  simulation  parameters,  Lee  (1994)  found  that 
for  a^(440)  in  the  range  of  0.01  to  1.0  m  *, 
a^2^a^i  *  -86  +  .16  ln(a^|),  normalized  rms  error:  17.2%, 

F  «  2.89  exp[-.505  tanh[.56  ln(a^/.043)]],  normalized  rms  error:  12.4%, 

02  »  14.17  -f-  .9  ln(a^;),  normalized  rms  error:  5.6%. 

Thus  whenever  a^i  is  determined,  the  full  spectrum  can  be  estimated  with  the  change 
of  the  curve  shape  considered,  at  least  to  first  order.  Instead  of  attempting  to  statistically  relate 
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a^i  with  ratios  of  R„,  it  is  possible  to  analytically  derive  it  if  N  is  equal  to  or  greater  than  5.  So, 
the  absorption  coefficient  of  pigments  as  well  as  the  total  absorption  coefficient  can  be  obtained 
when  there  are  R„  values  at  more  than  5  wavelengths. 

After  the  retrieval  of  and  a  from  measured  R„,  scalar  irradiance  at  depth  z  can  be 
calculated  (Eq.  4)  with  *  1.08  D^O)  a  (Gordon  1989),  where  Z>^0)  is  the  sub-surface 
distribution  function.  £^(0)  can  be  calculated  using  modifications  to  the  models  of  Gregg  and 
Carder  (1990),  Bishop  and  Rossow  (1991),  or  using  in  situ  mooring  or  drifter  PAR  data.  Thus, 
when  <f>„,  PAR^  and  u  are  known,  PP  for  any  depth  can  be  calculated  through  Eqs.  6  and  11. 

For  data  collected  from  the  Marine  Light-Mixed  Layer  (ML-ML)  study  (2rW/59°N)  of 
May  1991,  two  sets  of  calculations  were  made.  The  traditional  set  is  derived  using  Eqs.  6,  7, 
8  and  9  (the  "pigment"  method)  as  discussed  in  the  introduction.  The  other  is  derived  by  Eqs. 
6  and  11  (the  "absorption"  method),  where  the  absorption  coefficients  of  the  pigments  and  the 
total  absorption  coefficients  are  analytically  derived  just  from  measured  remote-sensing 
reflectance. 

In  both  calculations,  measured  PAR(0)  values  for  each  station  were  used  for  the  PAR{z) 
calculations  and  comparison  of  the  two  methods  because  most  of  the  stations  were  taken  during 
cloudy  weather.  27^(0)  is  approximated  as  1.2  (Platt  et  al.  1991)  in  the  "absorption"  method  as 
most  of  the  station  days  were  cloudy.  The  calculated  results  for  each  day  at  each  depth  were 
compared  with  the  measured  primary  production  values,  and  the  water-column-integral 
production. 
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Data  and  measurements 

Data  for  PAR{z),  PP{z),  and  RJk)  were  collected  on  an  ML-ML  cruise  from  May  17  to 
May  24,  1991  in  the  waters  south  of  Iceland  (21“W/59“N).  PAR  and  PP  measurements  were 
carried  out  on  floating  array  for  each  of  four  days:  May  17,  20,  22  and  24.  The  sampling  site, 
conditions,  wind,  mixing,  and  nutrients  are  presented  in  Marra  et  al.  (1994)  and  Plueddemahn 
et  al.  (1994).  In  general,  it  was  windy  for  the  stations,  and  the  euphoric  zone  was  well  mixed. 

Shipboard  measurements  -  PAR^z)  at  four  depths  (0,  2,  12.5  and  25  m)  was  monitored 
and  averaged  for  every  10  minute  interval  throughout  the  day  using  a  Biospherical  PAR  sensor 
attached  to  each  incubation  array  (Marra  et  al.  1994). 

Dawn-to-dusk  incubations  (17  hours)  with  four  replicates  were  carried  out  in  situ  at  each 
of  six  depths  (5,  10,  15,  20,  30  and  40  m)  chosen  to  span  the  euphoric  zone.  Primary  production 
(PF)  measurements  were  made  using  the  ‘‘‘C  technique  (Marra  et  al.  1994). 

Remote-sensing  reflectance,  R„(K)  -  Hyperspectral  remote-sensing  reflectance  R„(K)  was 
measured  for  the  noon  site  by  the  method  developed  by  (Carder  and  Steward  1985).  The 
instrument  was  the  Spectron  Engineering  spectral  radiometer  (model  SE-590).  It  has  253 
channels  and  covers  the  wavelength  range  from  370  -  1100  nm  with  a  spectral  resolution  of  7 
nm  using  a  200  ^m  entrance  slit.  With  this  instrument  the  water-leaving  radiance  and 
downwelling  sky  radiance  were  directly  measured,  downwelling  irradiance  was  measured  by 
viewing  a  standard  diffuse  reflector  (Spectralon,  ~  8.5%).  Reflected  sky  radiance  was  corrected 
by  the  method  of  Carder  and  Steward  (1985)  when  calculating  i?„(X). 
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Results  and  discussion 

Comparisons  of  calculated  and  measured  PAR  and  calculated  and  measured  PP  are 
presented  in  Figures  2a  -  2d  and  Figures  3a  -  3d,  respectively.  All  data  are  summarized  in  Table 
2  and  all  PP  data  in  Figure  4a,  with  Figure  4b  shows  the  results  when  v  is  0  (i.e.  no 
photoinhibition  considered). 

It  can  be  seen  that  the  "absorption"  method  yields  close  estimates  of  PAR{z),  while 

values  calculated  by  the  "pigment"  method  were  consistently  higher,  particularly  at 
depth.  One  reason  for  this  is  that  the  CZCS  pigment  algorithm  estimated  [chi  d\  as  much  as  a 
factor  of  5  lower  than  the  measured  surface  values  for  these  waters.  This  had  the  consequence 
that  Kj  values  (Eq.  9)  calculated  from  these  [C]  ([chi  fl]/0.75)  were  small.  This  difference  could 
be  the  result  of  an  incorrect  CZCS  algorithm  for  that  environment  (e.g.  Balch  et  al.  1989; 
Mitchell  and  Holm-Hansen  1991)  or  less  likely  that  there  were  substantial  errors  or  discrepancies 
in  the  measurements  of  R„  or  "sea  truth"  [chi  a].  Also,  the  empirical  relationship  between 
and  [C]  might  not  hold  for  these  waters.  However,  the  same  field  data  are  used  for  both 
methods,  so  models  and  algorithms  will  be  largely  responsible  for  differences  between  the 
methods. 

Calculated  primary  production  values  by  the  "absorption"  method  were  highly  comparable 
to  the  measured  values,  whereas  the  "pigment"  method  significantly  underestimated  PP, 
particularly  in  the  surface  waters.  The  average  difference  of  the  water-column-integrated  PP  is 
20%  between  the  measured  and  the  "absorption"  method,  while  it  is  61%  between  the  measured 
and  the  "pigment"  method  (Table  2).  The  P-  is  0.95  (n=24)  for  a  linear  regression  between  the 
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values  calculated  by  the  "absorption"  method  and  that  of  the  measurements,  with  a  slope  of  1.26 
and  +  32%  difference.  The  r^  is  0.85  (n=24)  between  the  values  calculated  by  the  "pigment" 
method  and  that  of  the  measurements,  with  a  slope  of  0.34  and  -  78%  difference  (see  Figure  4). 
These  results  indicate  that  there  is  about  a  factor  of  3  improvement  in  the  accuracy  of  PP 
calculation  from  remote  sensing  for  ML-ML  waters  using  the  "absorption"  method  suggested 
here. 

For  May  17  (Figure  3a),  the  calculated  PP  values  by  the  "absorption"  method  in  the 
euphotic  zone  were  generally  greater  than  the  measured  values.  Factors  that  can  account  for 
these  differences  include  possible  errors  in  measurement  and  R„  inversion,  grazing  effects  on 
the  measured  values,  and  either  the  photoinhibition  factor  used  was  smaller  than  the  "real" 
situation,  or  and/or  PAR^  used  was  higher  than  the  "actual"  values. 

For  May  20  and  May  24  (Figures  3b  and  3d),  the  PP  values  calculated  by  the 
"absorption"  method  were  very  close  to  those  measured,  which  suggest  that  our  assumed  values 
of  4>„,  PAR^  and  v  for  these  two  days  were  probably  close  to  reality. 

The  greatest  differences  between  measured  and  calculated  values  from  both  the  "pigment " 
and  "absorption"  methods  occurred  on  the  only  sunny  day  of  the  cruise,  May  22  (Figure  3c  i. 
where  measured  surface  layer  production  was  smaller  than  the  second  layer.  Other  than 
measurement  errors,  this  is  an  indication  of  increased  photoinhibition,  as  suggested  by  the 
calculated  values  of  the  "absorption"  method,  which  was  also  discussed  in  Carder  et  al.  (1994) 
and  Marra  et  al.  (1994).  However,  the  calculated  PP  values  by  the  "absorption"  method  are 
greater  than  the  measured  values  for  the  first  four  depths  (~40%).  Possible  reasons  for  this 
difference  might  be  an  overestimation  of  <f>„  and/or  an  underestimation  of  p  for  the  phytoplankton 
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in  this  environment,  since  4>„  and/or  v  might  be  a  function  of  light  history.  The  other  3  cruise 
days  were  overcast  with  an  expected  adaptation  of  the  plants  to  a  low  light  level.  The  sudden 
exposure  to  the  bright  light  probably  caused  extreme  photoinhibition  in  the  dark-adapted 
phytoplankton  population  with  the  resultant  overestimation  of  PP  by  the  "absorption"  method. 

The  PP  values  calculated  by  the  "pigment"  method  for  the  sunny  day,  however,  do  not 
show  a  photoinhibition  response  in  the  surface  waters  (Figure  3c).  One  reason  for  this  is  that 
not  only  varies  with  the  chlorophyll-specific  absorption  coefficient,  it  also  varies  with  the 
light  field  (Eq.  2).  An  a/  value  for  surface  water  is  not  necessarily  appropriate  for  the  whole 
water  column,  even  for  a  well-mixed  ocean  (Kishino  et  al.  1986).  So,  the  'a^,*  value  is  one  of 
the  major  sources  of  error  in  PP  calculations  by  the  traditional  method. 

For  the  waters  studied,  PP  values  calculated  by  the  "pigment"  method  were  lower  than 
measured  rates  by  a  factor  of  3.  Since  PP  values  calculated  by  the  "absorption"  method  were 

close  to  measured  values,  our  value  of  <{>„,  0.074  mol  C  (Ein  absorbed)  ',  is  apparently  close  to 

\ 

the  actual  value.  Therefore,  the  differences  between  calculated  and  measured  PP  for  the 
"pigment"  method  are  largely  due  to  the  estimation  of  [C]  and  the  product  of  p[C]a^*.  However, 
when  only  remote-sensing  data  are  available,  it  is  difficult  to  know  which  are  the  "correct"  u,* 
and  p[C]  values,  as  they  depend  upon  empirical  values.  What  makes  things  interesting  is  that  a 
factor  of  2  increase  in  both  [C]  and  will  make  calculated  PAR  and  PP  both  close  to  the 
measured  values.  The  pigment  data  in  Table  2,  however,  suggest  that  the  CZCS  algorithm 
accounts  for  most  of  the  error.  This  is  consistent  with  the  conclusion  of  Platt  et  al.  (1988)  that 
determination  of  biomass  by  remote  sensing  dominates  the  error  in  primary  production 
estimation.  So,  it  is  preferable  to  shift  the  primary  production  models  from  pigment- 
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concentration  base  to  pigment-absorption  base. 

Comparing  Eqs.  7  and  11,  the  main  difference  between  the  "absorption"  and  "pigment" 
methods  is  how  APP{z),  the  absorbed  energy  by  phytoplankton  pigments,  is  obtained  from 
remote-sensing  data.  In  the  traditional  method,  calculation  of  APP  depends  heavily  on  4 
numbers:  p.  A,  B,  and  a^.  When  there  are  only  remote-sensing  signals,  it  is  difficult  to  verify 
which  values  should  be  used  for  each  of  these  parameters,  although  each  of  them  could  be  tuned 
empirically  to  specific  regions  and  specific  seasons  (Platt  et  al.  1991).  Since  there  are  errors 
associated  with  each  number,  the  cumulative  error  in  APP  could  then  be  very  high,  even  if  R-t/? 
is  correct. 

In  the  APP  calculation  by  the  "absorption"  method,  however,  the  single  most  important 
parameter  is  analytically  derived  from  measured  R„.  Thus,  most  of  the  error  comes  from 
the  R„  inversion  process.  It  is  believed  that  absorption  coefficients  derived  from  R„  inversion 
have  an  accuracy  better  than  50%  for  "case  1"  and  "case  2"  waters  (Lee  1994).  This  means  that 
the  accuracy  in  the  APP  calculation  by  the  "absorption"  method  could  be  within  50%.  Values 
derived  using  the  traditional  method  were  only  1/3  of  the  measured  values,  however.  In  the 
"absorption"  method,  the  mid-step,  i.e.  deriving  [chi  a]  and  estimating  a  value  for  a^\  is 
eliminated.  The  result  is  highly  improved  accuracy  in  deriving  APP  and  PP. 

How  to  remotely  obtain  accurate  estimations  of  the  physiological  parameters  such  as  <t>„, 
PAR^  and  v,  however,  remains  a  challenge  for  the  calculation  of  PP  through  remote  sensing 
(Balch  et  al.  1992).  Including  light  history  and  surface  temperature  anomalies  (e.g.  see 
Kamykowski  and  Zentara  1986)  in  future  PP  models  may  reduce  this  deficiency.  Among  the 
three  parameters,  might  be  the  most  important  for  a  PP  model,  since  PP  is  proportional  to 


18 


the  value  of  <t>„.  Only  when  PAR(z)  is  very  high  is  there  a  strong  influence  on  PP  due  to  PAR^ 
and  If.  For  example,  for  p  =  0.01  (Ein/mVday)'‘,  is  0.90,  0.82,  0.67,  0.55  and  0.45  for 
PAR(z)  values  of  10,  20,  40,  60  and  80  Ein/mVday,  respectively.  This  means  that  by  dropping 
the  photoinhibition  term,  the  calculated  PP  will  increase  10%  to  about  55%.  Such  an  increase 
could  be  compensated  by  using  a  smaller  value,  but,  without  the  photoinhibition  term,  the 
increase  then  decrease  vertical  profile  of  PP  with  depth  would  not  be  simulated.  Without 
consideration  of  photoinhibition  for  the  waters  studied,  r^  decreased  from  a  value  of  0.95  to  0.93 
(n=24),  and  the  slope  increased  from  1.26  to  1.39  (see  Figure  4b).  This  suggests  that  the 
photoinhibiton  term  is  important  for  primary  production  model  applications  to  the  ML-ML  site. 


Conclusions  and  expectations 
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1.  Based  on  the  results  presented  here,  the  a^  simulation  works  very  well  in  modeling 
PAR(z)  based  on  PAR(0)  and  R„.  Combined  with  parameters  regarding  phytoplankton 
photosynthesis,  the  PP  values  calculated  by  the  "absorption"  method  were  close  to  the  measured 
ones  with  r^  =  0.95,  a  slope  =  1.26,  and  average  difference  of  32%.  This  means  a  great 
improvement  in  PP  estimation  accuracy,  where  the  modeled  values  were  three-fold  smaller  than 
the  measured  ones  using  the  traditional  method.  So,  it  is  preferable  to  shift  the  primary 
production  models  from  pigment-concentration  base  to  pigment-absorption  base.  These  results 
also  indicate  that  the  combined  PP  model. 


PP(z)  -  APP(z) 


PAR^^PARQl)^ 


-vPAKd 

f 
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works  fine  for  the  waters  studied. 

2.  It  is  not  necessary  to  know  [chi  d\  for  the  calculation  of  PP.  In  the  traditional  method, 
the  estimation  of  [C]  and  a^  separately  is  an  inherent  disadvantage  in  calculating  PP  based  upon 
satellite  or  aircraft  remote-sensing  measurements  of  ocean  color.  What  is  really  needed  for 
models  of  PP  is  the  absorbed  energy  by  phytoplankton,  and  the  physiological  parameters  of 
phytoplankton  photosynthesis. 

3.  Coupled  investigations  between  PP  and  R„  measurements  need  be  carried  out  widely 
to  generate  more  data  to  test  and  improve  the  method. 

4.  Cloud-covered  marine  environments  such  as  the  ML-ML  site  in  the  spring  cannot  often 
be  viewed  from  space.  They  can,  however,  be  measured  using  drifters  or  moorings  or  from 
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aircraft  flying  beneath  the  cloud  base.  Autonomously  piloted  aircraft  (Montgomery  et  al.  1994) 
are  becoming  available  that  could  provide  spectral  observations  of  the  ocean  in  the  future  to 
increase  our  knowledge  of  cloud-fast  environments.  Methods  such  as  developed  in  this 
contribution  can  provide  improved  estimates  of  primary  production  and  ultimately  pigments  if 
combined  with  appropriate  drifters,  moorings,  or  aircraft  and  sensor  packages. 

5.  For  space-based  measurements,  methods  must  be  pursued  to  remotely  estimate  the 
physiological  parameters,  perhaps  by  their  covariance  with  some  remotely  measured  variable 
such  as  sea  surface  temperature  anomalies,  wind-stress  history,  and  light  history. 
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Table  1:  Symbols  and  units. 

a(\)  Total  absorption  coefficient,  m  ';  a(\)  =  a^(\)  +  +  a^(X) 

a^g(X)  Absorption  coefficient  of  detritus  and  gelbstoff,  m  ‘ 

a^(X)  Absorption  coefficient  of  phytoplankton  pigments,  m'‘ 

a^’(X)  Pigment-specific  absorption  coefficient,  m^  (mg  chi  a)  ' 

a*(440)  and  a^(675),  respectively,  m  * 
a,„(X)  Absorption  coefficient  of  pure  water,  m'* 

APP  Absorbed  photons  by  phytoplankton,  Ein/mVday 

bh(K)  Backscattering  coefficient,  m'* 

[chi  a]  chlorophyll  a  concentration,  mg/m^ 

[C]  chlorophyll-like  pigment  concentration,  mg/m^ 

D/0)  Subsurface  downwelling  distribution  function 

Eg  Quantum  scalar  irradiance,  quanta/mVs 

Kj  Diffuse  attenuation  coefficient  for  downwelling  irradiance,  m'* 

Water-leaving  radiance,  W/m^ 

PAR  Photosynthetically  available  radiation  (integrated  from  4(X)  -  700  nm), 

Hn/mVday 

Pr4if(0)  Surface  PAR,  Ein/mVday 

PAR^  Value  of  PAR  where  <t>  =  0„/2,  Ein/mVday 

PP  Primary  production,  /xmol  C/1 

R„  Remote-sensing  reflectance,  =  LJE^,  sr* 


Quantum  yield,  mol  C  (Ein  absorbed)  ' 

Maximum  quantum  yield,  mol  C  (Ein  absorbed)  ' 
Wavelength,  nm 
fraction  of  [chi  a]  in  [C] 


Photoinhibition  factor,  (Ein  m'^day  ')  ' 


Table  2.  Result  summary. 
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Day 

May  17 

May  20 

May  22 

May  24 

P/i/?(0)  (Ein/mVday) 

38.27 

16.25 

65.72 

28.73 

^„(443)//f„(550) 

1.2 

1.6 

1.8 

1.8 

surface  [chi  a]  (mg/m*) 

2.9 
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1.5 

1.0 

CZCS  [chi  a]  (mg/m*) 
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0.34 

Mea.  integral  PP 
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Figure  legends 


Figure  la  -  lb.  Examples  for  the  a^(X)  simulation. 

Figure  2a  -  2d.  Depth  profiles  of  PAR  for  May  17,  May  20,  May  22  and  May  24, 

respectively. 

Open  triangle-down:  measured.  Filled  square:  estimated  by  the  "absorption"  method. 
Open  circle:  estimated  by  the  "pigment"  method. 

Figure  3a  -  3d.  Depth  profiles  of  PP  for  May  17,  May  20,  May  22  and  May  24, 

respectively. 

Symbols  as  in  Figure  2. 

Figure  4a  -  4b.  Comparison  of  all  measured  and  calculated  PP. 

In  4a,  photoinhibition  was  considered;  in  4b,  photoinhibition  was  not  considered 
(details  see  text). 

Symbols  as  in  Figure  2. 
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Abstract .  Using  bio-optical  estimates  of  gelb- 
stoff  and  a  few  platinum  measurements  of  dissolved 
organic  carbon  (DOCp^)  ,  a  budget  of  the  meridional 
flux  of  DOC  and  dissolved  organic  nitrogen  (DON) 
across  36°25'N  in  the  North  Atlantic  is  constructed 
from  previous  inverse  models  of  water  and  element 
transport.  Distinct  southward  subsurface  fluxes  of 
dissolved  organic  matter  (DOM)  within  subducted 
shelf  water,  cabelled  slope  water,  and  overturned 
basin  water  are  inferred.  Within  two  cases  of  a 
positive  gradient  of  DOCp^  between  terrestrial/ 
shelf  and  offshore  stocks,  the  net  equatorward 
exports  of  O2  and  DOCp^  from  the  northern  North 
Atlantic  yield  molar  ratios  of  2.1  to  9.1,  compared 
to  the  expected  Redfield  O2/C  ratio  of  1.3.  In  the 
first  case,  63%  of  the  apparent  oxygen  utilization 
demands  of  the  water  column  may  be  met  by  DOC, 
instead  of  only  14%  in  the  second  scenario,  pre¬ 
serving  a  role  for  falling  particles  in  the  sea. 
With  a  DOC/DON  ratio  of  10,  the  larger  net  south¬ 
ward  export  of  DON  across  36°25'N  balances  the 
postulated  net  northward  input  of  1.7  x  10^  kg  NO3 
s'^  of  unutilized  nitrate  within  the  Gulf  Stream. 
Without  an  enhanced  supply  of  DOM  from  the  shelves, 
a  zero  seaward  gradient  of  DOM  in  the  third  case 
suggests  that  none  of  the  poleward  nitrate  flux  is 
returned  southward  as  DON,  but  instead  a  net  pole- 
ward  flux  of  DON  prevails  as  well .  Our  present 
estimates  are  confounded,  however,  by  the  seasonal 
and  multiyear  variability  of  sinking  processes  in 
the  North  Atlantic.  Future  active  and  passive 
remote  sensors,  field  programs,  and  simulation 
models  must  now  discriminate  between  particulate 
and  dissolved  components  of  surface  color  signals 
to  verify  the  importance  of  both  continental 
margins  and  DOM  in  global  biogeochemical  cycles. 

1 .  Introduction 

Recent  measurements  of  dissolved  organic  carbon 
(DOC)  and  nitrogen  (DON)  in  sea  water,  involving 
varying  amounts  of  a  platinum  catalyst  [Suzuki  et 
al.,  1985,  1992;  Suglmura  and  Suzuki,  1988]  have 
rekindled  controversies  [Williams  and  Druffel, 
1988;  Jackson,  1988],  ignited  at  the  turn  of  the 
century  [Putter,  1909] ,  about  the  roles  of 
dissolved  organic  matter  (DOM)  in  the  ocean.  Prior 
large  estimates  of  DOC  in  the  ocean,  obtained  with 
different  methods  [Plunkett  and  Rakes tr aw ,  1955; 
Skopintsev  et  al.,  1966;  Gordon  and  Sutcliffe, 
1973;  Melnikov  and  Pavlov,  1978]  had  been  rejected 
in  favor  of  the  lower  values  obtained  with  ultra¬ 
violet  and  peroxodisulfuric  acid  techniques  of  wet 
combustion  [Armstrong  et  al.,  1966;  Menzel  and 
Vaccaro,  1964] ,  although  unmeasured  organic  acids 
could  account  for  differences  in  potentlometric  and 
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manometric  estimates  of  total  CO2  in  seawater 
[Bradshaw  and  Brewer,  1988] .  If  the  even  larger 
values  of  DOC  measured  with  the  platinum  technique 
are  correct,  they  raise  questions  about  prior 
estimates  of  regional  sources  and  sinks  of  the 
global  carbon  budget  [Tans  et  al.,  1990],  particu¬ 
larly  in  the  North  Atlantic,  which  is  the  subject 
of  our  analysis. 

Using  depth  profiles  of  Apparent  Oxygen  Utili¬ 
zation  (AOU)  and  the  platinum  estimates  of  DOC,  80% 
[Toggweiler,  1989]  to  100%  [Sugimura  and  Suzuki, 

1988]  of  the  respiration  in  the  aphotic  zone  is  now 
attributed  to  consumption  of  DOC,  rather  than  of 
falling  particles  in  the  western  Pacific.  This 
hypothesis  is  a  sharp  contrast  to  a  previous 
estimate  for  the  same  region  [Ogura,  1970] ,  which 
suggested  that  at  most  25%  of  AOU  demands  were  met 
by  a  DOC  source,  based  on  a  Redfield  molar  O2/C 
ratio  of  1.3  [Redfield  et  al.,  1963]  and  the  wet 
combustion  method  [Menzel  and  Vaccaro,  1964] . 
These  new  findings  imply  that  the  global  downward 
flux  of  DOC  may  be  equivalent  to  the  particulate 
rain  of  marine  carbon,  representing  50-67%  of  new 
production  in  simple  models  of  the  open  ocean 
[Toggweiler,  1989;  Bacastow  and  Maier-Reimer ,  1991; 
Najjar  et  al.,  1992], 

Local  invasion  of  most  of  the  deep  sea  by 
surface  DOM  must  occur  at  slower  rates  than  the 
settling  velocities  (-100  m  d"^)  of  large  particles 
[Alldredge  and  Gotschalk,  1989] ,  however.  At 
similar  rates  to  those  of  nitrate  infusion  up  the 
oceanic  nutricline  [Walsh,  1991] ,  for  example, 
diffusive  invasion  by  surface  DOM  at  a  rate  of 
1  cm^  s"^  is  equivalent  to  a  penetration  speed  of 
-37  m  yr"^  over  the  upper  100  m  of  the  water 
column.  At  this  rate,  it  would  take  at  least  27 
years  to  arrive  at  a  depth  of  1000  m;  eddy  diffu¬ 
sion  of  DOM  across  the  main  pycnocline  would  be 
even  slower.  Over  this  time,  lateral  advection, 
particularly  that  of  western  boundary  currents 
(>100  km  d'^)  ,  may  allow  pilgrimage  of  tropical  and 
subtropical  DOC  pools  to  polar  latitudes 
[Toggweiler ,  1989] ,  where  convective  mixing  may 
accelerate  the  descent  of  DOC. 

Otherwise,  photodegradation  [Keiber  et  al., 

1990]  and  bacterial  utilization  [Kirchman  et  al., 

1991]  of  DOC  in  surface  waters  may  return  most  of 
this  dissolved  form  of  photosynthate  back  to  the 
atmosphere  as  biogenic  CO2,  leaving  behind  a 
relatively  inert  component  of  DOC.  The  '^''C  age 
[Williams  and  Druffel,  1987]  and  C/N  ratio  [Jackson 
and  Williams ,  1985]  of  deep-sea  pools  of  DOM  deter¬ 
mined  by  UV-radiation  are  6000  years  and  25,  for 
example,  instead  of  the  C/N  ratios  of  5-7  [Sugimura 
and  Suzuki,  1988]  and  a  bomb-label  [Druffel  et  al.  , 

1989]  for  DOC  measured  with  the  platinum  method. 
The  more  recent  observations  imply  a  younger  and 
less  refractory  component  of  DOC  as  well,  at  least 
in  offshore  waters  [Kirchman  et  al.,  1991;  Suzuki 
and  Tanoue,  1991;  Benner  et  al.,  1992].  In  con¬ 
trast,  the  radiocarbon  age  of  near- surface  humic 
and  fulvic  acids,  comprising  20-50%  of  the  DOC 
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sampled  by  UV-oxidation  [Harvey  et  al.,  1983; 
Meyers -Schulte  and  Hedges,  1986;  Williams  and 
Draff  el,  1987],  is  about  4000  years  [Druffel  et 
al.,  1989] . 

The  identity  of  these  additional  labile  pools  of 
DOM,  if  and  whether  they  really  exist  in  soluble, 
colloidal  (0.005-0.050  ftm  diameter),  or  fine 
particle  (0.050-1.000  gn)  form,  remains  an  unre¬ 
solved  issue  [Benner  et  al.,  1992].  Their  sources 
(terrestrial,  coastal,  and  oceanic)  and  fates 
(photodegradation,  remineralization,  and  storage  as 
inert  forms  of  DOC)  are  also  open  questions,  which 
underscore  our  inadequate  understanding  of 
temporal -spatial  variance  in  the  biogeochemical 
cycles  of  carbon  and  nitrogen  [Walsh,  1991] . 

To  qualify  as  a  dynamic  storage  reservoir  in 
global  carbon  budgets  [Tans  et  al. ,  1990] ,  however, 
DOM  stocks,  like  those  of  dissolved  inorganic 
carbon  (DIG) ,  must  be  removed  from  surface  regions 
of  the  ocean  by  sinking  water,  if  not  by  falling 
particles ,  to  avoid  respiration  back  to  CO2  at  the 
sea-surface.  We  explore  this  mechanism  by  consid¬ 
ering  (1)  the  possible  DOC  stocks  in  shelf,  slope, 
and  basin  regions  of  the  North  Atlantic  Ocean  as 
estimated  with  both  platinum  (DOCp^)  and  ultravio¬ 
let  (D0C„„)  techniques;  (2)  their  potential  annual 
sinking  fluxes  in  relation  to  freshwater  and 
photosynthetic  supplies;  and  (3)  the  implied  net 
southward  transfer  of  subsurface  DOC  and  DON  in 
relation  to  other  meridional  fluxes  of  oxygen  and 
nitrate  across  36®25'N. 

2 .  Methods 

In  a  previous  study,  Rintoul  and  Wunsch  [1991] 
suggested  that  the  net  poleward  transfer  of  nitrate 
by  the  Gulf  Stream  System  across  36°25'N  (the 
location  of  this  section  is  shown  in  Figure  1) 
might  be  balanced  by  a  southward  transport  of  DON. 
Rintoul  [1988]  provides  zonally  averaged  estimates 
of  meridional  transport  of  water  over  17  vertical 
levels  (Table  1)  at  this  latitude,  based  on  an 
inverse  model,  which  conserves  both  mass  and  sili¬ 
cate  below  a  depth  of  750  m  and  has  a  reference 
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Fig.  1.  The  distribution  of  salinity  in  the  North 
Atlantic  at  a  depth  of  30  m  (after  Worthington 
[1976])  in  relation  to  (1)  a  trans -Atlantic 
hydrographic  section  at  36°25'N,  (2)  observations 
of  DOCpt,  at  46°30'N  (point  B)  and  59°30'N  (point  A), 
and  (3)  bio -optical  estimates  of  Gelbstoff  between 
Iceland  and  Nova  Scotia  (stations  1-9) . 


level  of  3000  dbar.  Using  their  results,  we  con¬ 
struct  budgets  of  the  possible  southward  export  of 
DOC  and  DON  from  the  North  Atlantic  at  four  of  the 
depth  Intervals,  over  which  sinking  water  may  exit 
this  basin  towards  the  equator.  The  net  exchange 
over  the  whole  water  column  is  then  obtained  by 
summing  the  DOCpt  D0C„v  transports  in  each  of 
the  17  levels  (Table  1). 

For  example,  as  a  result  of  heat  extraction  from 
the  Gulf  Stream  and  Influx  of  polar  water  within 
the  East  Greenland  Current,  a  precursor  of  North 
Atlantic  Deep  Water  (NADW)  is  formed  in  the  Iceland 
Sea  [Swift,  1984].  South  of  Denmark  Strait,  NADW 
sinks  to  depths  >4000  m  and  moves  south  as  part  of 
the  Deep  Western  Boundary  Current  (DWBC)  off  North 
America  [Broecker  et  al. ,  1991].  This  source  of 
sinking  water  is  a  likely  mechanism  for  the  export 
of  DOC  from  surface  waters  of  the  Arctic  Seas 
[Duursma,  1965],  providing  a  sequestration  time 
scale  of  thousands  of  years. 

Similarly,  mixing  of  Labrador  Current  Water 
(LCW)  with  the  western  boundary  current  to  the  east 
of  the  Grand  Banks  (note  the  packed  isopleths  of 
salinity  in  Figure  1)  results  in  about  4  Sv  of 
water  cabelling  to  depths  >2000  m  [^rauss  et  al., 
1990] ,  the  same  as  previously  estimated  for  sinking 
surface  waters  of  the  Labrador  Sea  [Sverdrup  et 
al.,  1942;  Wright,  1972].  Seasonal  convective 
overturn  of  additional  North  Atlantic  Current  (NAC) 
waters  farther  east  in  the  North  Atlantic  basin 
also  leads  to  surface  mixed  layer  depths  of  >400  m 
[Levitus,  1982]  in  the  northeast  end  of  the  sub¬ 
tropical  gyre,  where  Ekman  pumping  and  lateral 
induction  [Pedlosky,  1990;  Huang,  1990]  may  further 
sequester  DOC.  Finally,  some  of  the  southward  flow 
from  the  Greenland,  Iceland,  and  Labrador  Seas 
exits  the  Middle  Atlantic  Bight  (MAB)  at  Cape 
Hatteras  [Chapman  and  Beardsley,  1989] ,  where  the 
shelf  water  is  subducted  to  depths  >100  m  along  the 
northwest  wall  of  the  Gulf  Stream  [Lilllbridge  et 
al.,  1990]. 

The  Arctic  Ocean  is  analogous  to  the  head  of  a 
large  North  Atlantic  estuary  [Walsh  et  al.,  1989], 
where  74%  of  the  freshwater  discharge  (Figure  1) 
and  associated  terrestrial  DOC  loading  of  the 
northern  hemisphere  occur  [Walsh,  1988].  A  signif¬ 
icant  part  of  the  total  sequestration  of  DOC 
released  within  coastal,  slope,  and  basin  ecosys¬ 
tems  may  thus  take  place  along  the  northwestern 
edges  of  the  Gulf  Stream  System  in  the  North 
Atlantic.  This  could  occur  within  sinking  surface 
waters  via  NADW  formation,  LCW  cabelling,  basin 
convective  mixing,  and  shelf  water  subductlon,  part 
of  which  may  flow  south  across  36°25'N  at  respec¬ 
tive  depths  of  perhaps  3400-4300  m,  of  1600-3100  m, 
of  600-750  m,  and  of  100-350  m  (Table  1) . 

At  other  depths ,  currents  flow  north ,  of  course , 
such  that  the  net  exchange  of  water  across  36°25'N 
is  negligible  (Table  1) .  Is  the  net  exchange  of 
DOM  also  zero?  In  an  attempt  to  answer  this 
question,  we  estimate  the  possible  stocks  of  both 
DOCpt  each  depth  interval  along 
36°25'N. 

To  place  the  southward  transfer  of  DOM  within 
sinking  water  of  different  origins  in  the  context 
of  its  surface  supply,  we  also  estimate  the  possi¬ 
ble  excretory  release  of  DOM  by  various  groups  of 
plankton,  as  a  function  of  particulate  primary 
production  of  each  shelf,  slope,  and  basin  region 
[Platt  et  al.,  1991].  An  implicit  assumption  of 
this  analysis  is  that  oxidation  of  the  excreted  DOM 
withiii  sinking  waters  takes  place  at  depth,  such 


Walsh  et  al. :  Meridional  Fluxes  of  Dissolved  Organic  Matter 


15,627 


TABLE  1.  The  Estimated  Meridional  Fluxes  of  DOCpt  and  DOC„„ 
Across  36°25'N  in  the  North  Atlantic 


Depth 
Layer .  m 

Transport, 

Sv 

Water 

Type 

DOCpt(DOC, 

g 

Stock, 

D0Cpt(D0C„.„)  Flux, 

X  10®  kg  s"^ 

0-100 

15.2 

GS 

2.83 

(0.80) 

43.02 

(12.16) 

100-350 

-1.5 

shelf 

[20.00 

(10.80), 

[-30.00 

(-16.20), 

11.32 

(6.12), 

-16.98 

(-9.18), 

2.83 

(1.53)] 

-4.25 

(-2.30)] 

350-600 

2.2 

GS 

2.83 

(0.80) 

6.23 

(1.76) 

600-750 

-4.1 

NAC 

1.50 

(0.42) 

-6.15 

(-1.72) 

750-900 

2.4 

GS 

1.20 

(0.34) 

2.88 

(0.82) 

900-1200 

3.1 

AIW 

0.90* 

(0.25) 

2.79 

(0.78) 

1200-1400 

1.6 

GS 

1.20 

(0.34) 

1.92 

(0.54) 

1400-1800 

-2.1 

LCW/GS 

2.10 

(0.59) 

-4.41 

(-1.24) 

1800-2100 

-3.3 

LCW 

3.00 

(0.85) 

-9.90 

(-2.81) 

2100-2400 

-1.4 

LCW 

3.00 

(0.85) 

-4.20 

(-1.19) 

2400-2800 

-0.5 

LCW 

3.00 

(0.85) 

-1.50 

(-0.43) 

2800-3100 

-0.2 

LCW 

3.00 

(0.85) 

-0.60 

(-0.17) 

3100-3400 

1.1 

ACPW 

1.25 

(0.35) 

1.38 

(0.39) 

3400-3700 

-0.1 

NADW 

1.65 

(0.47) 

-0.17 

(-0.05) 

3700-4000 

-5.3 

NADW 

1.65 

(0.47) 

-8.75 

(-2.49) 

4000-4300 

-6.6 

NADW 

1.35 

(0.38) 

-8.91 

(-2.51) 

4300-4600 

-0.4 

AABW 

1.25 

(0.35) 

-0.50 

(-0.14) 

Total 

-1-0.1 

[-16.87 

(-12.50) , 

-3.85 

(-5.48), 

-H8.88 

(+1.40)] 

The  layer  transports  (Sv)  are  from  Rintoul  [1988]  and  the  DOCp^.uv 
concentrations  are  from  Bauer  [1991]  and  the  text.  The  water  types  are 
subducting  shelf  water  within  the  Gulf  Stream  (GS) ,  convecting  North  Atlantic 
Current  (NAC)  water,  cabelling  Labrador  Current  Water  (LCW) ,  and  sinking 
North  Atlantic  Deep  Water  (NADW) ,  Antarctic  Intermediate  Water  (AIW) , 
Antarctic  Circumpolar  Water  (ACPW) ,  and  Antarctic  Bottom  Water  (AABW) .  A 
negative  sign  denotes  southward  transport  of  water  and  DOCpt,  pv- 
*From  Bauer  [1991] . 


that  losses  of  photo -oxidation  and  bacterial  degra¬ 
dation  in  the  euphotic  zone  are  Ignored. 

3.  Results 

3.1.  Shelf  Waters 

We  first  use  aircraft  and  ship  bio -optical  data 
to  infer  possible  DOC  stocks  near  Cape  Hatteras 
(Figure  2),  where  perhaps  72%  of  the  particulate 
organic  carbon  (POC)  flux  of  the  Gulf  Stream  is 
derived  from  the  adjacent  South  and  Middle  Atlantic 
shelves  [Walsh,  1992].  In  contrast,  the  disparity 
in  half-lives  of  the  long-lived  DOC  and  the  short¬ 
lived  POC  pools  [Legendre  and  GosselLn,  1989] 
should  lead  to  much  longer  residence  for  DOC,  such 
that  some  of  the  DOM  off  Cape  Hatteras  may  have  a 
multiyear  origin,  both  from  waters  within  the  East 
Greenland  Current  [Chapman  and  Beardsley ,  1989]  and 
the  Mississippi  River  plume  [Atkinson  and  Wallace, 
1975].  This  situation  would,  of  course,  violate 
spatial  assumptions  of  linear  covariation  of  DOC, 
chlorophyll,  and  phaeopigments ,  needed  for  accurate 
satellite  algorithms  within  coastal  and  perhaps 
oceanic  waters  [Gordon  and  Morel,  1983]. 

Using  Coastal  Zone  Color  Scanner  (CZCS)  imagery, 
for  example,  the  mean  primary  production  of  North 
Atlantic  shelf  waters,  between  11°-70°N  (Table  2), 
is  estimated  to  be  732  g  C  yr'^  [Platt  et  al., 
1991]  ,  with  the  assumption  that  the  upwelling 


radiances,  sampled  by  the  satellite  radiometer, 
represent  mainly  chlorophyll.  This  is  not  always 
true ,  of  course ,  particularly  in  Case  2  coastal 
waters  [Morel  and  Prieur,  1977]  ,  where  humic  matter 
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Fig.  2.  A  prospectus  for  DOM  sequestration  within 
sinking  waters  of  the  North  Atlantic  during 
subduction  of  shelf  water  (area  1) ,  cabelling  of 
slope  water  (area  2),  convective  overturn  and 
lateral  induction  of  basin  water  within  NAC 
(area  3),  and  formation  of  NADW  (area  4). 
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TABLE  2.  Mean  Carbon  Fixation  (g  C  yr'^)  ,  Total  Photosynthesis 
(x  10^5  g  C  yr"i)  ,  Terrestrial  Input  (x  10^=  g  D0Cpt,uv  yr"^)  . 
Marine  Excretion  (x  10^^  g  DOC  yr'^)  ,  and  NAC  Contribution 
to  NADW  (x  10^^  g  DOCpt  „v  yt"^)  With  Respect  to  Sinking 
Waters  Within  Shelf,  Slope,  and  Basin  Areas 
(x  10®  km^)  of  the  North  Atlantic  Ocean 


Sinking 

Water 

Origin 

Latitude 

Mean  C 

Fixation 

Area 

Total 

Production 

Excretion/ 

Advection 

GS 

shelf 

11°-70°N 

732[366] 

2.5 

1.83[0.92] 

0.37[0.19] 

land 

11°-70°N 

- 

- 

0.35  Sv  X 
12(6.5)  g  m-3 

0.13(0.07) 

LCW 

slope 

38°-70'‘N 

329 

3.4 

1.12 

0.22 

basin 

38°-50°N 

319 

3.1 

0.99 

0.20 

NAG 

basin 

51°-70“N 

185 

5.4 

1.00 

0.20 

NADW 

shelf 

71°-80°N 

150 

0.7 

0.11 

0.02 

basin 

71°-80°N 

150 

0.5 

0.08 

0.02 

NAC 

71°-80°N 

- 

- 

2.00  Sv  X 

0.19(0.05) 

3(0.85)g  m‘® 

rder  et  al. 

,  1989]  . 

The  Gelbstoff, 

,  or  of 

CDOC  and 

marine  CDOC 

byproducts  o 

colored  dissolved  organic  carbon  (CDOC) ,  absorption 
coefficient  on  the  Alabama -Florida  shelves,  for 
example,  is  as  large  as  0.054  m'^  [Carder  et  al., 
1989]  ,  and  is  equivalent  in  absorption  to  a  chloro¬ 
phyll  concentration  of  1  mg  chi  m'®. 

More  traditional  ’■*C  estimates  of  carbon  fixa¬ 
tion  on  North  Atlantic  shelves,  between  11°-70°N, 
instead  suggest  a  mean  annual  production  of  159  g 
C  yr"^  [Walsh,  1988],  We  will  evaluate  these 
different  estimates  of  primary  production,  and 
associated  DOC  release,  after  estimation  of  the 
stocks  and  sinking  fluxes  of  DOM  of  different  shelf 
origins.  If  DOC  is  a  major  factor  in  the  contami¬ 
nation  of  color  signals  detected  by  passive  satel¬ 
lite  and  active  aircraft  remote  sensors  [Walsh, 

1988]  ,  what  are  the  implications  for  regional 
carbon  budgets? 

3.1.1.  DOC  stocks.  Away  from  coastal  regions, 
where  terrestrial  DOM  signals  are  the  strongest, 
the  abundance  of  marine  Gelbstoff,  and  therefore 
light  absorption,  generally  increases  with  the 
amount  of  annual  primary  production  [Carder  et  al.  , 

1989]  .  Within  the  Sargasso  Sea  during  December 
1984,  an  unproductive  time  of  year  [Menzel  and 
Ryther,  1960],  for  example,  the  ratio  of  the 
diffuse  attenuation  coefficients ,  k,j,  for  Gelbstoff 
or  CDOC,  and  kp,  for  the  particulate  matter  (chlo¬ 
rophyll  and  its  degradation  products) ,  at  441  nm 
ranged  from  0.3  to  0.7,  depending  on  the  optical 
model  [Topllss  et  al.,  1989].  This  range  is 
consistent  with  a  mean  value  of  0.6  at  440  nm  for 
the  respective  absorption  coefficients,  a,j/ap,  in 
the  Gulf  of  Mexico  [Carder  et  al . ,  1989],  since 
a  ~  k-cos  6,  where  6  is  the  subsurface  solar  zenith 
angle . 

At  December  shelf  and  slope  stations  in  the  MAB, 
k,j/kp  was  instead  a  mean  of  2.1  [Topliss  et  al., 
1989],  with  ratios  ranging  as  high  as  3.6,  the 
value  of  a,j/ap  observed  after  bloom  senescence  in 
coastal  waters  of  the  Gulf  of  Mexico  [Carder  et 
al. ,  1989] .  High  values  of  this  ratio  imply  excess 
amounts  of  colored  dissolved  organic  matter  (CDOM) 
on  the  continental  margins,  compared  to  steady 
state  conditions.  This  situation  results  from  an 
unknown  time  integral  of  both  terrestrial  influxes 


local  shelf  primary  production. 

The  winter  Topllss  et  al.  [1989]  data  can  be 
converted  to  estimates  of  DOCpt  by  first  estimat¬ 
ing  the  amount  of  CDOC.  Using  mass-specific 
absorption  coefficients  for  marine  humic  and  fulvic 
acids  [Carder  et  al. ,  1989],  an  average  cosine  of 
cos  $  »  0.82,  the  spectral  slope  of  0.013  nm"^  for 
the  winter  Topllss  et  al.  [1989]  data,  and  a  CDOC 
mixture  of  16%  humic  and  84%  fulvic  acids  [Carder 
et  al.,  1989],  we  compute  the  entries  of  Table  3 
from 

ad(A)  =  aa(450)  e'‘>-“®  («o-a)] 

where  3^(450)  =  0.031  CDOM,  and  CDOM  =  2  CDOC.  We 
obtain  means  of  2.2  and  0.17  g  CDOC  m‘®  due  to 
marine  humus  within  respective  shelf/slope  and 
basin  waters  (Table  3) . 

Is  such  a  10 -fold  winter  gradient  of  CDOC 
between  shelf  and  basin  waters  representative  of 
all  seasons,  and  does  it  reflect  similar  gradients 
of  total  DOC?  The  previous  wet  oxidation  tech¬ 
niques,  which  yield  full  recovery  of  humic  and 
fulvic  acids  (all  of  our  computed  CDOC  in  Table  3) , 
suggest  that  CDOC  constitutes  20-50%  of  the  total 
D0C„v  [Harvey  et  al.  ,  1983] ,  with  a  C/N  ratio  of  34 
[Heyers- Schulte  and  Hedges,  1986] .  However,  humic 
matter,  measured  by  the  platinum  method,  may 
account  for  only  6-14%  and  11-27%  of  the  total 
DOCpt,  since  Suzuki  and  Tanoue  [1991]  indicate  that 
DOCpt  may  be  3 . 54  x  DOCu.,  in  unproductive  areas  and 
1.85  X  DOCy.,  in  rich  regions. 

As  much  as  33%  of  the  DOCp^  may  instead  be 
reactive  polysaccharides  [Benner  et  al.,  1992]. 
Furthermore,  the  photo -chemical  production  of  low 
molecular-weight  organic  compounds  from  humic 
substances  in  natural  waters  [Kelber  et  al.,  1989, 
1990]  suggests  that  the  ratio  of  humic  (H)  and 
fulvic  (F)  acids  to  total  DOC  decreases  with  time 
in  surface  waters,  given  no  other  sources  or  sinks. 
This  process  is  perhaps  discernable  within  Harvey 
et  al.'s  [1983]  observations  for  depths  <  20  m  in 
the  Gulf  of  Mexico,  during  fertile  (spring)  and 
unproductive  (fall)  periods  of  the  year. 

Except  for  one  measurement  during  a  senescent 
bloom  off  Cape  San  Bias,  their  spring  data  [Harvey 
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TABLE  3 .  Mean  Pigment  Concentrations  (mg  m'^)  ,  Absorption  Coefficients  (m”^)  , 
and  Colored  Dissolved  Organic  Carbon  (g  m”^)  Within  MAB  Shelf/Slope  and 

Sargasso  Sea  Waters 


Chi  a 

Phaeo  a  ap(441) 

Ad  (441) 

CDOC 

DOCpt 

MAB  shelf/ 
slope 

1.84 

0.59 

0.072 

0.153 

2.20 

20.00 

Sargasso 

Sea 

0.35 

0.22 

0.025 

0.012 

0.17 

2.83 

These  values 

are  based 

on  both  Topliss  et  al. 

[1989]  measurements 

and  models , 

and  on  mass-specific  absorption  coefficients  for  a  mixture  of  15%  humic  and  85% 
fulvic  acids  [Carder  et  al.,  1989],  consistent  with  a  spectral  slope,  S;^,  of 
0.013  nm"^  for  the  absorption  by  colored  dissolved  organic  matter  (CDOM)  in 
December  1984.  Diffuse  attenuation  coefficients  were  converted  to  absorption 
coefficients  by  a  =  0.82  k,  while  CDOC/CDOM  =  0.5,  and  DOCpt  is  computed  from 
1.85  D0C^.„  in  coastal  waters  and  from  3.54  D0Cu„  in  oceanic  waters. 


et  al.,  1983]  can  be  expressed  by  the  regression 

DOC^^  =  0.30  +  0.78(H+F),  (2) 

whereas  the  fall  data  fall  along  a  line  expressed 
by 

D0Cu„  =  0.90  +  0.43(H+F)  (3) 

The  intercept  value  of  (3)  for  fall  conditions, 
i.e.,  the  portion  of  000^^,  which  does  not  covary 
with  (H+F) ,  is  threefold  that  of  (2)  for  spring 
waters.  Such  a  trend  may  reflect  photolysis  of 
some  of  the  light-absorbing  (H+F) ,  which  would  be 
reflected  in  a  smaller  CDOC/DOC  ratio. 

Accordingly,  we  compute  the  offshore  stock  of 
2.83  g  DOCpt  m"^  (Table  3)  in  the  Gulf  Stream  with 
a  CDOC/DOCpp  ratio  of  0.06,  and  the  shelf  stock  of 
20  g  DOCpt  m"^  with  a  ratio  of  0.11.  We  employ 
these  smaller  ratios  to  emphasize  both  photolysis 
of  CDOC  and  the  other  labile  fractions  of  DOM, 
captured  with  the  platinum  method,  since  our 
CDOC/DOCy.,,  ratio  is  0.20,  when  DOCpt  =1.85  D0C„.„  in 
coastal  waters,  and  0.21,  when  DOCpt  =  3.54  D0C„.,  in 
oceanic  waters  [Suzuki  and  Tanoue,  1991] .  Except 
for  the  shelf  estimates  of  Table  1,  all  other 
values  of  D0C^.„  at  different  depths  assume  oceanic 
DOCpt  =  3.54  D0C„^. 

Our  estimated  DOCp^  stock  for  offshore  surface 
waters  (Table  1)  approximates  the  previous  dry 
combustion  measurements  of  about  3.6  g  DOC  m"^  at 
the  edge  of  the  Nova  Scotian  shelf  and  2.4  g  DOC 
m"®  in  the  Sargasso  Sea  [Gordon  and  Sutcliffe , 
1973].  How  realistic  is  the  nearshore  estimate  of 
DOCpt  stocks  (Table  3)  ,  however?  Riverine  concen¬ 
trations  of  as  much  as  42.6  g  DOCpt  m‘®  (1.85  x 
prior  D0C„v  of  Wel-Bln  et  al.  [1983])  evidently 
suffer  rapid  dilution  on  the  adjacent  shelf  of  the 
East  China  Sea,  with  only  5.8  g  DOCp^  m“^  found  at 
mid-shelf  in  the  Yangtze  River  plume  and  3.1  g 
DOCpt  m“®  in  the  Kuroshio  [Suzuki  and  Tanoue,  1991] . 
We  thus  consider  a  range  of  three  DOCp^^uv  stocks 
for  shelf  waters  in  Table  1  and  compute  their 
offshore  fluxes  within  sinking  water,  for  compari¬ 
son  with  the  total  inputs  of  DOC  from  land  and  from 
release  by  the  marine  food  web  (Table  2) . 

3.1.2.  Sinking  water.  Containing  a  distinct 
diatom  assemblage  of  nearshore  species,  the  shelf 


water  at  Cape  Hatteras  is  subducted  to  depths 
>100  m  underneath  the  western  boundary  current 
[Lillibridge  et  al.,  1990].  Along  a  downstream 
distance  of  150  km,  for  example,  the  subsurface 
core  of  low  salinity  water  (and  associated  DOC 
stock)  descended  from  a  depth  of  -75  m  at  50  km 
seaward  of  the  Cape  Hatteras  shelf -break  to  -125  m 
during  an  October  1985  study  of  the  Gulf  Stream. 
Eastward  transfer  within  the  downstream  North 
Atlantic  and  Azore  Currents  over  months  to  years 
will  smear  such  land/shelf  DOC  signals,  of  course, 
like  that  of  radium  228  [Sarmiento  et  al.,  1990], 

We  estimate  the  possible  sinking  loss,  however, 
by  assuming  a  surface  width  of  low- salinity  waters 
of  shelf  origin  to  be  perhaps  15  km  along  the 
northwest  wall  of  the  Gulf  Stream  at  37°N.  Over  a 
depth  interval  of  at  least  100  m  at  a  speed  of  1  m 
s"'^,  the  transport  within  the  first  15  km  of  the 
nearshore  Gulf  Stream  would  be  1.5  Sv  of  water. 
This  transport  is  similar  to  that  of  1  Sv  estimated 
for  the  combined  water  fluxes  from  the  South  and 
Middle  Atlantic  shelves  at  the  100-m  isobath  off 
Cape  Hatteras  [Walsh,  1992]. 

If  subduction  were  to  occur  all  year,  continued 
storage  of  a  combination  of  terrestrial  and  shelf 
DOC,  like  radium  228  and  O2  [Sarmiento  et  al., 
1982,  1990] ,  may  occur  within  the  main  pycnocline 
of  the  subtropical  gyre  of  the  North  Atlantic. 
Here,  arguments  about  local  oxygen  utilization  and 
primary  production  rates  off  Bermuda  [Jenkins  and 
Goldman,  1985;  Platt  and  Harrison,  1986]  may  be 
confounded  by  these  lateral  supplies  of  organic 
matter,  albeit  DOC  rather  than  non- sinking  POC 
[Martin  et  al.  ,  1987]  .  Depending  upon  transit  time 
and  subsurface  respiration  losses,  the  southward 
dispersion  of  the  shelf  DOC  signals  may  be  a 
seasonal,  or  a  multiyear,  process. 

To  employ  a  suitable  averaging  scheme  across 
these  time  scales,  we  construct  a  zonal  mean  of  the 
possible  meridional  transport,  over  17  depth 
intervals,  of  DOCp^  across  36°25'N  latitude.  A 
southward  transport  in  Table  1  of  1.5  Sv  of  water 
over  a  depth  interval  of  approximately  100-350  m  (a 
density  interval  of  26.4  to  26.8  a  from  Roemmich 
and  Wunsch  [1985])  is  obtained  by  averaging  over 
the  6000-km  length  of  this  section  (Figure  1) . 
Such  a  southward  transport  of  water  and  oxygen  at 
this  depth  interval  [Rintoul,  1988]  is  assumed  by 
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us  to  also  represent  the  fate  of  shelf  and  land 
DOC,  subducted  at  the  surface  (Figure  2)  along 
-40°N  between  20°-70°W  [Huang,  1990] ,  with  possible 
equatorward  transfers  of  4.25  to  30  x  10^  kg  DOCp^ 
s'^  into  the  upper  pycnocline  of  the  subtropical 
gyre  (Table  1) . 

Although  the  vertical  structure  of  the  age 
of  DOCpt  in  the  water  column  off  Bermuda  has  not 
yet  been  sampled  with  the  same  resolution  as 
Table  1,  early  results  [Bauer,  1991]  confirm  our 
assumptions  of  a  recent  source  of  DOC  within  the 
upper  pycnocline.  At  a  depth  of  about  500  m,  a 
subsurface  minimum  in  the  radiocarbon  age  of  DOCp^ 
was  observed  (about  -27511o),  indicating  a  relatively 
young  origin.  In  contrast,  a  A^'‘C  value  of  -450?tt 
was  measured  for  a  DOC  stock  of  0.9  g  DOCp^ 
(Table  1)  within  the  oxygen  minimum  layer  near  a 
depth  of  850  m  [Bauer,  1991] . 

3.1.3.  Excretory  supply.  If  we  first  assume 
the  combined  land  and  shelf  DOC  stock  to  be  20  g 
DOCpt  m'^,  the  southward  DOC  flux  along  the  north¬ 
west  wall  of  the  Gulf  Stream  might  be  as  high  as  30 
X  10^  kg  DOCpt  (Table  1),  or  0.9  x  10“  g  DOCpt 
yr'^  on  an  annual  basis.  This  is  presumably  an 
upper  bound,  since  it  is  almost  twice  the  larger 
estimate  of  land  and  shelf  supply  of  DOC  to  the 
North  Atlantic  shelves  (Table  2)  .  A  smaller  DOC^,^ 
shelf  stock  of  10.8  g  m'^  and  the  same  water 
transport  of  1.5  Sv  yield  instead  an  export  of  16.2 
X  10='  kg  DOC,,^  s'l  (Table  1),  or  0.5  x  10“  g  D0C„^ 
yr’^;  this  is  still  larger  than  the  estimated 
inputs  of  terrestrial  and  marine  DOCu^. 

A  very  large  total  primary  production  of  1.83  x 
10=-^  g  C  yr'=-  over  a  shelf  area  of  2.5  x  10®  km^ 
[Platt  et  al.,  1991],  for  example,  and  a  phyto¬ 
plankton  excretion  rate  of  20%  of  this  carbon 
fixation  yield  a  marine  excretory  input  of  0.37  x 
10=-®  g  DOC  yr'^-  (Table  2) .  A  total  freshwater 
discharge  of  1.1  x  lO*-®  m®  yr"'-  to  the  North  Atlan¬ 
tic  shelves  between  11°-70°N  [Walsh,  1988]  and  a 
mean  terrestrial  DOC  concentration  of  12  g  DOCp^  m'® 
(1.85  X  6.5  g  DOCy^  m'®  of  Spitzy  et  al.  [1991]) 
suggest  a  possible  land  input  of  either  0.7  x  10“ 
g  DOCuv  yt'"-,  or  1.3  X  lO"-®  g  DOCpt  yt"^.  If  these 
estimates  of  photosynthesis,  excretion,  and  terres¬ 
trial  supplies  of  DOC  are  all  correct,  the  combined 
sources  sum  to  0.44  x  lO-'®  g  DOC^,.,,  yr“=^  and  0.50  x 
10“  g  DOCpt  (Table  2)  . 

Mean  excretion  rates  of  shelf  phytoplankton 
range  from  5%  (n=21)  of  primary  production  in  the 
South  Atlantic  [Thomas,  1971],  17%  (n=838)  in  the 
Middle  Atlantic  [Thomas  et  al.  ,  1979]  ,  to  38%  (n=4) 
in  the  North  Atlantic  [Choi,  1972]  Bights,  suggest¬ 
ing  an  average  of  20% .  Such  a  shelf  DOC  release 
rate  is  similar  to  respective  means  of  22-23%  for 
phytoplankton  excretion  of  DOC  in  slope  (n=4)  and 
basin  (n-5)  regions  of  the  western  North  Atlantic 
[Thomas,  1971;  Choi,  1972].  When  these  data  are 
pooled  with  the  results  of  14  other  regional 
studies  [Sharp,  1977] ,  a  grand  mean  excretion  rate 
of  23%  is  obtained. 

Sharp  [1977]  suggested  that  the  high  excretion 
rates  of  phytoplankton  might  be  experimental  arti¬ 
facts,  but  analysis  of  more  recent  data  suggests  a 
mean  rate  of  13%,  which  meets  less  than  half  of  the 
substrate  demands  of  coexisting  bacterial  popula¬ 
tions  [Baines  and  Pace,  1991] .  Thus  far,  we  have 
Ignored  releases  of  DOC  by  zooplankton  [Jumars  et 
al. ,  1989]  and  bacterioplankton  [Brophy  and 

Carlson,  1989]  ,  thereby  underestimating  the  biotic 
supply  of  DOC.  An  excretory  release  rate  of  20% 
approximates  the  numerous  observations  of  phyto¬ 


plankton  excretion  in  the  MAB  [Thomas  et  al.  ,  1979] 
and  allows  for  some  heterotrophic  input  of  DOC  as 
well  [Baines  and  Pace,  1991]. 

We  believe  that  the  satellite  estimate  of  a 
shelf  primary  production  of  732  g  C  m"^  yr'"-  [Platt 
et  al. ,  1991]  may  be  too  large,  however,  because  of 
coastal  stocks  of  DOC,  sensed  as  chlorophyll 
biomass  by  the  CZCS.  If  the  remotely  sampled  DOC 
concentrations  constitute  half  of  the  color  signal 
of  coastal  waters  [Hochman,  1992]  ,  and  the  specific 
photosynthetic  rate  (yr""')  remains  the  same,  the 
actual  primary  production  might  be  366  g  C  m"^  yr’=- 
(Table  2).  A  smaller  excretory  input  of  0.19  x 
10“  g  DOC  yr""-  would  follow  (Table  2)  ,  supporting 
less  of  a  shelf  export  of  DOC  to  be  subducted 
beneath  the  Gulf  Stream.  We  thus  consider  the  case 
of  an  intermediate  seaward  DOC  flux,  estimated  by 
aircraft  sensors ,  and  then  a  third  scenario  of  no 
difference  between  shelf  and  basin  stocks  of  DOC. 

Active  remote  sensors  on  aircraft  [Hoge  and 
Swift,  1981]  sample  at  10®  faster  rates  than 
shipboard  surveys  [Walsh  et  al.,  1988],  such  that 
the  laser- induced  measurement  of  a  fourfold  range 
in  Gelbstoff  concentrations,  between  coastal  and 
offshore  waters  in  December  [Topliss  et  al.  ,  1989]  , 
may  be  more  accurate  than  the  spatial  assessment  of 
DOC  abundance  from  a  few  ship  stations.  According¬ 
ly,  we  estimate  that  the  shelf  stock  of  DOM  might 
instead  be  an  intermediate  11.32  g  DOCpt  m"®, 
fourfold  the  basin  concentration  (Table  1).  In 
this  case,  the  offshore  fluxes  of  DOM  would  be  17 
X  10®  kg  DOCpt  and  9.2  x  10®  kg  DOC^y  s'® 
(Table  1).  Over  a  year,  these  possible  sinking 
fluxes  of  0.3  X  10®®  g  DOCyy  yr"®  -0.5  x  10®®  g  DOCpt 
yr'®  represent  all  of  the  combined  terrestrial  and 
shelf  supplies  of  0.26-0.32  x  10®®  g  DOC  yr"® 
(Table  2). 

Finally,  the  no -DOC -gradient  scenario  yields 
offshore  fluxes  of  2.3  x  10®  kg  DOC^y  s"®  and  4.2  x 
10®  kg  DOCpt  (Table  1),  or  0.07-0.13  x  10®®  g  DOC 
yr"®,  i.e.,  only  27-41%  of  the  coastal  supplies  of 
DOC  at  a  reduced  annual  rate  of  photosynthesis  of 
366  g  C  m"^  yr"®  (Table  2) .  The  offshore  decline  of 
primary  production,  from  732  [366]  to  329  g  C  m"^ 
yr"®  (Table  2) ,  may  reflect  less  contamination  by 
DOC  of  the  satellite  estimates  of  chlorophyll 
biomass  in  slope  regions.  Within  basin  regions 
between  38°-50°N  and  51°-70°N,  the  respective  rate 
of  photosynthesis  is  even  smaller,  from  319  to  185 
g  C  m"®  yr"®  (Table  2)  ;  we  thus  assume  that  the 
estimates  of  primary  production  in  slope  and  basin 
waters  [Platt  et  al.,  1991]  are  correct  for  the 
purposes  of  our  carbon  budgets . 

3.2.  Slope  Waters 

3.2.1.  DOC  stocks.  We  next  consider  how  much 
surface  DOC  stocks  might  be  cabelled  downward  as  a 
mixture  of  the  Labrador  and  North  Atlantic  Currents 
to  depths  of  1400-2400  m  on  the  Newfoundland  slope 
[iCrauss  et  al.,  1990].  Of  the  90-100  Sv  transport 
of  the  Gulf  Stream  at  70°W  [Johns  et  al.  ,  1989]  and 
55°W  [Richardson,  1985],  part  recirculates  to  the 
west,  part  moves  southeast  as  the  Azores  Current 
(Figure  2) ,  and  about  33  Sv  proceeds  northeastward 
as  the  NAC  [Krauss  et  al.,  1990]. 

With  movement  towards  the  northeast  at  0. 5-1.0 
m  s"®,  between  the  subarctic  and  subtropical  gyres 
of  the  North  Atlantic  [Smith  et  al .  ,  1990],  this 
downstream  western  boundary  current,  east  of  the 
Grand  Banks,  should  still  have  a  Gulf  Stream  DOC 
content  of  2.83  g  DOCpt  m"®  (Table  1),  assuming  a 
half-life  for  the  labile  DOC  of  about  25-50  years 
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[Kelber  et  al.,  1990;  Toggweiler,  1990;  Bacastow 
and  MaLer-Reimer ,  1991]  .  We  then  need  an  estimate 
of  the  DOC  content  of  the  coastal  end  member  of 
this  mixing  pair. 

The  actual  source  regions  of  shelf  DOC  entrained 
at  Cape  Hatteras  may  extend  from  north  of  Greenland 
{Chapman  and  Beardsley ,  1989]  to  at  least  the 

Mississippi  River  [Atkinson  and  Wallace,  1975]  as 
a  result  of  longshore  coherence  of  flow  and  the 
long  half-life  of  some  DOC  components.  In  contrast 
to  the  Mississippi  and  Hudson  Rivers,  where  urban 
DOC  fluxes  constitute  as  much  as  90%  of  the  daily 
carbon  loading  of  8-10  g  D0Cu.„  m"^  [Mueller  et  al., 
1976;  Segar  and  Berberian,  1976;  Thomas  et  al.  , 
1979;  Leenheer,  1982],  the  DOC  content  of  the 
relatively  pristine  St.  Lawrence  River  is  now  less 
than  half,  about  2-4  g  DOC^,,  m‘^  [Pocklington  and 
Tan,  1983] .  We  might  thus  expect  the  additional 
DOC  loading  of  the  western  North  Atlantic  from 
subarctic  shelf  waters  (Figure  2),  partially 
derived  from  pristine  Arctic  Rivers ,  to  be  less 
than  at  mid- latitudes . 

Accordingly,  we  assume  that  the  2  Sv  of  shelf 
water,  draining  the  Canadian  Archipelago  via  Davis 
Strait  [Chapman  and  Beardsley,  1989;  Walsh  et  al. , 
1989]  contains  6  g  DOCpt  m'^  (1.85  x  3.24  g  DOC^^  m'® 
[Pocklington  and  Tan,  1983]).  A  mixture  of  2  Sv  of 
this  shelf  water,  with  6  g  DOCpt  m'^,  and  33  Sv  of 
prior  Gulf  Stream  water,  with  2.83  g  DOCpt 
yields  a  recirculating  NAC  transport  (Figure  2)  of 
35  Sv,  containing  3  g  DOCpt  m'^  (Table  1) .  Part  of 
this  augmented  NAC  sinks  in  the  Labrador  Sea  as 
cabelled  slope  water,  and  part  sinks  in  the  eastern 
Atlantic  as  convected  basin  water  (Figure  2) . 

3.2.2.  Sinking  water.  At  least  4  Sv  of  water 

cabells  along  the  Newfoundland  slope  to  join  the 
DWBC  at  depths  >2400  m  [Krauss  et  al. ,  1990] .  This 
estimate  is  similar  to  prior  transport  calculations 
of  sinking  surface  water  in  this  region  [Sverdrup 
et  al.,  1942;  Wright,  1972].  Over  a  year,  continu¬ 
ous  sinking  of  this  amount  of  water  would  yield  an 
annual  slope  DOC  sequestration  of  0.1  x  10“^^  g  DOCu^ 
yr"^  to  0.4  x  10^^  g  DOCp^  yr"^,  i.e.,  half,  or 

double,  the  estimated  excretory  supply  of  DOC 
within  just  slope  regions  over  38°-70°N  (Table  2). 

3.2.3.  Excretory  supply.  Over  the  northern 
slope  area  of  3.4  x  10®  km^,  a  total  photosynthesis 
of  1.12  X  10^®  g  C  yr'^  [Platt  et  al.,  1991]  and  a 
20%  excretion  rate  yield  an  excretory  release  of 
0.22  X  10^®  g  C  yr"^  within  the  nearshore  mixing 
member  (Table  2).  Over  a  similar  basin  area  of  3.1 
X  10®  km^,  between  38°-50°N,  more  primary  production 
of  0.99  X  10^®  g  C  yr"^  yields  another  excretory 
input  of  0.20  X  10^®  g  C  yr'^  (Table  2).  Their  sura 
balances  the  advective  loss  of  DOCpt  within  4  Sv  of 
sinking  Labrador  Current  Water  (LCW) . 

If  the  DOCpt  measurements  are  artifacts ,  how¬ 
ever,  as  a  result  of  organic  contaminants  of  the 
catalysts  or  poor  blanks,  for  example,  less  biotic 
supply  of  DOC  is  required  to  balance  sinking  losses 
of  D0C„.„  within  LCW.  The  sequestration  of  0.1  x 
10^®  g  DOCpv  yr"^  within  4  Sv  of  LCW  only  requires 
the  supply  of  dissolved  photos3mthate  over  a 
smaller  area  of  -1.6  x  10®  km^,  i.e.,  along  the 
Labrador  slope  between  50°-60°N.  On  the  other 
hand,  if  all  of  the  southward  transport  of  6.8  Sv 
of  LCW  across  36°25'N  (Table  1)  represents  just  one 
year  of  ventilation  of  surface  waters ,  and  the 
DOCpt  stocks  are  not  artifacts,  the  associated 
equatorward  flux  of  DOC  would  be  0.6  x  10^®  g  DOCpt 
yr"'^,  i.e.,  150%  of  the  sum  of  slope  and  basin 

supplies  of  DOC  (Table  2) . 


3.2.4.  Southern  Ocean  supply.  In  terms  of  a 
meridional  flux  of  DOC  across  36°25'N,  we  thus 
assume  that  5.4  Sv  of  LCW,  over  a  depth  range  of 
1800-3100  m,  represents  a  recent  influx  from  the 
Labrador  Sea.  The  relict  DOC  stocks  at  1400-1800 
m  are  instead  assumed  to  be  a  mixture  of  LCW  and 
Gulf  Stream  waters ,  in  which  the  lower  values  of 
DOC  represent  respiration  losses  of  multiyear 
inputs  from  Antarctic  source  regions ;  specifically 
within  Antarctic  Intermediate  Water  (AIW  of 
Table  1) .  Note  that  the  DOC  pools  of  Antarctic 
Circumpolar  Water  (ACPW)  and  of  Antarctic  Bottom 
Water  (AABW)  are  estimated  to  be  1.25  g  DOCp^  m"® 
(0.35  g  D0Cp.„  m"®) ,  somewhat  higher  than  those  of 
the  oxygen  minimum  layer  of  AIW  [Bauer,  1991] . 

These  additional  stocks  of  Southern  Ocean  origin 
are  assumed  to  have  initial  surface  values  of  1.5  g 
DOCpt  m"®  (0.42  g  DOC^^  m"®)  ,  half  that  estimated  for 
Arctic  waters,  to  be  consistent  with  observations 
of  k,j  in  the  Greenland,  Barents,  and  Bellingshausen 
Seas  [Mitchell,  1992].  Based  on  302  bio-optical 
measurements  in  these  Seas,  the  diffuse  attenuation 
coefficient  due  to  dissolved  substances,  k,j,  at 
441  nm  is  0.036  m‘^  in  the  Antarctic,  compared  to 
0.063  m"®^  in  the  Arctic.  This  situation  may  result 
from  greater  terrestrial  loading  [Mitchell,  1992], 
and/or  larger  marine  photosynthetic  [Walsh,  1989] 
release,  of  DOC  in  the  Arctic. 

3.3.  Basin  Waters 

To  explore  the  basin  (>2000  m)  sinking  losses  of 
DOC,  we  estimate  how  much  might  be  sequestered  in 
the  North  Atlantic  by  two  convective  processes : 
(1)  overturn  and  lateral  induction  of  NAC  waters 
into  the  subtropical  gyre,  south  of  50°N  [Huang, 
1990],  and  (2)  overflow  of  NADW  [Swift,  1984] 
across  the  600-m  sill  of  Denmark  Strait  at  65°N. 

3.3.1.  MAC  transfers .  Before  transformations 
of  the  NAC  into  the  Irminger,  Norwegian,  and 
Portugal  Currents  (Figure  2) ,  seasonal  convection 
occurs  en  route  in  the  northern  North  Atlantic, 
with  deepening  of  the  surface  mixed  layer  [Levitus, 
1982]  .  North  of  bO^N,  for  example,  where  the  depth 
of  the  surface  mixed  layer  in  February  can  exceed 
700  m  [McClain  et  al . ,  1990],  convective  overturn 
of  part  of  the  NAC  may  result  in  lateral  induction 
of  an  equatorward  transport  of  at  least  4.1  Sv, 
which  could  reach  a  depth  interval  of  600-750  m 
across  36°25'N  (Table  1). 

With  a  NAC  stock  of  3  g  DOCpt  m‘®,  this  transfer 
would  amount  to  a  southward  flux  of  12.3  x  10®  kg 
DOCpt  t>r  0.4  X  10®®  g  DOCpt  over  12  months. 
If  this  shallow  portion  of  the  basin  export  of  DOC 
from  the  northern  North  Atlantic  does  occur  year 
round,  it  represents  twice  the  estimated  annual 
release  of  DOC  within  surface  waters  between  51°- 
70°N  (Table  2) . 

A  small  carbon  fixation  of  185  g  C  m"®  yr'®  over 
5.4  X  10®  km®  and  a  20%  excretion  rate  yield  a 
release  of  0.2  x  10®®  g  DOC  yr'®  (Table  2).  Such 
low  rates  of  photosynthesis  [Thomas,  1971]  may  lead 
to  larger  percentages  of  DOC  excretion  by  oceanic 
microflagellates  [Wolter,  1982],  perhaps  even 
during  spring  blooms  of  this  region  [Kirchman  et 
al.,  1991].  A  larger  excretion  rate  of  40%,  for 
example,  would  allow  a  biotic  influx  of  0.4  x  10®® 
g  DOC  yr‘®. 

An  alternative  hypothesis  is  a  seasonal  decline 
of  photosynthesis  and  of  DOC  release  during  light- 
limited  periods  of  the  year,  when  initiation  of  the 
southward  transfer  of  overturned  NAC  occurs.  The 
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mean  depth  of  the  surface  mixed  layer  at  59°30'N, 
20°30'W  deepens,  for  example,  from  45.1  m  during 
May-October  to  432.2  m  during  November-April  (based 
on  the  AT/Az  criterion  of  0.5°C  [Levitus,  1982]). 
Finally,  a  third  scenario  would  be  to  discard  the 
DOCpt  estimate  of  NAG  stocks  in  favor  of  an  annual 
southward  NAG  flux  of  0.11  x  lO’-^  g  DOG^.,  yr‘^. 

Since  Duursms  [1965]  suggested  that  there  might 
be  a  seasonal  influx  of  surface  Arctic  DOG  stocks 
to  deep  waters  via  sinking  NADW  as  well,  we 
implement  the  second  hypothesis  in  Table  1.  We 
represent  this  process  by  a  50%  reduction  of  the 
NAG  stock  to  1.5  g  DOGp^  m"^  (Table  1)  to  reflect 
lower  winter  values ,  yielding  an  annual  export 
across  36°25'N  of  0.2  x  10^^  g  DOGpt  yr"^. 
Validation  of  this  hypothesis  of  seasonal  DOG 
release  is  seen  in  the  higher  summer  values  of 
DOGpt  within  surface  waters  between  47°-59°N 
(Figure  3) .  A  similar  logic  is  invoked  for  DOG 
export  within  sinking  NADW. 

3.3.2.  NADW  transfers .  A  continuous  tongue  of 
relatively  high  DOG,  about  0.7  g  DOG  m"^  with  early 
techniques  [Duursma,  1965],  could  be  traced  from 
surface  to  bottom  waters  south  of  Iceland,  within 
NADW  during  April  1958,  but  not  in  the  following 
September.  More  recent  data  suggest  intermediate 
values  of  0.8 -1.2  g  DOG^.,,  m"^  at  the  surface  of  the 
Lincoln  Sea  [Gordon  and  Cranford,  1985],  compared 
to  prior  observations  of  2-3  g  DOG  m"^  in  the 
Ganadian  Basin  [Melnikov  and  Pavlov,  1978]. 
Assuming  a  3.54  underestimate  of  D0Gp.„  in  these 
unproductive  regions  (Table  2)  yields  a  possible 
mean  stock  of  3.5  g  DOGpt  m'^  within  Arctic 
precursors  of  NADW  in  the  Greenland  and  Iceland 
Seas . 

If  we  further  assume  that  all  5  Sv  of  the 
transport  of  the  East  Greenland  Gurrent  (3  Sv  of 
polar  origin,  2  Sv  of  Atlantic  origin  [Aagaard  and 
Reed,  1987])  are  converted  to  NADW  above  Denmark 

DOC  (g  m'®) 


Fig.  3.  The  vertical  structure  of  DOGpt,  as 
measured  by  Y.  Suzuki,  during  May- June  1989  in 
temperate  and  subarctic  regions  of  the  eastern 
North  Atlantic  (see  Figure  1  for  the  locations  of 
these  regions) . 


Strait,  we  obtain  the  same  estimate  for  NADW 
formation  as  Worthington' s  [1976]  calculation:  4  Sv 
through  Denmark  Strait  and  1  Sv  across  the  Iceland- 
Scotland  ridge .  Such  a  mixture  of  the  NAG  and  the 
East  Greenland  Gurrent  end  members  yields  a 
possible  DOG  stock  for  NADW  of  3.3  g  DOGp^  m'^  (0.9 
g  DOG^.^,  m"®)  .  Over  12  months,  the  NADW  export  of 
DOG  from  the  Iceland  and  Greenland  Seas  might  then 
be  either  0.15  x  10^^  g  D0Gu.„  yr"^,  or  0.52  x  10^^  g 
DOGpt  i.e.,  ignoring  seasonality  in  the 
release  and  sequestration  of  DOG. 

Once  more,  this  estimate  is  an  upper  bound, 
because  a  mean  carbon  fixation  of  150  g  G  m'^  yr'^ 
[Ellertsen  et  al.,  1989],  over  both  the  0.7  x  10® 
km^  shelf  of  the  Barents  Sea  [Walsh,  1988]  and  the 
0.5  x  10®  km^  basin  of  the  Greenland  Sea  [Carmack 
and  Aagaard,  1973],  only  sums  to  0.19  x  10^®  g  G 
yr“^  (Table  2)  .  An  excretion  rate  of  20%  then 
provides  a  combined  shelf/basin  influx  of  at  most 
0.04  X  10^®  g  DOG  yr'^  within  local  waters.  The 
supply  of  2  Sv  of  NAG  from  the  south,  with  a  stock 
of  3  g  DOGpt  contribute  0.19  x  10^®  g  DOGpt 
yr'i  to  NADW  (Table  2). 

There  is  presumably  no  net  import  of  DOM  from 
the  Arctic  Ocean  via  Fram  Strait,  however,  since 
the  same  63  mg  DON  m‘®  is  found  [Kattner  and 
Becker,  1992]  under  the  ice  in  the  East  Greenland 
Gurrent  (Transpolar  Drift  origin)  and  in  open 
waters  of  the  West  Spitsbergen  Gurrent  (NAG 
origin) ,  whose  transport  is  also  about  5  Sv  [Walsh 
et  al.,  1989].  The  biogenic  and  advective  sources 
of  DOG  within  sinking  NADW  may  thus  sum  to  0.09  x 
lOi®  g  DOGu^  yr'i,  or  0.23  x  lO^®  g  DOGpt  y^'S  i.e., 
half  that  of  the  above  DOG  exports . 

By  summing  a  constant  flux  over  the  year,  we  may 
have  again  overestimated  the  annual  sequestration 
of  DOG  within  sinking  waters  of  the  northern  North 
Atlantic.  Accordingly,  we  assume  a  seasonal  supply 
of  DOG  within  NADW  (see  Figure  3  for  the  spatial 
analogue),  like  that  of  NAG,  with  a  winter  estimate 
of  1.65  g  DOGpt  (Table  1),  when  NADW  is  presum¬ 
ably  formed.  As  in  the  case  of  LOW,  we  also  assume 
that  only  5.4  Sv  of  NADW  transfer,  over  depths  of 
3400-4000  m  along  36°25'N,  represent  the  annual 
carbon  loading  of  0.28  x  10^®  g  DOGpt  y^'^ 
(Table  1). 

3.3.3.  Southern  Ocean  transfers .  The  addition¬ 
al  6.6  Sv  of  entrained  NADW,  over  depths  of  4000- 
4300  ra,  instead  represent  a  multiyear  accumulation 
of  both  Antarctic  and  Arctic  sources ,  with  an 
intermediate  stock  of  1.35  g  DOGpt  ^  (Table  1). 
Before  we  turn  to  the  implications  of  these  meridi¬ 
onal  fluxes  over  the  whole  water  column  at  36°25'N, 
we  first  address  how  realistic  our  estimates  of  the 
DOG  stocks  might  be.  An  analysis  is  made  of  recent 
data  taken  during  the  1989  and  1991  spring  blooms 
of  the  North  Atlantic . 


4.  Gorroboration 
4.1.  Vertical  Structure 

During  May  18-19  and  June  5-6,  1989,  Y.  Suzuki 
measured  DOGpt  concentrations  at  several  stations 
[Slagle  and  Helmerdinger ,  1991]  near  59°30'N, 
20°25'W  and  46°30'N,  17°45'W  (locations  A  and  B  of 
Figure  1) ,  To  provide  sufficient  depth  resolution, 
we  plotted  his  data  from  two  adjacent  stations  at 
each  location  in  Figure  3,  which  show  possible 
signatures  of  cabelled  slope  water  and  overturned 
basin  water. 
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The  measured  DOC  stocks  within  the  upper  25  m  of 
the  water  column  in  region  A  (59°30'N)  yield  a  mean 
(n=ll)  of  1.9  g  DOCpt  m'® ,  compared  to  a  mean  (n=10) 
of  2.4  g  DOCpt  m”^  in  region  B  (46°30’N),  where  a 
maximum  of  3.0  g  DOCp^  m'®  was  found.  These  values 
at  the  lower  latitude  are  equivalent  to  our  esti¬ 
mates  of  2. 8 -3.0  g  DOCpt  Gulf  Stream  waters 
and  LCW  along  36°25'N  (Table  1).  Note  that  the 
subsurface  values  of  DOCpt  the  higher  latitude 
are  almost  half  those  of  surface  waters ,  reflecting 
increased  excretion  during  the  spring  bloom 
[Kirchman  et  al. ,  1991]. 

Preservation  of  the  vertical  structure  of  DOC 
partially  depends  upon  latitude,  i.e.,  the  extent 
to  which  local  vertical  mixing  smears  the  DOC 
tracers  introduced  at  the  surface .  Region  A  at 
59°30'N  has  a  sufficiently  deep  winter  mixed  layer 
[levitus,  1982]  to  allow  uniform  stocks  of  1.2  g 
DOCpt  over  depths  of  200  to  500  m  (Figure  3) . 
Fall  overturn  of  this  part  of  the  summer  water 
column  at  59°30'N  would  indeed  lead  to  a  depth- 
averaged  stock  of  1.26  g  DOCpt  a  mean 
density  of  27.37  cr. 

If  the  overturned  basin  waters  at  59°30'N  were 
to  flow  south  along  the  27.3-27.5  a  isopycnals 
after  Ekman  pumping  [Huang,  1990]  ,  these  DOC  stocks 
would  then  occur  at  deeper  depths  of  700-900  m  near 
both  46°30'N  and  36°25'N,  where  we  assumed  a  stock 
of  1.5  g  DOCpt  (Table  1).  The  vertical 
structure  of  DOC  at  46°30'N  appears  to  provide 
similar  validation  of  our  hypothesis  of  LCW  seques¬ 
tration  of  excretory  products . 

The  experimental  precision  of  the  DOCpt  measure¬ 
ments  is  ±  1-5%,  depending  upon  sea  state  [Slagle 
and  Helmerdinger,  1991] ,  such  that  the  10%  subsur¬ 
face  increment  of  DOC  at  depths  of  1500-2000  m  in 
region  B,  but  not  in  region  A  (Figure  3),  may 
reflect  cabelling  LCW.  Since  LCW  and  NADW  both 
travel  south  mainly  in  the  DWBC  of  the  western 
Atlantic  Ocean,  as  Indicated  by  tritium  dispersal 
and  radiocarbon  ages  [Broecker  et  al.,  1991],  we 
would  expect  their  signals  to  be  weakest  in  the 
northeastern  Atlantic  along  20°W.  Indeed  at 
greater  depths  of  2700  m  in  regions  A  and  B,  the 
respective  stocks  of  0.8  and  1.4  g  DOCp^  m”®  are 
equivalent  to  our  assumed  DOCpp  content  of  ACPW  and 
AABW  (Table  1) . 


4.2.  Horizontal  Structure 

A  transect  of  nine  stations,  from  site  A,  south 
of  Iceland,  to  Nova  Scotia  during  May- June  1991 
(Figure  1) ,  provides  an  assessment  of  the  various 
values  of  shelf  DOC,  assumed  to  be  subducted  under 
the  Gulf  Stream  at  depths  of  100-350  m  (Table  1). 
The  surface  values  of  CDOM  absorption  at  440  nm 
show  a  linear  trend  with  salinity,  inshore  of  the 
Labrador  Current  (Figure  4) .  They  are  threefold 
higher  in  summer  waters  of  the  Nova  Scotian  shelf 
(Table  4),  compared  to  either  site  A,  or  the 
hibernal  Sargasso  Sea  (Table  3) .  The  shelf  obser¬ 
vations  of  0.052-0.078  m"'-  are  similar  to  that  of 
0.063  m"^  in  the  Arctic  [Mitchell,  1992],  but  half 
the  absorption  values  observed  [Topliss  et  al.  , 
1989]  in  MAB  shelf  waters  (Table  3). 

The  spectral  absorption  coefficients  of  these 
CDOM  data  were  converted  [Carder  et  al.,  1989]  to 
estimates  of  humic  and  fulvic  acids  (Table  4) . 
Equation  (2)  for  a  spring  bloom  situation  of  less 
photolysis  was  then  used  to  convert  (H+F)  to  DOCu^, 
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Fig.  4.  The  relation  of  Gelbstoff  absorption  at 
440  nm  to  surface  salinity  between  Iceland  and 
Nova  Scotia  (see  Figure  1  for  locations  of 
stations  1-9)  . 


TABLE  4.  Absorption  Coefficients  (a^  440)  due  to  DOM  Within  the  Upper  Meter  at 
Stations  of  the  1991  Transect  From  Iceland  to  Nova  Scotia 


Station 

440,  m'^ 

S;^,  nm"'- 

H+F,  g  m'3 

D0Cu„,  g  m‘® 

DOCpt,  S 

1 

0.0256 

-0.0150 

1.24 

1.26 

2.27 

2 

0 . 0344 

-0.0133 

1.12 

1.17 

2.11 

3 

0.0333 

-0.0146 

1.47 

1.44 

2.59 

4 

0.0325 

-0.0150 

2.05 

1.89 

3.40 

5 

0.0327 

-0.0156 

2.23 

2.03 

3.65 

6 

0.0603 

-0.0141 

2.90 

2.55 

4.59 

7 

0.0519 

-0.0159 

4.05 

3.44 

6.19 

8 

0.0619 

-0.0160 

4.84 

4.05 

7.29 

9 

0.0775 

-0.0165 

4.90 

4.10 

7.37 

The  samples  were  filtered  with  preflushed  0.2  pm  pore  diameter,  Nucleopore 
pads,  and  measured  in  10  cm  cells  of  a  Cary  2200  spectrophotometer.  The  spectral 
slope  of  the  semi-log  regression  between  310  and  490  nm  is  denoted  as  S;^,  while 
the  humic  (H)  and  fulvic  (F)  acid  concentrations  are  estimated  from  a^  450,  S;^, 
and  the  equations  provided  by  Carder  et  al.  [1989]  .  Our  estimates  of  D0C„y  are 
then  derived  from  the  H+F  values,  as  discussed  in  the  text. 
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yielding  values  of  1.16  to  1.51  g  DOC^^  m'®  for 
stations  1-3  at  site  A  (Table  4) .  With  our  eutro- 
phic  DOCpt/DOCuv  ratio  of  1.85,  the  estimated 
surface  DOC  stocks  at  site  A  in  1991  become  2.15- 
2.79  g  DOCpt  (Table  4),  compared  to  observations 
of  1.50-1.97  g  DOCpt  m"®  during  1989  (Figure  3). 

While  these  offshore  observations  are  from 
different  years,  they  are  from  the  same  site  in  the 
same  season,  providing  a  check  on  our  hypotheses  of 
various  gradients  of  DOC  between  shelf  and  basin 
regions  of  the  North  Atlantic.  Our  blo-optlcal 
measurements  Indicate  a  rapid  increase  of  DOM  on 
the  Grand  Banks,  with  perhaps  3.5-fold  greater 
values  of  DOCp^  on  the  outer  Nova  Scotian  shelf 
than  at  site  A  (Table  4) .  Recall  that  the  aircraft 
estimates  of  Gelbstoff  fluorescence  also  suggested 
a  fourfold  winter  gradient  of  DOC  between  coastal 
and  oceanic  waters  of  the  MAB  [Topliss  et  al .  , 
1989].  Furthermore,  previous  shipboard  observa¬ 
tions  of  DOM  fluorescence  indicate  a  threefold  to 
sixfold  gradient  of  stocks  between  basin  waters  at 
60°N,  20°W  and  either  the  East  Greenland  Current, 
or  shelf  waters  just  west  of  Scotland  [Duursma, 
1965]  . 

5.  Conclusions 

What  then  are  the  consequences  of  our  many 
assumptions?  We  begin  with  the  implications  of  the 
depth- Integrated  DOCp^  transports.  The  net  DOCpp 
transfer  of  the  whole  water  column  sums  to  a 
southward  export  across  36°25'N  of  16.87  x  10®  kg 
DOCpt  s'^,  or  1.4  X  10®  kmol  DOCpt  s'®,  in  the  case 
of  the  largest  value  of  shelf  DOCpt  stocks 
(Table  1) .  The  same  flow  fields  and  the  observed 
oxygen  concentrations  along  36'’25'N  yield  a  net 
southward  flux  of  2.9  x  10®  kmol  O2  s'®  [Rintoul, 
1988] . 

If  we  assume  a  molar  O2/C  consumption  ratio  of 
1.3  for  bacterial  degradation  of  subsurface  carbon 
substrates  [Redfield  et  al.,  1963],  the  ratio  of 
2.1  for  the  above  02/DOCpt  net  fluxes  suggests  that 
63%  of  the  AOU  demands  across  36'’25'N  would  be  met 
by  sinking  DOC,  and  37%  by  falling  POC,  if  all  of 
the  O2  were  consumed.  Furthermore,  a  DOC/DON  ratio 
of  10  from  dry  combustion  [Gordon  and  Sutcliffe , 
1973] ,  which  is  an  intermediate  value  between  that 
of  5  [Sugimura  and  Suzuki,  1988]  and  15  [Benner  et 
al.,  1992]  from  platinum  catalysts,  suggests  a  net 
equatorward  flux  of  1.7  x  10®  kg  DON  s'®.  This 
fortuitously  balances  the  postulated  poleward 
export  of  1.7  x  10®  kg  NO3  s'®  across  36°25'N 
[Rintoul,  1988]  of  unutilized  nitrate  within  the 
Gulf  Stream. 

The  smaller,  fourfold  gradient  of  DOCp^  between 
nearshore  and  offshore  regions  suggests  instead  a 
net  southward  flux  of  only  3.85  x  10®  kg  DOCpt 
(Table  1),  or  0.32  x  10®  kmol  DOCpt  ^^d  an 
02/DOCpt  ratio  of  9.1.  With  this  second  scenario 
and  the  Redfield  ratio,  86%  of  the  oxidizable 
carbon  substrate  may  be  in  particulate  form,  with 
just  14%  as  DOC.  Using  the  same  DOC/DON  ratio,  0.4 
X  10®  kg  DON  s'®  may  cross  36°25'N  towards  the 
equator,  offsetting  only  24%  of  the  poleward  flux 
of  NO3  in  this  case. 

If  we  choose  a  larger  C/N  ratio,  say,  15  [Benner 
et  al. ,  1992]  or  25  [Jackson  and  Williams,  1985], 
then  just  9-15%  of  the  northward  flux  of  Gulf 
Stream  nitrate  would  be  returned  as  DON  in  the 
second  scenario.  Larger  values  of  shelf  DOM,  or 
greater  oceanic  rates  of  photos)mthesis  and  excre¬ 


tion,  would  then  be  required  to  support  Rintoul  and 
Wunsch' s  [1991]  h)q)othesis.  In  contrast,  a  zero 
seaward  gradient  of  DOCpt  leads  to  a  net  poleward 
transfer  of  8.88  x  10®  kg  DOCpt  Iri  the  third 
case  (Table  1) ,  such  that  no  DON  would  move 
equatorward. 

Even  if  the  DOCpt  estimates  are  artifacts,  the 
net  DOCty  fluxes  show  a  similar  pattern  (Table  1) . 
Because  of  the  different  DOCpt/DOCu^  ratios  used  for 
shelf  and  oceanic  waters  (1.85  versus  3.54),  the 
second  DOCyv  scenario  reflects  an  eightfold  gradi¬ 
ent  of  D0Ct„,  while  the  third  case  poses  a  twofold 
gradient.  In  this  third  scenario,  no  DOCp^,  nor 
DON,  are  exported  south  from  the  northern  North 
Atlantic  (Table  1) . 

In  contrast,  the  second  case  yields  a  net 
southward  export  of  5.48  x  10®  kg  DOC^^  s'® 
(Table  1),  with  an  02/D0Cuv  ratio  of  6.3.  This 
calculation  suggests  that  79%  of  the  AOU  demands 
might  be  met  by  POC,  and  21%  by  DOCp.,,  similar  to 
Ogura’ s  [1970]  conclusions.  Only  32%  of  the 
proposed  northward  influx  of  NO3  would  then  be 
returned  south  as  DON,  leaving  open  the  possibility 
of  an  equatorward  flux  of  ammonium,  if  nitrifica¬ 
tion  is  slow  and  N2  is  not  lost  to  the  atmosphere 
north  of  36°25’N. 

Further  quantification  of  our  present  estimates 
are  confounded,  however,  by  seasonal  and  inter¬ 
annual  variations  of  the  sinking  processes,  by 
uncertainties  in  the  amount  of  primary  production 
and  excretion  of  DOC,  and  by  accurate  assessment  of 
spatio-temporal  changes  of  DOM  stocks.  Reduced 
subduction  beneath  the  Gulf  Stream,  increased 
coastal  photosynthesis,  reduced  CDOC/DOC  ratios, 
and  increased  release  of  DOC  within  coastal  habi¬ 
tats  from  the  marine  food  web  and/or  anthropogenic 
loadings  would  all  favor  (1)  larger  stocks  of  shelf 
DOC,  and  (2)  a  larger  role  of  DOC  in  oxygen  con¬ 
sumption  within  the  deep  sea.  The  opposite  trends 
of  these  processes  would  lead  to  smaller  stocks  and 
roles  of  shelf  DOC  in  global  carbon  cycles . 

We  do  conclude  that  some  shelf  export  of  DOC, 
with  a  positive  gradient  between  coastal  and 
oceanic  stocks,  as  well  as  falling  particles,  are 
required  to  balance  carbon,  nitrogen,  and  oxygen 
budgets  of  the  North  Atlantic,  based  on  one  snap¬ 
shot  of  hydrographic  properties  [Roemmich  and 
Wunsch,  1985] .  Seasonal  subduction,  cabelling,  and 
convection  of  surface  waters  may  initiate  annual 
storage,  in  sinking  waters  of  the  North  Atlantic, 
of  perhaps  0.9-1. 2  x  10®®  g  DOC  yr'®  (Table  1). 
Depending  upon  the  different  estimates  of  primary 
production  and  DOC  stocks  of  coastal  waters ,  the 
amount  of  DOM  supplied  by  rivers  and  the  shelf  food 
web  may  constitute  29-43%  of  the  proposed  basin¬ 
wide  sequestration  of  carbon  and  nitrogen. 

The  next  generation  of  biophysical  models 
involving  satellite  color  imagery  must  improve  upon 
more  primitive  ones  [Walsh  et  al. ,  1987;  Wroblewski 
et  al.,  1988;  Wroblewski,  1989;  Ishizaka,  1990; 
Platt  et  al.,  1991;  Gregg  and  Walsh,  1992]  by 

(1)  first  distinguishing  between  the  dissolved  and 
particulate  components  of  the  radiance  signal,  and 

(2)  by  providing  resolution  of  the  fates  of  the 
terrestrial  and  marine  sources  of  CDOC.  Future 
bio-optical  research  must  now  also  focus  on  coastal 
waters  of  high  and  middle  latitudes ,  as  well  as 
open  ocean  regimes,  to  provide  eventual  validation 
data,  from  follow-on  satellite  and  aircraft  color 
sensors,  for  our  crude  mathematical  hypotheses  of 
carbon  and  nitrogen  cycling  in  the  North  Atlantic. 
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Model  for  the  interpretation 
of  hyperspectrai  remote-sensing  refiectance 


Zhongping  Lee,  Kendall  L.  Carder,  Steve  K.  Hawes,  Robert  G.  Steward, 
Thomas  G.  Peacock,  and  Curtiss  O.  Davis 


Remote-sensing  reflectance  is  easier  to  interpret  for  the  open  ocean  than  for  coastal  regions  because  the 
optical  signals  are  highly  coupled  to  the  phytoplankton  (e.g.,  chlorophyll)  concentrations.  For  estuarine 
or  coastal  waters,  variable  terrigenous  colored  dissolved  organic  matter  (CDOM),  suspended  sediments, 
and  bottom  reflectance,  all  factors  that  do  not  covary  with  the  pigment  concentration,  confound  data 
interpretation.  In  this  research,  remote-sensing  reflectance  models  are  suggested  for  coastal  waters,  to 
which  contributions  that  are  due  to  bottom  reflectance,  CDOM  fluorescence,  and  water  Raman  scattering 
are  included.  Through  the  use  of  two  parameters  to  model  the  combination  of  the  backscattering 
coefficient  and  the  Q  factor,  excellent  agreement  was  achieved  between  the  measured  and  modeled 
remote-sensing  reflectance  for  waters  from  the  West  Florida  Shelf  to  the  Mississippi  River  plume.  These 
waters  cover  a  range  of  chlorophyll  of  0.2-40  mg/m^  and  gelbstoff  absorption  at  440  nm  from  0.02-0.4 
m"k  Data  with  a  spectral  resolution  of  10  nm  or  better,  which  is  consistent  with  that  provided  by  the 
airborne  visible  and  infrared  imaging  spectrometer  (AVIRIS)  and  spacecraft  spectrometers,  were  used  in 
the  model  evaluation. 


Introduction 

The  use  of  the  power  law  of  spectral-radiance  ratios^-^ 
to  measure  pigment  concentrations  requires  that  the 
water-leaving  radiance  be  largely  determined  by  varia¬ 
tions  in  the  pigment  concentration,  with  all  other 
optical  constituents  covarying  with  this  quantity. 
The  method  works  quite  well  for  the  open  ocean,  or 
case  1  waters,^  in  part  because  the  water-leaving 
radiance  of  open  ocean  waters  is  hardly  affected  by 
bottom  reflectance,  land  runoff,  or  suspended  sedi¬ 
ments.  Although  aeolian  dust  may  be  carried  by 
winds  to  the  open  ocean,^  the  dominant  effect  of  the 
particulates  may  still  derive  from  ph3d;oplankton.® 
The  power-law  approach  can  be  much  less  accurate 
for  estuarine  and  coastal  areas,®  however,  because 
many  of  the  optical  constituents  are  independent  of 
phytoplankton  concentrations.  In  these  areas,  the 
water-leaving  radiance  includes  not  only  parts  that 
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are  due  to  elastic  scattering  by  water  molecules, 
phytoplankton  detritus,  suspended  particulates,  and 
bottom  reflectance,  but  also  parts  that  are  due  to 
inelastic  scattering  of  colored  dissolved  organic  mat¬ 
ter  (CDOM)  fluorescence  and  water  Raman  scattering. 
Thus,  changes  in  ocean  color  resulting  from  sus¬ 
pended  sediments  or  dissolved  organic  matter  may  be 
falsely  interpreted  as  changes  in  pigment  concentra- 
tion.®-’' 

An  approach  to  address  these  problems  is  to  mea¬ 
sure  the  light  field  and  analytically  separate  the 
different  spectral  contributors.  Optical  models  have 
been  developed  for  the  subsurface  irradiance  reflec¬ 
tance,’^’®  but  satellites  measure  the  radiance  leaving 
the  water  surface.  The  water-leaving  radiance  is 
governed  by  two  distinct  parts:  the  solar  input  and 
in-water  properties.  When  the  remote-sensing  reflec¬ 
tance  (i?rs)  is  defined  as  the  ratio  of  the  water-leaving 
radiance  L„(A)  to  the  above-surface  downwelling  irra¬ 
diance  £d(0+,  \),  J?rs  will  be  independent  of  the  inten¬ 
sity  of  the  solar  input.  Models  have  been  suggested 
by  Carder  and  Steward,®  Gordon  et  and  Peacock 
et  al.  to  explain  the  measured  R^s,  but  in  these  works 
no  contributions  from  CDOM  fluorescence,  water 
Raman  reflectance,  or  bottom  reflectance  were  in¬ 
cluded.  Also,  in  these  works®’^^  an  arbitrary  Q  fac- 
tor^^  was  used. 

For  water-depth  measurements^®’^^  or  bottom- 
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feature  mapping,  the  diffuse  attenuation  coeffi¬ 
cient^^  or  an  unclear  effective  attenuation  coeffi¬ 
cient^®  is  usually  used.  For  improvement  of  the 
interpretation  of  such  measurements,  an  explicit 
expression  is  necessary  for  attenuation  of  the  reflec¬ 
tance  term  that  represents  bottom  reflectance.  In 
this  study,  hyperspectral  remote-sensing  reflectance 
for  waters  from  the  West  Florida  Shelf  to  the  Missis¬ 
sippi  River  plume  was  measured  and  modeled  with  a 
derived  Q  factor  and  the  addition  of  bottom  reflec¬ 
tance,  CDOM  fluorescence,  and  water  Raman  scatter¬ 
ing. 


In  the  interpretation  of  the  measured  on  the  left 
side  of  Eq.  (3),  each  component  on  the  right  side  of 
Eq.  (3)  is  expressed  by  the  optical  properties  of  the 
water. 

A.  Remote-Sensing  Reflectance  of  the  Water  Column 

For  a  homogeneous  water  body  consider  a  wavelength- 
independent  factor  I  as  the  influence  of  the  air-sea 
interface  on  water-leaving  radiance.  Then  the 
which  is  due  to  eleastic  scattering  in  the  water 
column,  can  be  described  in  terms  of  values  just  below 
the  interface  as 


Theory 

The  upwelling  radiance  leaving  the  ocean  is  a  compli¬ 
cated  mix  of  signals  caused  by  many  components. 
The  major  contributions  arise  from  the  following: 
absorption  by  molecules  and  particulates,  elastic  scat¬ 
tering  by  molecules  and  particulates,  and  bottom 
reflectance  in  shallow  waters.  Inelastic  scattering 
processes  (e.g.,  water  Raman  scattering  and  fluores¬ 
cence  of  CDOM,  chlorophyll,  and  phycoerythrin)  are 
also  contributors.  Works  by  Carder  and  Steward® 
and  Gordon^®  dealing  with  chlorophyll  a  (chi  a) 
fluorescence  have  been  reported,  but,  because  the 
fluorescence  efficiency  varies  by  an  order  of  magni¬ 
tude,®  this  term  is  not  considered  in  the  present 
model.  As  peak  chlorophyll  fluorescence  occurs  in  a 
narrow  band  centered  around  685  nm,^'^  its  absence 
from  the  model  is  clearly  seen  when  one  compares 
values  of  measured  and  modeled  i?rs  (685)  at  higher 
chlorophyll  concentrations.  Perturbations  between 
modeled  and  measured  curves  at  ~  580  nm  and  ~  685 
nm  will  be  considered  in  the  discussion  relative  to 
fluorescence  that  results  from  phycoerythrin  and  chi 
a,  respectively. 

It  is  assumed  that  the  water-leaving  radiance  Lu,(X) 
is  dominated  by  the  following  four  components: 
elastic  scattering  from  molecules  and  particles  L„“’(\), 
bottom  reflectance  LJ’{X),  CDOM  fluorescence  LJ{\), 
and  water  Raman  scattering  LJ^{X).  It  is  also  as¬ 
sumed  that  to  the  first  order  (single  scattering  and 
quasi-single  scattering^®)  the  water-leaving  radiance 
can  be  expressed  as 

LJ\)  =  L„“(\)  +  +  LJ(\)  +  L/(X).  (1) 

The  symbols  and  definitions  used  in  this  paper  are 
summarized  in  Table  1.  Wavelength  dependence  is 
included  here,  but  for  convenience  it  is  not  included 
in  the  following  discussions  except  when  necessary 
for  clarity. 

Remote-sensing  reflectance  is  defined  as 


R 


W 

rs 


(4) 


where  7  =  t+t_/n^.  For  a  zenith  Sun,  a  nadir¬ 
viewing  instrument,  and  a  calm  surface,  t  ~  0.98  and 
I  ~  0.533  (Ref.  12)  because  n  (1.341)  varies  only 
slightly  when  the  salinity  changes,  i®  For  larger  so¬ 
lar  zenith  angles  and  foam-covered  seas,  t  will  be 
lower.  In  our  work,  we  specifically  avoid  foam- 
covered  seas. 

The  subsurface-irradiance  reflectance  from  the  wa¬ 
ter  column  is  defined  as  the  ratio  of  the  subsurface 
upwelling  irradiance  R„“'(0~)  to  the  subsurface  down- 
welling  irradiance  Rd(0~)  by 


R“' 


Rr(0~) 

E,{0-)  ’ 


(5) 


and  Austin^®  has  related  ^^“'(0")  and  Lu“'(0“)  through 
the  Q  factor  by 


^  LA0-)  ’ 

so  Rrs'"  can  be  expressed  as 

7 

'“'rs  ' 


(6) 

(7) 


For  irradiance  reflectance  R",  Gordon  et  al.^°  devel¬ 
oped  a  series  relation  with  the  Monte-Carlo  method 
by 


R“  =  S 

m=0 


\a  +  65 

This  equation  was  simplified®^  ®®  to 

h 


R“ 


0.33 


(8) 


(9) 


R 


rs 


Ew 

EdiO^)  ■ 


(2) 


Breaking  this  equation  into  contributions  from  the 
various  mechanisms  listed  in  Eq.  (1),  we  have 

Ets  =  -Rrs“  +  EJ  +  EJ  +  Rrs^.  (3) 


for  values  of  bb/a  up  to  ~0.25.  The  constant  0.33 
actually  varies  slightly  with  the  solar  zenith  angle.®®’®^ 
Because  this  paper  deals  with  remote-sensing  reflec¬ 
tance,  which  is  less  influenced  by  the  fraction  of 
forward  scatter  that  upwells  at  solar  zenith  angle 
>0°,  Rfs  is  not  as  sensitive  to  the  Sun  angle,  as  is 
irradiance  reflectance.  Thus,  we  retain  0.33  as  a 
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Table  1 .  Definitions  and  Units  of  Variables' 


Variable 

Units“ 

Definition 

a 

m-i 

Total  absorption  coefficient:  a^,  +  Ug  +  Op 

m-i 

Gelbstoff-absorption  coefficient 

dp 

Particle-absorption  coefficient 

div 

m-i 

Pure  water  absorption  coefficient 

b 

m-i 

Scattering  coefficient 

hb 

m-i 

Backscattering  coefficient 

^bm 

m-i 

Backscattering  coefficient  of  molecules 

^bp 

m-i 

Backscattering  coefficient  of  particles 

m-i 

Water  Raman  scattering  coefficient 

C 

m-i 

Beam  attenuation  coefficient:  a  +  b 

Downwelling  distribution  function  just  below  the  surface 

Vertically  averaged  downwelling  distribution  function 

Ed 

Wm-2 

Downwelling  irradiance 

Eo 

Wm-2 

Total  scalar  irradiance 

Eoi 

Wm-2 

Downwelling  scalar  irradiance 

■fi'ou 

Wm-2 

Upwelling  scalar  irradiance 

Wm-2 

Upwelling  irradiance  from  water  column  only 

H 

m 

Water  depth 

k 

Radiance  attenuation  coefficient 

j 

rad 

Subsurface  solar  zenith  angle 

Kd 

Downwelling  diffuse  attenuation  coefficient 

Ku 

Upwelling  diffuse  attenuation  coefficient 

Lu'^ 

W  m"^  sr“^ 

Upwelling  radiance  from  water  column  (elastic  only) 

Liw 

W  m~^  sr'^ 

Total  water-leaving  radiance 

LJ 

W  m'2  sr“i 

Water-leaving  radiance  from  bottom  reflectance 

LJ 

W  m"2  sr"^ 

Water-leaving  radiance  from  CDOM  fluorescence 

T  R 

Wm-2sr-i 

Water-leaving  radiance  from  water  Raman  scattering 

T  w 

W  sr'^ 

Water  leaving  radiance  from  water  column  (elastic  only) 

n 

Refractive  index  of  water 

Q 

sr 

Ratio  of  irradiance  to  radiance 

R 

Irradiance  reflectance 

■^rs 

sr“^ 

Remote-sensing  reflectance 

sr"^ 

Remote-sensing  reflectance  from  bottom  reflectance 

RJ 

sr"i 

Remote-sensing  reflectance  from  CDOM  fluorescence 

RrJ^ 

sr"^ 

Sum  of  Rr/  and  Rrs^ 

R  ^ 

iVrs 

sr"i 

Remote-sensing  reflectance  from  water  Raman 

sr'i 

Remote-sensing  reflectance  from  water  column  (elastic  only) 

R<^ 

Irradiance  reflectance  from  water  column  (elastic  only) 

Air-sea  surface  transmittance  for  upwelling  radiant  flux 

Air-sea  surface  transmittance  for  downwelling  radiant  flux 

Tie 

Wsr-i 

Intensity  of  the  inelastically  scattered  light 

yb 

Exponent  for  particle  backscattering  coefficient 

yQ 

Exponent  for  particle  Q  factor 

a 

rad 

Scattering  angle 

p 

m*'  sr“^ 

Volume  scattering  function 

<|) 

rad 

Azimuth  angle 

7 

Irradiance  ratio  of  skylight  to  sunlight 

K 

m“^ 

Quasi-diffuse  attenuation  coefficient:  a  +  bi, 

Xx 

nm 

Excitation  wavelength  for  inelastic  scattering 

Quantum  efficiency  of  CDOM  fluorescence 

>!< 

m“*  nm“^ 

Inelastic  scattering  coefficient 

p 

Bottom  albedo 

0 

rad 

Zenith  angle 

“Blank  entries  denote  dimensionless  quantities. 


constant  and  embed  any  sun-angle  influence  within 
the  Q  factor  expression  (see  Appendix  A). 

The  total  backscattering  coefficient  65  includes  two 
components:  backscattering  by  molecules  6bm  and 
particulates  6bp-  The  total  absorption  coefficient  a 
includes  contributions  that  are  due  to  pure  seawater 
absorption  a^,  gelbstoff  or  CDOM  absorption  Og,  and 
particulate  absorption  Op.  Inserting  these  into  Eq. 


(9)  and  taking  I  ~  0.533,  we  can  write  Brs“’  as 


0.176 


^bm  ^ 


bp 


G,,,  ~h  Gtr  G„ 


Q 


(10) 


Eq.  (10)  pertains  to  optically  deep  water.  When 
optically  shallow  water  is  encountered,  scattering 
media  and  backscattered  signals  are  reduced  because 
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of  the  short  water  column.  We  consider  that  the 
subsurface  consists  of  two  parts  coming  from  two 
layers:  one  from  the  layer  above  the  bottom  and  one 
from  the  layer  below  the  bottom.  Then  the  subsur¬ 
face  coming  only  from  the  upper  layer  can  be 
obtained  by  a  reduction  of  the  optically  deep  expres¬ 
sion  by  an  amount  equivalent  to  the  contribution  of 
the  missing  water  column  below  depth  H.  Thus  for 
shallow  waters  with  depth  H  and  a  totally  absorbing 
bottom,  is  approximated  by 


Rr 


0.176  h\ 
a 


{Ka  +  ^TJdz 


(11) 


ing  function  for  inelastic  scattering  as 


Tijo.,  A.) 

E{KW’ 


(15) 


where  X)  is  the  intensity  of  the  inelastically 
scattered  light  at  scattering  angle  a,  dV  is  the  scatter¬ 
ing  volume,  E{'K^)  is  the  irradiance  of  the  excitation 
beam,  and  a  is  the  angle  between  the  excitation  beam 
and  the  output  photon  directions. 

If  Pie  is  considered  to  be  isotropic,  then  remote¬ 
sensing  reflectance  resulting  from  CDOM  fluores¬ 
cence  and  water  Raman  can  be  expressed  (see  Appen¬ 
dix  B  for  details)  as 


where  z  is  positive  downward  from  the  surface. 

If  we  define  the  quasi-diffuse  attenuation  coeffi¬ 
cient  as  K  =  a  -I-  feft,  then  Kd  ~  and  ~ 

Dd  and  are  the  distribution  functions  for  the 
downwelling  and  upwelling  light  fields,  and  D^jBci  ~ 
2  according  to  Gordon  Thus 

0.176 

^  [1  -  exp[-3{D^}K7f],  (12) 

where  {0^]  is  the  vertically  averaged  downwelling 
distribution  function  and  (D^)  =  1.08Z)d(0).2s  DjiO) 
is  the  downwelling  distribution  function  just  beneath 
the  surface,  and  0^(0)  =  l/cos(y),  with  an  error  of 
less  than  3%.^® 

B.  Remote-Sensing  Reflectance  Resulting  from  Bottom 
Reflectance  firs'’ 

Assume  that  the  bottom  is  a  Lambertian  reflector 
with  bottom  albedo  p,  then  R^J’  can  be  approximated 
as 


,  0.533 

RJ  - - P  exp[-{[D,]K  +  k)H],  (13) 

TT 

where  k  is  the  effective  attenuation  coefficient  for  the 
radiance  from  an  extended  Lambertian  source.  How 
k  relates  to  the  quasi-diffuse  attenuation  coefficient  k 
is  not  well  understood.  Heuristically,  it  should  be  a 
value  between  the  beam  attenuation  coefficient  c  and 
the  quasi-diffuse  attenuation  coefficient  k.  Taken 
from  the  Monte-Carlo  simulations  for  a  totally  diffuse 
light  source,25  k  is  approximately  1.4k  to  1.7k  for  kH 
in  the  range  of  0. 5-4.0.  As  an  average  in  this  work, 
k  =  1.5k,  as  is  used  by  Marshall  and  Smith.^®  Then 
Eq.  (13)  becomes 

RJ  =  0.17p  exp[-(1.5  +  {I>d})KH].  (14) 


C.  Remote-Sensing  Reflectance  Resulting  from  CDOM 
Fluorescence  and  Water  Raman  Scattering  and  firs” 

In  general,  these  terms  are  due  to  inelastic  scattering 
(indicated  by  the  subscript  ie)  by  CDOM  molecules 
and  water  molecules.  We  defne  the  volume  scatter¬ 


RJ 


=  0.072  I 


kx  (^g(K)Rd(O  !  K 


X,  ^  [2a(X)  +  a(Xj]Ed(0  ,  X) 

I  I 

exp  -s  In - 


R  ”  =  0.072 


A 


dX,. 


[2o(X)  -I-  a(Xj]Err(0  ,  X) 


(16) 

(17) 


Field  Measurements 

From  1990  to  1993,  measurements  of  optical  proper¬ 
ties  for  case  1  and  case  2  waters,  which  include  waters 
from  the  West  Florida  Shelf  to  the  mouth  of  Missis¬ 
sippi  River,  were  taken.  Case  2  waters  are  those 
that  contain  optical  materials  that  are  not  derived 
from  phytoplankton,  in  addition  to  phytoplankton 
and  phytoplankton-derived  materials.®  Table  2  sum¬ 
marizes  the  field  data  for  selected  stations.  Figure  1 
shows  the  station  locations  in  the  Gulf  of  Mexico. 
For  each  station,  remote-sensing  reflectance  Rrs  and 
surface-water  particulate  absorption  were  mea¬ 
sured.  For  the  1993  stations,  a  long-path  (50  or  100 
cm)  spectrophotometer  was  used  to  measure  Ug. 

For  Rrs  we  directly  measured  the  upwelling  radi¬ 
ance  above  the  sea  surface  and  downwelling  sky 
radiance  using  a  Spectron  Engineering  spectral  radi¬ 
ometer  (Model  SE-590)  following  the  method  of 
Carder  We  measured  downwelling  irradiance 

above  the  sea  surface  with  the  SE-590  by  viewing  a 
Spectralon  diffuse-reflection  calibration  panel.  Re¬ 
mote-sensing  reflectance  values  were  determined  by 
removal  of  the  reflected  skylight  from  the  upwelling- 
radiance  values®’®'^  and  division  of  the  result  by  the 
downwelling-irradiance  values.  For  each  station, 
three  sets  of  measurements  were  taken,  and  an 
averaged  Rrs  spectrum  was  derived,  with  coefficients 
of  variation  much  less  than  5%. 

For  Up,  we  measured  a  surface-water  sample  imme¬ 
diately  after  its  collection  following  the  method  devel¬ 
oped  by  Mitchell  and  Kiefer.®®  Briefly,  for  each 
water  sample,  ~  1000  mL  (with  variations  according 
to  the  clarity  of  the  water  sample)  were  filtered 
through  Whatman-type  GF/F  glass-fiber  filters. 
Hyperspectral  optical  densities  of  the  sample  pad  and 
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Table  2.  Station  Locations  and  Water  Depth 


Station 

Latitude 

Longitude 

Time 

Date 

Bottom  Depth 

Mod.  Bottom  Depth 

STOl 

27°27'  N 

82‘"55'  W 

10.5“ 

4  Mar  1990 

~  14  m 

13.7  m 

ST02 

27°20'  N 

83°03' W 

13.0“ 

4  Mar  1990 

~25  m 

25  m 

ST03 

27n2'N 

83°11'W 

14.9“ 

4  Mar  1990 

~  35  m 

36  m 

ST08 

28°48'  N 

91°30'W 

08.5* 

12  Apr  1993 

OD'* 

OD 

STIO 

28°15' N 

91°30'W 

14.0* 

12  Apr  1993 

OD 

OD 

ST27 

29°32'-N 

85°47'  W 

09.2* 

19  Apr  1993 

OD 

OD 

ST12 

28°52'  N 

89°33'  W 

10.8“ 

5  Jun  1993 

OD 

OD 

ST14 

28°48'  N 

90°02'  W 

16.1“ 

5  Jun  1993 

OD 

OD 

ST19 

27°34'  N 

83°20’  W 

09.5“ 

8  Jun  1993 

~33  m 

35  m 

“Eastern  standard  time. 

>  ^Central  daylight  time. 

“Eastern  daylight  time, 
“^ater  is  optically  deep  (OD). 


a  wet  blank  pad  were  measured  with  the  Spectron. 
The  optical-path-elongation  factor  p  was  calculated 
with  Eq.  (2)  in  Bricaud  and  Stramski.^®  The  optical 
density  measurements  were  repeated  three  times 
with  no  significant  variation  noted  among  them. 

At  the  1990  stations,  gelbstoflf  absorption  was 
derived  from  surface-layer  values  determined  with 
a  Biospherical  Instruments  MER-1048,  through  the 
use  of  the  expression  Ug  =  cos(y)  -  -  a^.  At 

the  1993  stations,  was  measured  with  50-cm  or 
100-cm  path-length  instruments,  respectively,  after 
the  sample  was  filtered  through  0.2-\im  pore-diam¬ 
eter  Gelman  Supor-200  filters. 

Model 

For  the  modeling  of  measured  spectra,  values  for 
Oui  and  6bm  were  already  known,®®  Up  was  measured, 
and  Ug  was  measured  or  derived  from  spectra. 
What  needs  to  he  considered  is  how  Q,  p,  H,  t|,  s, 
\o,  and  CT  change  for  different  environments. 


Because  Raman  is  a  type  of  molecular  scattering, 
is  considered  to  have  a  wavelength  dependence 
similar  to  that  of  the  water-molecule  scattering  coeffi- 
cient,®!  i.e.,  a  function  of  X,"^.  Because  6j;(488)  = 
2.6  X  10"*  m~i  (Marshall  and  Smith®®),  thus  6ie(X,.)  = 
2.6  X  10~^(488/X,c)‘‘,  and  the  frequency  shift  for  water 
Raman  scattering  was  fixed  at  3350  cm“i  as  an 
average  from  Collins  et  al.^^  The  incoming,  total 
downwelling-irradiance  spectrum  was  measured  with 
a  Licor  1800  spectral  irradiance  meter®^  from  300- 
850  nm.  Then  it  is  straightforward  to  calculate 
with  Eq.  (17)  when  the  total  absorption  spectrum  is 
known. 


As  can  be  seen  from  Eq.  (16),  there  are  at  least  four 
variables  [ti(Xi),  s,  Xq,  and  a]  needed  to  calculate  the 
remote-sensing  reflectance  that  results  from  CDOM 


Longitude 


Fig.  1.  Station  locations  in  the  Gulf  of  Mexico. 
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fluorescence  when  the  total  absorption  coefficient  is 
known.  From  lab  measurements  of  CDOM  fluores¬ 
cence  for  the  West  Florida  Shelf  experiments,^^  the 
quantum  efficiency  'ri(\j;)  was  between  ~0.5%  and 
~1.5%,  and  it  was  generally  rather  constant  for 
different  excitation  wavelengths.  For  the  log-nor¬ 
mal  expression  in  R^/,  the  shape  factor  s  was  10,  X.o  ~ 
(.95Xj:  -  45),  and  o  ~  (195  -  of  which  were 

quite  constant  for  the  different  stations.^^  Based  on 
these  measurements,  q  =  1.0%  was  used  for  all  the 
stations  except  ST03,  where  q  =  1.5%  was  measured. 
At  ST12,  ST14,  and  STOl  where  significant  terrig¬ 
enous  CDOM  was  present,  q  =  0.5%  was  used. 


If  we  consider  the  upwelling  radiance  of  the  water 
column  Lu‘"  as  consisting  of  two  parts — one  L„'"  that 
is  due  to  scattering  by  molecules  and  one  L/  that  is 
due  to  scattering  by  particles — with  the  assumption 
that  to  the  first  order  (for  single  scattering  and 
quasi-single  scattering^®)  the  sum  of  and  L/  gives 
Lu’",  then  Eq.  (10)  can  be  written  as 


Rr 


0176  /febm  M 

\Qm  Qpl 


(18) 


in  which  and  Qp  are  the  Q  factors  for  molecules 
and  particles,  respectively,  and  are  defined  as 


This  value  depends  not  only  on  the  optical  properties 
of  the  water  body,  but  also  on  the  water  depth  and  the 
bottom  albedo.  In  the  modeling  work,  the  water 
depth  was  based  on  the  Provisional  Chart  No.  1003®^ 
for  the  Gulf  Coast,  and  the  bottom  albedo  was  based 
on  earlier  measurements  of  bottom  samples  from  the 
region  that  had  near-shore  values  of  0. 1-0.2  (used  for 
STOl)  and  offshore  values  from  0.4-0.5  (used  for 
ST02  and  ST03).  Figure  7  (below)  shows  examples 
of  those  albedo  spectra.  The  quasi-diffuse  attenua¬ 
tion  coefficient  k  is  assumed  to  be  equal  to  the 
absorption  a.  For  ST19,  an  albedo  value  of  0.1  is 
required,  which  suggests  that  the  bottom  might 
contain  more  heavy  minerals  or  grass  at  that  site. 
Direct  bottom-albedo  measurements  are  lacking  at 
individual  stations  and  are  needed  for  a  wide  variety 
of  bottom  types. 


R  ^ 

"rs 

When  using  Eq.  (10)  to  model  the  measured  Rrs"",  we 
need  to  know  Cw,  Og,  Op,  6bm,  hbp,  and  Q.  Values  of  a„, 
and  febm  are  already  known.®®  Wfiien  Ug  and  Op  are 
measured,  only  6bp  and  Q  for  different  water  bodies 
emd  solar  zenith  angles  are  required.  The  particu¬ 
late  backscattering  coefficient  6bp  has  been  considered 
to  be  a  spectral  function  of  for  offshore  waters  or 
to  be  spectrally  constant  for  near-shore  waters.^  ®® 
For  the  factor  Q,  however,  only  a  few  measure¬ 
ments  exist,  and  its  values  have  been  reported  from 
3.2  to  12.®  Q  has  been  taken  to  be  approximately  4.7 
and  spectrally  constant  from  440-550  nm,®®  although 
Kirk®'^  gives  Q  as  ~4.9,  and  Gordon  et  suggest 
a  value  of  ~3.4.  For  many  studies,  Q  is  often 
arbitrarily  chosen  as  a  spectral  constant.®’^^’^®  From 
measurements  by  Davis,®®  however.  Carder  et  al.^ 
found  that  Q  is  not  spectrally  constant  for  the  1990 
stations,  and  there  was  a  trend  for  Q  to  increase  with 
wavelength  (an  inverse  trend  compared  with  6bp)- 
Recently,  Morel  and  Gentili®®  published  Monte-C8U"lo 
simulations  of  Q  for  a  variety  of  water  types,  but  they 
provide  no  explicit  expression  for  Q  as  a  function  of  a 
and  6b  or  p.  So,  for  modeling  for  a  region  where 
spectral  6bp  and  Q  need  to  be  considered,  at  least  four 
parameters  are  needed  based  on  Eq.  (10):  two  for 
6bp(X)  and  two  for  Q(\). 


Qm 


L„-(0-)  ’ 


V(o-)' 


(19) 


For  Qm  to  the  first  order  an  estimate  can  be  made 
based  on  the  phase  function  and  illumination  geom¬ 
etry.  For  a  given  illumination  geometry,  the  shape 
of  the  radiance  distribution  within  the  water  is 
determined  primarily  by  the  volume  scattering  func¬ 
tion  through  single  scattering.  For  example,  Gor¬ 
don®®  suggested  that  a  single-scattering  approxima¬ 
tion  can  be  used  to  specify  the  variation  of  R  with  the 
solar  zenith  angle,  and  Krk®®  used  single  scattering 
to  describe  the  average  cosine.  We  combined  the 
approach  used  by  Jerlov^®  to  provide  an  estimation  of 
radiance  and  irradiance  with  Sun  angle  and  depth, 
Austin’s  definition  of  the  Q  factor,^®  and  the  volume 
scattering  function  of  water  molecules  given  by  Mo- 
reF^  to  calculate  the  Qm  for  sunlight  (see  Appendix  A 
for  details).  The  results  can  be  approximated  with 
the  following  simple  function: 

Qm^'^^U)  »  5.92  -  3.05  cos(i).  (20) 

With  the  assumption  that  the  Qm  attributable  to 
skylight  is  approximately  3.14,  the  effective  Qm  for  a 
mixture  of  sunlight  and  skylight  is  given  by  (Appen¬ 
dix  A) 


Qm  — 


1  -I-  7(\) 


1  -I-  y(k) 


Qr^) 

3.14 


r.Qm^'^^U), 


(21) 


if  we  define  7(X)  =  and  calculate  7(X) 

using  the  model  developed  by  Gregg  and  Carder.'*® 
Model  results  of  Qm®™(j)  centered  around  3.3  for 
environments  studied  in  this  contribution  and  are 
shown  in  Table  3,  and  the  calculated  is  consis¬ 
tent  with  Morel  and  Gentili’s  Monte-Carlo  results.®® 
Because  we  know  the  volume  scattering  function 
for  neither  the  total  water  sample  nor  the  particles, 
6bp  and  Qp  cannot  be  independently  estimated. 
However,  because  6bp  has  been  considered  a  function 
of  6bp(400)(400/X)>®  as  in  Smith  and  Baker,®®  we  may 
also  consider  Qp  to  be  a  function  of  Qp(400)(X/400)y®. 
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Table  3.  Parameters  for  Each  Station  ■< 


Station 

j 

X 

Y 

ag(440) 

ap(440) 

O^mea(440) 

[chi  a]  (mg/m^) 

STOl 

35° 

3.3 

0.0090 

1.5 

0.082 

0.040 

0.045 

1.05 

ST02 

26° 

3.2 

0.0020 

2.4 

0.042 

0.036 

0.035 

0.61 

ST03 

27° 

3.2 

0.0010 

2.4 

0.034 

0.028 

0.026 

0.70 

ST08 

43° 

3.7 

0.0062 

0.3 

0.31 

0.28 

0.295 

STIO 

17° 

3.0 

0.0009 

1.8 

0.059 

0.023 

0.021 

ST27 

CO 

3.5 

0.0011 

1.7 

0.078 

0.029 

0.034 

ST12 

o 

CO 

3.3 

0.0029 

0 

0.42 

1.40 

1.12 

38.58 

ST  14 

21° 

3.1 

0.0066 

0 

0.38 

1.21 

1.16 

20.26 

ST  19 

CO 

3.5 

0.00061 

1.9 

0.023 

0.022 

0.013 

0.22 

“Blank  entries  indicate  that  no  data  are  available. 


Then  6bp/ Qp  can  be  combined  and  modeled  as 

^  ^  6bp(400)  iyb+yq) 

%  "  Qp(400) 

=  Z(-)  .  (22) 

where  X  and  Y  are  two  unknowns  determined  for 
specific  particulate  suites  and  solar  illumination  situ¬ 
ations  if  Eq.  (18)  is  inverted. 

After  calculating  R^s^,  RJ,  and  RJ’,  only  X  and  Y 
remain  unknown.  The  modeled  7?^“  arid  the  re¬ 
sidual  of  jRps  -  -  Rr/  -  i2rs*  Were  matched  to 

derive  Zand  Y  with  a  predictor-corrector  approach  to 
modeling,  as  in  Carder  and  Steward.® 

Results  and  Discussion 

Using  the  methodology  described  above  (see  Model) 
we  modeled  R^s  for  case  1  and  case  2  waters,  which 
include  (1)  the  West  Florida  Shelf  waters,  with  shal¬ 
low,  gelbstoff-rich  coastal  waters,  and  (2)  Gulf  of 
Mexico  waters,  with  phytoplankton  blooms  in  the 
Mississippi  River  plume  (S  >  17%).  As  examples. 
Figs.  2-5  show  the  detailed  model  components  for  Rrs. 
and  Fig.  6  shows  the  results  of  all  the  listed  stations. 
Table  2  provides  the  station  locations  as  well  as  the 
measured  and  modeled  water  depths  of  the  shallow 
stations.  Table  3  lists  the  model  parameters  j, 


X,  Y,  cig(440)  and  ap(440)  along  with  the  measured 
values  of  ap(440)  and  the  chi  a  concentration  values 
for  each  station.  Table  4  details  the  fractional  contri¬ 
butions  that  Rrs'")  Rrs’^,  Rts^,  ^nd  R,.J’  make  to  the 
measured  i?rs  at  440  and  550  nm. 

It  can  be  seen  from  Figs.  2-6  that  excellent  fittings 
were  achieved  between  the  measured  and  modeled  Rj.^ 
for  all  the  stations  except  the  spectral  region  near  685 
nm,  where  chi  a  fluorescence  is  present  in  the  field 
data.  The  overall  averaged  difference  between  the 
measured  and  modeled  is  2%,  which  is  well  within 
the  measurement  accuracy.  The  that  is  required 
by  the  model  is  within  15%  of  the  measured  Op  except 
near  the  Mississippi  River  at  ST12  (25%),  with  the 
highest  chi  a  concentration  (38.6  mg/ m®).  The  aver¬ 
age  difference  between  the  measured  and  required  Op 
is  8.9%  (6.9%  when  ST12  is  excluded).  The  maxi¬ 
mum  15%  or  25%  difference  can  perhaps  be  explained 
by  the  accuracy  involved  in  the  method  of  Up  measure¬ 
ment  because  of  the  |3  factor,  which  varies  signifi¬ 
cantly  among  species.®®  This  may  be  especially  im¬ 
portant  for  ST  12,  which  was  near  the  Mississippi 
River  mouth  where  the  heavy  load  of  sediments  and 
minerals  might  cause  additional  uncertainty.  Also, 
the  influence  of  the  horizontal  and  vertical  structures 
of  the  waters  increases  for  mesotrophic-eutrophic 
waters,  so  patchiness  can  affect  accuracies  in  the 
more  hypertrophic  waters.  Finally,  the  low  signal 
obtained  for  the  upwelling  radiance  measurements  at 


wavelength  (nm) 

Fig.  2.  Measured  versus  modeled  for  ST27. 


wavelength  (nm) 

Fig.  3.  Measured  versus  modeled  for  ST14. 
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wavelength  (nm) 

Fig.  4.  Measured  versus  modeled  R^s  for  STO 1 . 

ST12  made  the  R^s  calculation  sensitive  to  corrections 
for  reflected  skylight. 

From  Table  3,  the  ratio  of  0^(440)  to  0^(440)  is 
highly  variable,  with  a  range  from  ~  0.3  to  ~  3.0,  and 
the  X  value  does  not  covary  with  the  pigment  concen¬ 
tration  of  chi  a  for  the  waters  studied.  This  illus¬ 
trates  that  the  model  works  well  over  a  wide  range  of 
conditions  and  also  suggests  why  the  power-law 
algorithm  does  not  work  well  for  coastal  waters. 
The  highest  X  value,  0.0090  m“i  sr'^,  was  at  the 
shallow,  mesotrophic  waters  at  STOl,  which  suggests 
a  high  influence  of  detritus  and  suspended  sediments. 
Brisk  northwesterly  winds  suspended  sediments  in 
the  shoal  regions  to  the  east  and  north  of  the  station, 
and  sediment  and  detritus  likely  were  transported  by 
the  ebb  tidal  currents  from  Tampa  Bay  to  the  study 
site.2'^ 

The  Y  values  were  generally  within  the  range  of 
0-2.4  for  the  waters  reported  here.  This  range 
might  be  interpreted  to  be  partially  due  to  6bp  and 
partially  due  to  Qp.  For  the  particle  backscattering 
coefficient  6bp  the  wavelength  exponent  is  in  the  range 
0-3.0  for  a  range  of  particle  sizes  (e.g.,  bacteria,'*® 
ph37toplankton  cells,*'*  and  coccoliths*®).  From  Eq. 
(22),  if  y  =  2.4  andjQ  =  1.0,  thenyb  =  1.4,  which  is 
within  the  0-3.0  range  reported  elsewhere. 

Model  results  at  ST03  suggest  relatively  higher 
CDOM  fluorescence  and  water  Raman  influences. 


400  500  600  700 


wavelength  (nm) 

Fig.  5.  Measured  versus  modeled  R^s  for  ST19. 
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wavelength  (nm) 

Fig.  6.  Measured  (solid  curve)  versus  modeled  (dashed  curve)  R^^ 
for  all  the  stations  except  STOl. 

because  a  higher  t|  (1.5%)  was  encountered.®®  This 
value  is  ~  3  times  greater  than  the  value  suggested  by 
Spitzer  and  Dirks*®  for  terrigenous  CDOM.  If  we 
exclude  this  station,  more  than  90%  of  the  water¬ 
leaving  radiance  is  accounted  for  by  the  sum  of  the 
elastic  scattering  from  molecules,  particles,  and  the 
bottom.  This  is  consistent  with  the  reports  of  Mar¬ 
shall  and  Smith®®  and  Stavn,*®  because  the  waters 
studied  in  their  reports  were  clearer.  Water  Raman 
scattering  makes  more  of  a  contribution  when  the 
water  is  clear,  and  i?rs^/^rs  will  typically  have  a 
higher  value  at  550  nm  than  at  440  nm. 

It  is  interesting  that  the  ratio  i?rs(440)/Rrs(550)  did 
not  vary  widely  (within  15%)  because  of  inelastic 
scattering  (see  Table  4).  Among  stations  without 
bottom  influence,  differences  in  the  ratio  were  within 
~  10%,  which  suggests  the  spectral-radiance  ratio  is 
effective  for  most  deep  waters  without  consideration 
of  CDOM  fluorescence  and  water  Raman  scattering. 
But  it  is  obvious  that  as  the  bottom  influence  in¬ 
creases  the  usefulness  of  the  power-law  algorithm 
decreases.  Also,  the  power-law  algorithm  cannot 
distinguish  between  the  absorption  of  CDOM  and 
that  of  pigments.  Note  also  that  Rrs“'(490)/iirs(490) 
values  as  low  as  0.77  were  determined  (not  explicitly 
shown),  which  suggests  that  great  care  must  be  taken 
when  one  is  interpreting  remote-sensing  curves  for 
the  intermediate  wavelength  at  shallow  coastal  sta¬ 
tions. 

At  the  optically  shallow  stations  (STOl,  ST02, 
ST03,  and  ST19),  the  model-derived  depths  were 
within  approximately  10%  of  the  chart  depths  with¬ 
out  consideration  of  any  tidal  influence  (typically 
<  0.5  m).  This  demonstrates  a  potential  to  use  this 
model  to  survey,  e.g.,  by  aircraft  overflights,  dramatic 
changes  in  shelf  bathymetry  that  can  occur  as  a  result 
of  major  storm. 

For  STOl,  ST02,  ST03,  ST08,  ST12,  and  ST14,  the 
general  agreement  between  the  modeled  and  mea¬ 
sured  Rj-s  values  are  very  good,  with  small  differences 
near  580  nm,  where  the  measured  Rj.^  >  modeled  R^s- 
Other  than  the  modeling  error,  there  are  at  least 
three  possible  reasons  for  this:  (1)  bottom-albedo 
uncertainty,  (2)  phycoerythrin  fluorescence,*'^  and  (3) 


Table  4.  Optical  Component  Contributions°  to  R„ 


Station 

R  “ 

-O-rs 

/■^rs 

Rr.f^/R 

rs 

R 

rs^/^TS 

ii,s(440)//? 

s(550) 

440  nm 

550  nm 

440  nm 

550  nm 

440  nm 

550  nm 

Measured  R^s 

Corrected  Rrs 

STOl 

MM 

0.02 

0.10 

0.87 

0.92 

0.91 

0.04 

0.12 

1.48 

1.73 

0.86 

0.81 

0.10 

0.12 

0.02 

0.04 

1.87 

2.00 

ST08 

0.98 

0.99 

0.03 

0.02 

0.49 

0.48 

STIO 

0.92 

0.93 

0.08 

0.07 

1.64 

1.60 

ST27 

0.90 

0.94 

0.08 

0.07 

1.34 

1.29 

ST12 

0.90 

0.97 

0.09 

0.03 

0.31 

0.27 

ST14 

0.97 

0.99 

0.03 

0.01 

0.32 

0.30 

ST19 

0.92 

0.92 

0.06 

0.07 

0.03 

0.01 

2.46 

2.36 

=  i?rs^  +  iJrs^  and  corrected  R^s  =  measured  R^s  —  Rrs^  —  RrJ’-  Blank  entries  indicate  that  no  data  are  available. 


water-absorption  coefficient  uncertainty.*^  A  spec¬ 
trally  constant  bottom  albedo  was  used  for  the  shal¬ 
low  stations.  Earlier  measurements  of  bottom  al¬ 
bedo  (Fig.  7)  did  display  some  spectral  dependence, 
but  these  types  of  changes  could  not  provide  the 
sharp  increase  and  then  decrease  with  wavelength  in 
i?rs  required  for  the  measured  and  modeled  R^s  curves 
to  converge.  Also,  there  was  no  bottom  contribution 
to  i?rs  at  ST08,  ST12,  or  ST14.  More  realistic  expla¬ 
nations  include  the  lack  of  a  term  for  phycoerythrin 
fluorescence  or  the  differences  between  the  water- 
absorption  coefficients  in  this  spectral  region  re¬ 
ported  by  Smith  and  Baker^®  and  Tam  and  Patel.^® 
Further  study  is  required  to  resolve  this  issue.  The 
differences  between  the  measured  and  modeled  R^s 
curves  near  685  nm  are  expected  because  no  term  is 
included  in  the  model  to  describe  the  chi  a  fluores¬ 
cence. 

Summary 

Contributions  to  the  water-leaving  radiance  spectra 
for  a  variety  of  waters  were  attributed  to  elastic 
shattering  by  water  molecules,  suspended  particles, 
and  bottom  reflectance,  and  to  inelastic  scattering  by 
water  Raman  and  CDOM  (or  gelbstoff)  fluorescence. 
Inelastic  scattering  by  pigments  was  not  considered. 
For  optically  deep  water,  remote-sensing  reflectance 


wavelength  (nm) 

Fig.  7.  Bottom-albedo  spectra  for  near-sbore  (dasbed  curve)  and 
offshore  (solid  curve)  sediments,  measured  for  samples  on  earlier 
cruises  and  retrieved  from  a  grab  sampler  or  by  divers. 


of  the  water-column  part  (elastic  scattering  only)  Rj-s'" 
was  simulated  as  follows: 


R 


W 

rs 


0.176 

—  4-  1 

d—Y] 

+  Og  +  Up 

V/J 

(23) 


where  5bm  is  known,  can  be  estimated,  and  A  and  Y 
are  spectral  constants.  averaged  approxi¬ 

mately  3.3,  and  Y  was  less  than  2.4  for  the  waters 
considered.  For  optically  shallow  waters,  the  expres¬ 
sion  for  bottom  reflectance 


RJ  *  0.17p  exp[-(1.5  +  {Dj})aH],  (24) 

works  well  for  the  shallow  waters  that  we  considered. 
Together,  the  water-column  term  and  the  bottom- 
reflectance  term  accounted  for  more  than  90%  of  the 
total  remote-sensing  reflectance. 

Close  agreement  between  modeled  and  measured 
Rrs  was  achieved  for  all  stations  when  all  scattering 
mechanisms  mentioned  above  (both  elastic  and  inelas¬ 
tic)  were  included.  The  ratio  a^(440)/ap(440)  cov¬ 
ered  a  range  from  ~  0.3  to  ~  3.0,  which  indicates  the 
wide  usefulness  of  the  model.  For  contributions 
other  than  those  from  the  water  column,  as  much  as 
23%  of  i?rs(490)  is  attributable  to  water  Raman, 
CDOM  fluorescence,  and  bottom  reflectance  for  an 
optically  shallow  (25  m)  station.  For  the  power  law 
of  spectral  radiance  ratio,  most  error  comes  from 
reflected  bottom  radiance  for  coastal  waters.  For 
most  deep  waters,  the  power-law  algorithm  can  be 
used  without  the  correction  of  CDOM  fluorescence 
and  water  Raman  with  little  error  if  the  optical 
properties  covary  with  chlorophyll,  because  the  inelas¬ 
tic  effects  cover  the  whole  range  from  400-600  nm 
Avith  generally  less  than  10%  contributions  to  the 
water-leaving  radiance. 

The  Up  required  by  the  model  is  generally  within 
15%  of  the  measured  a^,  with  an  average  difference  of 
8.9%  (6.9%  when  ST12  is  excluded).  This  suggests  a 
method  to  remotely  measure  the  pigment-  and  gelb- 
stoff-absorption  coefficients,  although  derivation  of 
chi  a  concentration  will  depend  on  knowledge  of  the 
specific  absorption  coefficient  for  a  region.  Also  the 
model-derived  bottom  depths  for  optically  shallow 
waters  are  within  10%  of  the  chart  depths,  which 
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suggests  its  possible  use  to  remotely  measure  bottom 
depth  for  the  shelf  waters. 

Appendix  A:  Simple  Estimation  of  Qm{J, 

If  we  follow  the  method  described  by  Jerlov,^®  the 
subsurface  upwelling  irradiance  X.)  caused 

by  subsurface  sunlight  X)  can  be  obtained 

for  single  scattering  and  quasi-single  scattering^®  by 


Appendix  B:  Remote-Sensing  Reflectance  for  CDOM 
Fluorescence  and  Water  Raman 

For  z  positive  downward  from  the  surface,  0  the 
zenith  angle,  and  cj)  the  azimuthal  angle,  with  the 
consideration  of  isotropic  Pie,  to  the  first  order,  the 
inelastic  radiance  L„  ie  in  the  direction  0  and  the 
upwelling  irradiance  ie  at  depth  z  resulting  from 
the  depth  interval  dz  are  simplified  to 


X)  =  x) 


sec(y)exp[-cz  sec(y)] 
c 


=  1 


dz 

cos(0) 


(Bl) 


f"  p(«,X)sin(e) 

^  Jo  L/2 


(Al) 


where  a  is  the  scattering  angle,  0  is  the  zenith  angle, 
(j)  is  the  azimuth  angle,  and 


cos(a)  =  -cos(0)cos(7)  -t-  sin(0)sin(_;)cos(4)). 

The  upwelling  radiance  from  the  nadir  resulting  from 
this  situation  is 


plrr  —j,  X) 

X  exp[-cz  sec(y)].  (A2) 

Recalling  Austin’s  definitions^  for  the  Q  factor  E^jEu, 
we  get 


r  p(a,  X) 

P(Tr  -  J,  X) 
cos(J)  +  1 

(A3) 

Because  p(a,  \)  for  water  molecules  is  given  by  Mo¬ 
rels  and  its  angular  distribution  is  considered  to  be 
wavelength  independent  for  the  visible  region,  the 
above  equation  can  be  simplified  to 


~  5-92  -  3.05  cos(y).  (A4) 


As  backscattered  skylight  also  contributes  to  the 
upwelling-radiance  field,  its  influence  on  the  actual 
Qm(j,  X)  needs  to  be  considered.  Defining  the  ratio 
between  the  subsurface,  downwelling  sky  irradiance 
and  solar  irradiance  to  be  y(X)  and  assuming  that  the 
Q  factor  that  is  due  to  skylight  is  ~  3.14,  we  have 


Qm(j> 


1  +  7(\) 


1  +  y(X) 


3.14 


(A5) 


and 


d^u,ie(2,  0,  A.)cos(0)sin(0)d0 

tt/2 

=  2t7  j  Pie(X;„  X)EJz,  XjdX^dz,  (B2) 

where  EJz,  Xj  is  the  scalar  irradiance  at  depth  z,  and 
Eo(z,  Xj  =  Eod(z,  X^)  +  Egjz,  X^)  =  D(i(l  -I-  2R(X^)) 
\c)-  Consider  that  R{X^)  is  small  (<0.05)  and 
DujDdi  ~  2)^°  is  independent  of  depth,  the  subsurface 
irradiance  resulting  from  the  inelastic  scattering  for  a 
deep  water  column  is 


EudO  ,  X)  «  2tt 


I 

vXv 


^)Ed{0  ,  Xj 

2k(X)  +  k(XJ 


(B3) 


If  one  defines  Qie  as  the  Q  factor  for  the  inelastic¬ 
scattering  field,  then  the  subsurface  upwelling  radi¬ 
ance  caused  by  inelastic  scattering  is 


•^'u,ie(9  ,  X) 


211 

Qi 


Pie(^x>  ,K)  ,, 

2k(X)  +  k(XJ 


(B4) 


The  inelastic  total  scattering  coefficient  v|j(Xj,,  X) 
(m*  V nm)  is  defined  as 


=  Pie(«,  X^,  X)da).  (B5) 

‘'4ir 

Because  Pie(a,  X^,,  X)  is  considered  isotropic,  then 

X)  =  4iipJX,„  X).  (B6) 

According  to  the  definition  of  remote-sensing  reflec- 
tsmce  and  with  Eqs.  (B4)  and  (B6),  we  have 


■^rs,ie(k) 


If 

2Qie  ,  X 


il>(x^,  X)Erf(0  ,  XJ 
[2k(X)  +  k(XJ]£,(0-,  X) 


dX,.  (B7) 


So,  for  molecular  scattering  and  when  j  and  7(X)  are  For  CDOM  fluorescence,  if  'Ti(X,r)  is  defined  as  the 
known  for  any  station,  Qm{j,  X)  can  be  estimated.  quantum  efficiency  for  the  emission  band  excited  by 
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then®’^® 


'n(^^x)  = 


f  X  k) 

Ik  a,{K) 


(B8) 


(|;(\j.,  k)  can  be  characterized  by  a  log-normal  curve,®® 
so 


vl,(X„  \)  = 


\A 


exp 


(B9) 


in  which 


A  = 


(BIO) 


where  iq(Xj:),  Xqj  s,  and  a  may  vary  with  the  type  of 
CDOM  and  K. 

In  general,  bb  ^  a  for  most  oceanic  waters,®^  so  k  is 
close  to  a,  and  based  on  the  calculation  for  chi  a 
fluorescence  made  by  Gordon^®  the  Qie  factor  for 
inelastic  scattering  is  ~3.7.  Then  combining  Eqs. 
(B7)  and  (B9),  the  remote-sensing  reflectance  caused 
by  CDOM  fluorescence  can  be  reduced  to 


Rj{k)  «  0.072 


K  cig(KWd{0  .  K) 

k  [2a(\)  -t  a(X,)]£rf(0-,  k) 


exp  -s  In 


d\^. 


(BID 


Unlike  broadband  (~100  nm)  CDOM  fluorescence, 
the  water  Raman  emission  has  a  half  bandwidth  of 
approximately  20  nm.®®  If  this  bandwidth,  is  omit¬ 
ted,  i.e.,  a  narrow  Raman  emission  is  assumed,  the 
inelastic-scattering  coefficient  i[)(X;„  k)  for  water  Ra¬ 
man  can  be  related  to  the  Raman  scattering  coeffi¬ 
cient  as 


ilj(\^,  \)d\^  =  buiK), 


and  from  Eq.  (B7),  with  k  ~  a,  the  remote-sensing 
reflectance  for  water  Raman  is 


Rj{k)  =  0.072 


bR(K)Ej(0  ,  K) 
[2a(\)  -t-  a(X;,)]E<i(0“,  X) 


(B12) 
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